一 
ee 


[i 


aa 1 。 


第 一 部 分 “ 变 分 法 和 有 限 元 法 


第 一 章 引 论 


先 谈 一 点 变 分 法 的 历史 。 

变 分 学 历史 上 第 一 个 重要 问题 是 生 辆 提出 的 。 他 研究 了 所 谓 
“水 桶 问题 ”， 大 概 是 这 样 。 一 水 桶 以 赔 定 角速度 绕 对 称 轴 族 
转 ， 间 桶 中 水 面 旦 什么 形状 时 所 受 姐 力 最 小 ? 他 把 这 个 河 题 归结 
为 : 求 函 数 y= f(z)， 使 积分 


(2 FCF CR) 
QCFKzI=-| dr 


取景 小 值 。 这 问题 已 经 是 典型 的 变 分 问题 ， 但 牛顿 当年 用 来 解决 
问题 和 的 方法 还 不 是 变 分 学 中 的 基本 方法 。 

导致 变 分 法 建立 的 著名 问题 是 瑞士 数学 家 约翰 、 贝 努 里 在 
1696 年 提出 的 ， 即 “最 速 降 线 问题 ”。 此 问题 发 表 于 当年 6 月 号 
的 《教师 学 报 》 上 。 问 题 一 提出 来 ， 立 即 吸 引 了 当时 世界 上 许多 
最 卓越 的 数学 家 的 注意 。 在 该 杂志 1697 年 5 月 号 上 刊登 了 牛顿 ， 
沫 布 尼 兹 、 洲 出 达 和 贝 努 里 兄弟 吕 的 解法 。 他 们 殊途同归 ， 用 各 
不 相同 的 方法 得 到 同一 解答 。 答 案 是 有 趣 的 ， 原 来 “最 速 降 线 2” 
居然 是 摆 线 ! 这 些 解 题 方 法 中 昔 含 着 天 才 的 思想 ， 就 像 铀 矿 中 储 
焉 着 铀 -235 一 样 。 但 是 只 靠 这 些 方法 ， 变 分 法 是 创立 不 起 来 的 。 
正如 只 依靠 古代 的 极限 思想 无 法 建立 微分 法 一 样 。 那 么 ,. 谁 来 清 
除 杂 质 ， 炼 出 真 金 呢 ?历史 安排 了 大 数学 家 尤 拉 ! 他 是 创造 各 种 
数学 方法 前 大 师 ， 他 创立 的 方法 几乎 总 带 有 明基 性 和 一 般 许 。 他 
在 1734 年 解决 了 更 广泛 的 最 速 降 线 问题 ， 但 他 对 自己 当时 用 的 方 
法 还 不 满意 。 他 开始 寻找 解决 这 类 问题 的 一 般 方 法 ， 这 种 方法 终 


由 “ 风 努 里 兄弟 ” 指 James Besraoulli 和 他 的 弟弟 Joha Bernoulli。 
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于 被 他 找到 了 ， 变 分 法 也 就 建立 了 。 尤 拉 在 1736 年 的 论文 中 指 
出 :要 使 积分 


Qcry= {Preyr) ,y's) dz (1) 


取 极 值 ， 肾 数 y(x) 必 须 满足 

oF -3 )=0 

Ov dr\ ody’ 
方程 《2 ) 就 是 著名 的 “ 尤 拉 方程 ”。 不 过 ， 许 多 重大 数学 成 果 
的 第 一 个 证 明 往 往 失 之 繁 融 ， 尤 拉 这 次 也 不 例外 。 后 来 法 国 大 数 
学 家 拉 格 兰 日 改进 了 尤 拉 的 证 明 ， 使 之 完全 符合 数学 分 析 的 精神 
而 且 十 分 简洁 ， 他 在 1755 年 把 这 个 证 明 告 诉 了 尤 拉 。 至 此 ， 变 分 
学 作为 一 个 新 的 数学 分 支 算 是 形成 了 。. 

刚才 谈 到 的 是 变 分 法 古典 部 分 的 内 容 ， 它 所 研 性 的 主要 问题 

可 以 归结 为 ， 让 适当 的 函数 集合 内 选择 函数 y= f(z)， 使 积分 
《1 ) 取 极 值 。 解 决 这 一 问题 又 归结 为 解 尤 拉 方 程 《2》 。 看 起 
来 问题 似乎 并 不 复杂 ， 办 法 也 很 平常 ， 不 过 解 一 个 微分 方程 黑 
了 。 其实 不 然 ， 我 们 依靠 这 种 方法 ， 就 能 用 统一 的 数学 程序 来 解 
决 自然 界 和 其 他 方面 的 千差万别 的 问题 ， 就 能 用 奇妙 的 变 分 原理 
解释 无 数 的 自然 现象 。 


《2) 


参考 文 献 
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$1 泛 函 


在 这 里 不 介绍 泛 函 的 精确 定义 ， 因 为 那样 做 需要 比较 多 的 现 
找 数学 知识 ， 势 必 昂 大 圈子。 实际 上 ， 用 几 个 例子 由 可 以 说 明 
“ 泛 函 ”这 一 概念 的 实质 。 
例 1 设 阔 数 y= 斤 z) 在 闭 区 间 50 , 1 上 连续 ， 且 f(XY) 守 
0 ， 则 微分 
QLA=) fs)dz 


是 曲线 2 = f(x) 与 直线 y= 0 之 闻 的 面积 了 9。 显然 ， 对 任何 一 个 
在 [0 ，1 J 上 上 连 续 且 非 负 的 函数 y 二 了 (zx) 来 说 ，Q[ 记 有 唯一 确定 
的 值 与 它 对 应 。 比 方 说 ， 当 


J{z)=x 
则 
QA)I= | zur=l 
当 fc) 
则 Qrf(z= 人 = 有 
0 3 


现在 我 们 回忆 一 下 一 元 函数 的 定义 。 对 于 两 个 变量 和 2， 
衬 在 某 数 集 刀 中 取 值 ，2 在 基数 集 闵 中 取 值 。 如 果 任 取 忆 中 基 
值 ，y 都 有 唯一 确定 的 值 与 它 对 上 应， 那么 我 们 就 称 变量 x 为 自 变 
量 ， 变 量 2 为 自 变量 * 的 耳 数 。 我 们 常 说 两 数 5 的 值 由 自 变量 


TD 若 折 四 尺 了 和 坪 为 于 ， 则 QE 让 是 页 璀 线 y = 了 tx)， 有 曲线 # = 8， 直 线 
工 - 6 太 直 站 = 也 转 成 的 图 形 的 面积 。 


上 


的 值 唯一 确定 。 函 数 定 义 是 大 家 熟悉 的 ， 其 实质 是 “唯一 确定 ? 
四 个 字 。 再 返 观 例 1 ， 当 函数 扰 z) 给 定 了 ，QLCAz)] 的 值 就 " 唯 
一 确定 ”了 ， 这 意味 着 什么 ? 这 说 明 我 们 可 以 把 变量 QC 了 (7z)] 
看 成 是 “函数 (zx) 的 汕 数 ”。 这 QLI(z)] 也 记 为 985[f1]， 在 此 ， 
f 类 似 自 变 量 ，QQ 是 了 的 函数 。 

例 2 我 们 把 一 切 在 [cs，p] 工 连续 的 函数 ( 指 实 函数 ， 以 
下 辐 此 》 的 集合 记 为 CTLa，b)。 考 虑 积分 


MC#I=) flr)ds 


任 取 一 个 在 Ca ,29 上 连续 的 函数 近 z)，M 有 唯一 确定 的 值 与 它 对 
应 。 邮 可 视 为 了 的 函数 。 比 例 和 例 1 比较，f 的 容许 范围 广 些 。 

我 们 看 出 在 M 和 了 之 间 ， 有 一 种 国 数 关系 ， 但 这 是 推广 了 的 
函数 关系 ， 也 就 是 下 文中 要 讲 的 “省 国 ”。 原 来 所 谓 “ 证 函 ”， 
不 过 是 更 广泛 意义 下 的 函数 关系 轩 了 。 

例 5 把 大 z) 在 [a， 细 上 的 零点 个 数 宁 记 作 六 < 门 ， 则 任 取 
一 全 在 [ae，pb] 上 有 定义 的 实 函 数 态 z), 立 有 确定 的 值 与 它 对 应 。 
如 取 La，Bb] 为 [0 ，2z]， 则 

N[Leosr]= 2 
Csinxj= 3 
NCz?~—1)= 1 


从 以 上 志 个 例子 中 ， 我 们 可 以 概括 出 如 下 的 重要 概念 ， 把 具 
备 某 种 性 质 的 函数 的 集合 记 作 厂 。 对 于 集合 了 中 的 任何 亢 数 六 = ) 
〈 即 对 任何 了 E D) ， 变 量 包 都 有 唯一 确定 的 值 与 它 对 应 ， 那 么 
变量 CO 叫做 依赖 于 项 数 放 zz) 的 泛 画 ， 记 为 
Q=Qrfz)] 或 Q=Qr 门 @ 
这 是 泛 函 概念 的 一 个 粗浅 的 介绍 ， 对 于 这 本 书 来 说 ， 这 种 说 明 已 
经 够 用 了 。 


中 我 们 把 m 重 零点 算 可 个 检点 。 加 (x) = (zz - 3) 在 [1，4J 上 于 有 两 个 零点 。 
加 ”QLF (TD 也 记 作 Qcy tx) 3， 在 此 y= y(tz) = f(z)。 


§ 2 


1. 泛 肖 QEy(x)2 中 ,ylXx) 的 变 分 

对 于 泛 函 QLy(X)],y(x) 是 集合 局 中 任何 元 素 。 加 果 ytx) 由 
yo(T) 变 成 y,(z)， 则 YCX) 一 yo《2) 思 做 y(z) 在 yo(x) 上 的 变 分 
(yo?) 及 yi.(7) 均 属于 D) ， 记 作 


bs 


分 


Ov=y(rT)— yr) (3) 
我 们 以 后 常用 
6y=y(7)— yz) (4) 
这 是 指 在 y(Cz) 上 的 变 分 。 
2， 连续 泛 画 


对 于 活 冰 QCy(x)J 而 言 ， 如 果 当 yl) 的 变 分 6y 充分 小 时 ， 
局 的 改变 量 可 以 任意 小 ， 那 么 就 称 汉 区 Q[Lv(z)j 是 连续 的 。 例 
如 ， 对 于 泛 通 
Qey(z)=| yr)de 


而 言 ， 当 v(x)E CLa，6] (参看 本 章 § 1 的 例 2〉，QCy(z)] 
有 定义 。 对 任 给 的 :>0， 只 要 
£ 
i ly (zx) — yo 7)| 人 
(CCzZ) 及 zz) 均 属于 c[a，2]) 
则 
| QEyCz)I— QFyo #5) yx) yo 2) dz 


< 人 ty -yd < /de 


可 见 QCy(z)]=| v(x)dz 是 连续 泛 函 。 


3. 线性 泛 函 
如 果 泛 函 QI5Y(z 轨 与 y(X) 的 关系 是 线性 的 ， 出 就 是 说 ， 马 


要 B 二 


满足 以 下 条 件 

1) QIcy(z)] 二 CQR[Ly(z)] {c 为 任 总 常数 

2) QIy(x)+ a 
这 时 ， 称 QIvy(z)] 为 线性 泛 函 。 例 如 


Qrytz)={ v(xr)ar 
是 线性 泛 函 。 事 实 上 


QTcw (lr)+ cvs)I=) Cowlz) t+ ey rt) dr 


一 of vila)dr+ cef yal)dr= ciQKwtCzD)] 
+ C2 QL ya 2T)) 


但 是 
b 
Qcry(z)={ yd 


不 是 线性 泛 孙 ， 这 一 点 请 读者 自己 验证 。 

4. 泛 次 的 变 分 

在 上 文中 已 经 说 明 “ 泛 熙 ”是 函数 这 一 概念 的 排 广 。 现 在 我 
们 还 要 指 志 ， 所 请 “ 泛 孙 的 变 分 ”下 是 “函数 的 微分 ”这 一 概念 
的 推广 。 如 果 明 和 白 了 这 一 点 ， 您 就 懂得 了 “ 变 分 ”概念 的 实质 。 
为 了 把 问题 说 清楚 ， 我 们 有 必要 先 复 习 一 下 “微分 ”这 一 概念 。 

现在 我 提 一 个 看 来 非常 简单 的 问题 ， 什 么 是 沪 数 y 二 /(z) 的 
微分 ? 

如 果 您 马上 回答 说 ， 孙 数 的 微分 就 是 

dy=f{’(r}dr 

虽然 我 不 能 说 您 短 钳 了 ， 但 是 您 并 没有 涪 清 想 “ 答 分” 的 实 有 
要 是 您 在 谈 到 微分 定义 时 ， 只 能 说 出 上 击 这 和 铝 话 来 ， 那 么 她 我 直 
言 ， 您 并 没有 学 好 数学 分 析 。 

究竟 什么 是 函数 的 微分 昵 ? 我 们 知道 尊 数 关系 &= 大 z) 往 往 


时” 此 二 末 件 可 咱 一 个 式 子 来 表达 ， 即 
QQkclaitzy+casyatz)] aciaryitirItesQty:0zr)) 对 任何 常数 circ 成立。 


了 


是 非常 复杂 的 ， 因 此 函数 六 zz) 的 增 量 Ay = 上 zz+Az) 一 rz) 也 
往往 是 十 分 复杂 的 。 例 如 ， 当 
fw) 一 Sim 
这 和 时 函数 的 卉 量 是 
Ar 


Ay=sin(7 + AL) ~ sin = 2cos( 汪汪] < ) Sin 本 


涩 > 已 项 ， 用 此 式 二 Az 的 值 可 得 Avy 的 值 ， 俱 Azy 与 Az 之 间 的 函 
数 关 系 是 非 线性 的 。 在 数学 分 析 中 曾经 指出 ， 如 果 函 数 ! 一 三 z) 
全 给 定点 + 处 有 和 导数 f(r)， 这 针 


. Ay re 
lim A (rz) 


态 和 一 操作 
于 是 
A (zr) ta ( 当 A2-> 0 时 ， a-> 0) 
所 以 
Ay=/{’(x)ArtoaAr=dytaAr 《5) 


行 细 观察 此 式 ， 我 们 发 现 函 数 的 增 量 由 两 项 相好 而 得 。 第 一 项 是 
dy 一 1 (7)Ax， 当 I 国定 时 f(z) 是 常数 ， 盾 以 dy 是 与 Ax 成 比 
例 的。 通常 说 dy 是 At 的 线性 通 数 。 第 二 项 是 aAx， 当 Ax-> 0 时 
a 0 ,所 以 
CA 0 (Ar—2=0) 
可 见 c.Az 是 比 Az 高 险 的 无 穷 小 蝶 ， 即 
CN OfA7) {《 当 Ar 一 0) 
所 以 
Ay= 7 (TIAT+OCAr) C6} 
可 见 。 dy 二 了 (xz)Az 是 明 数 增 量 的 主要 部 分 ,也 就 是 说 f(x)Ax 
也 是 Ay 的 线性 部 分 ， 又 是 Ay 的 主要 部 分 ， 即 所 衣 “ 线 性 主要 部 
分 ”。 我 们 把 范 数 增 量 的 线性 主要 部 分 (x)Azx 叫做 通 数 的 微 
分 。 这 样 讲 ， 才 说 明了 微分 概念 的 实质 。 
由 此 可 见 ， 当 |Azx| 充分 小 的 时 候 ， 我 们 就 可 用 微分 dy= 


8 +。 


7 了 《x)Ax 作 为 丧 量 Ay 的 近似 信人 了 。 再 来 看 y = sinz 这 个 函数 ， 就 
有 如 下 近似 式 
Ay=sin(z AX)— sinresf’(r)Ar=cosrz Ar 
用 cosxAz 作 为 Ay 的 近似 值 ， 这 就 方便 多 了 。 可 以 说 ， 用 线性 
关系 来 逼近 非 线性 关系 ， 这 是 微 积分 学 的 基本 思想 1 
有 了 了 以 上 的 准备 ， 我 们 就 可 以 讲 诈 函 的 变 分 了 。 先 看 一 例 ， 


对 于 泛 酋 QCy(z)] 一 | .yw(z)dz， 它 的 增 量 可 表 为 
AQ= QIy.(7)I— QLy(r)) = QCy(r)+ oy — QLy(7)) 
=|. [Cyr) + OY dx -| y (rd 


=) Cy) taytz)6y+ (09) dz -| ya)dz 


-| 2y(x) Oydx + Cy):dz 


在 此 y=y(7) ~ y(r) 
可 见 ， 此 泛 函 名 的 增 量 AQ 由 两 项 相 加 而 得 。 将 第 一 项 记 为 


zx(zysyaz=TCwGz)],6o 


当场 数 VCz) 国 定时 ， 了 fy(Cz),5%] 是 关于 8 的 线 竹 泛 函 。 这 是 
入 为 对 任何 常数 C 而 言 ， 有 


Try(z),C5og= 人 22(z)Cisdz 


一 C| .avsyaz =CTry(r),6y) 


TEy(z), BY.+ 0y] = 人 2y(2)(6Y, + 6ys)dr 


| 2yCz)6y dr+ | 2g(z)6yadz 
=7LY(r), SY) + TLy(z), dy 
我 们 再 来 观察 第 二 项 | (5y)*dz。 在 此 5y 一 y.(z) 一 y(z)， 其 中 


人。 


y(T) 是 已 给 定 的 测 数 ，y.(x) 是 任 到 的 函数 ，y(z) 及 y1(z) 均 属 


于 Cfa,6b] 
若 max iy(z)—-y(z)| =max lSy| 一 一 0 
Ce Bs 
由 5 
jevyaz| < max (3y):(b—a) 
5 [| 
可 知 
人 szadz 
ee 一 一 > 小 
max |G2| 


这 就 是 说 ， 当 max 16y| 一 0 时 ，| (58y)*dz 是 比 max 6y| 高 价 
的 无 穷 小 量 ， 不 妨 记 为 
| swyaz= oSy) 
于 是 
AQ=7Try(r),5y]+0(6y) 
这 个 公式 与 《6 ) 式 何等 相似 ! 它 指出 一 个 事实 , 即 泛 函 


Q ra(z) 一 全 we(z)az 的 增 时 可 分 解 为 两 部 分 ， 第 一 部 分 是 


4G2 的 线 忻 证 图 ， 第 二 部 分 是 比 62 更 高 阶 的 无 穿 小 量 。 这 种 情 形 
与 函数 增 基 的 状况 多 么 相似 ! 知 司 的 读者 立刻 会 想到 ，“ 这 第 一 
部 分 !Zfy(z), 7] 是 不 是 泛 函 Q 的 变 分 呢 ? 正 是 这 样 ! 现在 我 们 
就 来 狐 述 “函数 的 变 分 ”的 定义 。 

定 交 ”对 于 泛 光 QCy(x)]， 给 yl(r) 以 增 量 5y ( 即 yCz) 的 变 
分 ， 参 看 本 章 8 2 中 的 第 (4) 式 ) ， 则 泛 函 @ 有 增 量 AQ=Q 
[ftz)+ 6y] 一 QCy(XY)J。 如 果 人 QQ 可 表 为 

AQ=7TIy(7),6yI+ Bly(r), dy] (7) 
在 此 ，7 [sg(z), 0 对 62 而 言 〈 当 yz) 给 定 ) 是 线性 泛 函 面 
CE2tz) .62 -> 0 (sy(z) 国 定 ，max |8y| 一 >0) 


max lul 


那么 ，7T[Ly(z),6yJ 称 为 泛 函 的 变 分 ， 记 作 5Q@。 可 见 泛 画 Q[y 


- 10 。 


《7 站 的 变 分 5Q 是 @ 的 增 景 的 “线性 主要 部 分 ”。 初 学 者 应 将 
“ 泛 函 的 变 分 ”与 “函数 的 微分 ”这 两 个 定义 作 仔 细 的 比较 ， 直 
到 您 懂得 了 后 者 是 前 者 的 推广 为 止 。 
再 举 一 例 。 将 在 [0 , 1] 上 有 定义 的 汞 数 的 集合 记 作 刀 了 。 任 
取 y(x) ED，y*(《0 ) 的 值 是 唯一 确定 的 。 所 以 QEy(x)]==y? (0) 
是 定义 让 吕 上 的 一 个 泛 函 。 当 y(z) 取 定 ， 求 此 泛 函 Q@ 的 变 分 68。 
显然 
AQIIz +t SyI— QLyCr)I=Cy(r)+ GUN- 
— yt0)=y(0) + 2y C0) :6 0) + Ov(0) 
~—y:(0)=2y(0)6y(0) + [Oyt0) 
在 此 ，6w(0)==6y 1).c= 二 Ey(7) 一 y(7T)Y,.o 
观察 AQ=2y(0)6y(0)-+ (5y(0))*， 可 以 看 出 
2y(0)6y(0) = 二 TCy(z),6y1 是 关于 6 的 线性 泛 汪 又 当 
max 164] 一 > 0 时 ， 


(Sy(0))” (max Si 了 
maxlSy| ~ max 16yl 


可 见 2y(0)678(0) 是 A 包 的 线性 主要 部 分 ， 也 就 是 Q 的 变 分 。 即 
6Q=2y(0)5y(0) 
注意 ， 对 于 闻 数 u(x) 一 y*(z)， 它 的 微分 是 
du=2y(x)dzx 
此 微分 式 和 上 面 的 变 分 式 是 多 么 相似 ! 这 使 我 们 想到 泛 函 的 变 分 
与 范 数 的 导数 也 许 在 放 算 方面 有 联系 。 事 情 果然 如 此 。 现 在 来 介 
绍 一 条 常用 的 定理 。 它 药 大 意 是 ， 当 泛 栈 GUIy(z7] 的 变 分 存在 ， 
这 时 我 们 固定 y(x*) 和 6y=y1(7) 一 y(z)， 考 虑 
Qfry(r)}+ady] 
显然 ， 它 是 变量 4 的 甬 数 ， 可 记 为 
QTy(r) tady]= (a) 
不 但 如 此 ， 而 和 且 隙 数 g(a) 在 4 = 0 处 的 导数 ww (0) 就 等 于 泛 郴 
QLy(7X)] 的 变 分 。 现 在 来 证 明 此 定理 。 
定理 ”如 果 泛 六 QEfy(x)] 的 变 分 6Q= 二 7 了 [y(xX),6y] 存在 ， 
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那么 此 变 分 等 于 函数 p(c)=QIy(zrz)+c6y) 的 导 函 数 在 = 0 处 
的 值 (y(z) 及 56y 均 国定 》。 
证 明 ， 当 y(z) 及 6y 均 固定 时 ， 泛 聊 Q[y(x)4a6y]=$(a) 


是 4 的 函数 ， 且 
2 一休 人 
= QFy(r) -Fady ~ QEy(z)] 
一 > 中 代 


因为 汉 范 Qiy(x)] 的 变 分 存在 ， 所 以 
QLy(r) + oy — QIy(r)] 
一 了 [oz Go 二 OLD GD 
于 是 
QLy(r) +adyl— QUy(r)) 
Ty(r),oy + piy(r),ad6y] 
当 y{z) 固 定 ，7TLy(z),6y 对 6y 而 言 是 线性 府 函 。 所 以 
TIy(z),ady)=aT Ly(z), dy) 


9 lim Po) — 0) 
这 


这 一 了 1 

= lim Laz),2983++ PCy(7), 08] 
一 四 人 人 

fi A 


=T[y(x), Syi+ lim EY 


但 是 ， 由 公式 《7》 可知 


人 {7 ,OH. y ! 
1 C 当 y(z) 固 定 ，max 16y| 一 >0) 
lim PLy(x) ,05y) in PLYy(x), ady) .max lady 
“0 2 ao ma |Gyl C 


Su max lzsy 
Bw | le | 


| = max |6y| 有 界 ， 有 . 


2 4« 

“p00)=TIy(r), Sy) (8) 

5. 泛 杂 的 极 值 

设 泛 函 Q[Ev(x)] 在 曲线 y(xz) 的 集合 上 有 定义 ， 曲 线 y(z) 沟 
位 于 平面 上 的 有 界 区 域内 【不妨 认 时 
为 8B 是 和 矩形) ， 如 图 1 所 示 。 

如 果 泛 前 RCy(x)] 在 任何 与 y。 
《z) 充 分 接近 的 曲线 y(z) 上 的 值 均 不 
天 于 CCzyo(Cz)]， 亦 即 泛 遂 的 增 量 
AQ=QIy(7)I ~ QIy(r)IED 
时 《y(z) 及 yol7X) 这 些 曲 线 均 位 于 有 
界 区 域 内) ， 就 说 泛 函 CE 
在 曲线 yu(z) 上 取 强 极 大 和 值 .关于 泛 函 
QCytz)] 在 曲线 yo(x) 上 取 强 极 小 
值 ， 有 类 似 定义 。 

为 了 后 面 的 需要 ， 我 们 还 得 把 “曲线 yo(z) 的 邻 域 ? 这 一 概 
念 介绍 一 下 。 我 们 把 适合 不 等 式 

ly(x)— yx)| Se 

的 一 切 幅 线 y(Cz) 的 集合 叫做 曲线 wz) 的 一 邻 域 电 。 在 此 我 们 还 
加 上 曲 钱 ylz) 及 yolz) 均 连续 且 均 位 于 有 界 区 域 BB 内 的 要 求 。 这 是 
《点 的 分 域 ”这 一 概念 的 推广 。 于 是 泛 函 极 值 定义 可 令 述 如 下 ， 

若 在 曲线 va(z) 的 某 一 个 e- 邻 域 中 任 取 一 曲线 VCz)， 总 有 

AQ=QEy(r)I— QLy(z)IE0 

则 称 泛 函 QULY(zJ] 在 晶 线 yo(z) 上 取 强 极 大 值 。 类 此 ， 有 泛 函 包 
Cy(z)) 在 曲线 yo(z) 上 取 强 极 小 值 的 定义 。 

不 过 ， 在 变 分 学 中 用 得 更 多 的 是 另 一 种 曲线 邻 域 概念 和 以 它 
为 基础 的 另 一 种 泛 函 极 值 定义 。 当 曲线 y(z) 属 于 曲线 yo(z) 的 这 
种 邻 域 时 ， 不 但 y(z) 和 yo(z) 的 值 很 接近 ， 而 且 它 们 的 导数 值 也 
很 接近 。 具 体 说 来 是 这 样 ， 首 先 设 yo(z) 是 包含 在 号 域 中 的 C: 类 
曲线 名 ， 凡 满足 条 件 


中 又 称 为 Yo《Z}) 的 堆 阶 5 一 邻 碱 ， 见 下 文 。 
回 当 yu(z) 有 1 阶 过 续 导数 ， 称 y,(z) 为 C1 类 昌 并 。 


1 
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1) ly(x)~—y(r) ee 
2) ly’(z)—yi(r) <e 
3) y(z) 也 是 8 域内 的 C' 类 曲线 
的 曲线 的 集合 ， 叫 做 曲线 yo(x) 的 1 阶 e- 邻 域 。 
类 此 ， 可 定义 曲线 yolx) 的 天 阶 s- 邻 域 。 
而 byo(z) 的 0 阶 8- 邻 域 好 上 文中 首先 说 到 的 那 种 邻 域 。 
对 泛 函 Q[y(z)] 来 说 ， 如 果 Q[yolz)] 是 QLyCx)] 在 yo 人 (x) 
的 菜 个 1 阶 - 邻 域 中 的 最 大 信 ， 则 称 泛 汪 QLy(z)] 在 yo(z) 取 弹 
极 大 值 。 关 于 泛 函 QLv(x)] 在 yoltzx) 取 弱 极 小 值 有 类 似 定 义 。 
弱 极 大 值 与 弱 极 小 值 合 称 “ 弱 极 值 ”。 在 下 文中 我 们 由 讨论 
泛 函 的 弱 极 值 ， 所 以 我 们 以 后 就 把 “ 弱 极 值 ” 简 称 为 “ 极 值 ”， 
这 不 至 于 和 前 面 那 种 强 极 值 定义 发 生 混 清 。 
说 了 这 人 么 多 的 定义 ,应 当 让 读者 轻松 一 下 了 。 还 是 举 个 例子 。 
吧 。 
例 对 于 泛 函 Q[y(z)] 一 人 TPRC3Jds 


在 满足 y(0) 二 0，y(1)= 1 且 有 一 阶 连 续 导 数 的 函数 类 中 ， 求 
使 QQ 取 极 小 值 的 曲线 yo(X)。. 

熟悉 数学 分 析 的 读者 一 举 即 知 QLy(z)]=| YiTVTCEJaz 
是 曲线 y(x) 玖 长 度 。 可 见 此 问题 实际 上 是 求 过 点 4(0,0) 和 8B 
{1,1) 的 和 类 得 线 中 最 短 的 一 条 ， 显 然 ye(z) 是 连接 4、 忆 两 点 
的 二 线段 。yo,(x) 的 表达 式 请 您 自己 写 出 来 吧 。 

也 许 您 还 不 满意 ， 还 要 求 不 用 这 种 几何 方法 来 解决 问题 。 这 
种 想法 很 好 ， 因 为 一 般 说 来 ， 随 便 取 一 个 泛 画 QLy(z)?,， 未 必 
有 明显 的 几何 意义 。 奢 么 ， 如 何 求 解 极 值 问题 呢 ? 办 法 在 以 下 各 
节 中 介绍 。 


$ 3 最 简单 泛 函 的 尤 拉 方 程 
尤 拉 方 程 是 变 分 法 中 的 基本 方程 。 在 本 节 中 我 们 对 形 如 


Qry(r)’ Ss F(z,y(r) yz)ds 


的 江东 把 它 呀 做 最 简单 尚 ) 来 导出 此 方程 。 
1. 预备 定理 
预备 定理 ”如果 国 数 ”= 上 六 zz) 在 ia,2] 上 连续， 又 


| rezo(z)az= 0 C9) 


对 任何 具有 如 下 性 质 的 函数 ?7(z) 成 立 。 这 些 性 质 是 ， 
1》7(+) 在 [a , 65] 上 有 连续 导数 ; 
2) nta)= 0=7( 5b); 
3) I2Czil 之 els 是 任意 给 定 的 正 
数 )。 
那么 ， 丫 数 1 (x) 在 [ a , 5] 上 恒 为 零 。 
证 明 车 有 zo€ta,b), 使 fi) 
0 ， 就 会 得 到 矛盾 。 论 证 方法 是 这 样 的 ， 


因为 (z) 在 re 点 连续 ， 所 以 有 正 数 3 存 ”了 条“ | 
在 ， 当 [Iz 一 zo 之 6 时 ，f/(z)> 0 成 立 。 ee 
现在 作 函 数 %(z)， 合 它 满足 图 2 


1) 当 XE€ (zro 一 ,ro+G) 时 ， pl(r)> 05 

2) 当 zE[a,zo 一 9] 或 zE[zo+Tg, pb] 时 ，&YCz) 一 0 

3) az) 在 [oa ,5 上 有 连续 导数 。 
w(xz) 的 图 象 如 图 2 所 示 。 观 此 图 可 以 想见 这 种 丽 数 %Cz) 是 可 以 
作出 的 ， 而且 是 很 多 的 〈 关 于 这 一 点 的 详细 叙述 见 下 一 音 ) 。 又 
取 7(Zz) 王 407z)， 显 然 儿 zz) 也 具有 上 述 三 条 性 质 且 对 任何 
之 0 而 言 ， 内 要 适当 选取 常数 4， 就 可 使 17(z)| 人。 对 于 这 
种 ?7(z)。 积 分 


人 Pr)n(z)dz= {fr) nr)dr> 0 


这 与 (9 ) 式 发 生 矛 盾 。 所 以 , 当 x El(a ,6b),f(x) 不 大 于 零 。 同 
法 可 以 证 明 ， 当 工 EC(a ,5)，f(zx) 不 小 于 零 ,可 见 当 rT E (a,6b) 
时 ，f(z) 三 0。 又 f(x) 在 z= 二 a ,bb 处 都 连续 ， 所 以 f(o)=0 
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二 f(b5)。 于 是 ， 当 x Ef[a,b】 

f(r)=0 
证 毕 。 

2. 预备 定理 的 精确 论证 

细心 的 读者 会 发 现 ， 上 面 的 论证 还 有 不 严格 之 处 。 那 就 是 函 
数 节 Cr) 的 作法 。 您 说 从 图 象 可 以 想见 这 种 %(z) 是 很 多 的 ， 我 就 
不 信 ， 请 您 拿 出 一 个 满足 那么 多 条 件 的 $8(z) 给 我 看 一 看 ! 这 种 说 
法 不 是 故意 抬杠 ， 而 是 表现 了 一 种 科学 精神 。 在 本 节 中 ， 我 们 来 
答复 这 个 问题 。 本 来 ， 我 们 只 要 作出 一 个 满足 上 一 节 中 所 说 的 三 
个 条 件 的 函数 %( 7 ) 就 行 了 。 一 般 的 变 分 法 教 本 差不多 都 是 这 样 
敌 的 。 但 是 ， 在 本 节 我 们 将 作出 满足 更 多 条 件 的 函数 &(z) 来 ! 这 
种 (x) 不 仅 满足 上 述 三 个 条 件 , 而 且 还 有 任意 阶 连 续 导数 。 它 在 
现代 数学 舞台 上 所 演 着 一 个 断 重 要 的 角色 。 如 果 您 今后 攻读 “ 广 
义 函 数 ”、“ 有 限 元 理论 ”这 些 正在 蓬勃 发 展 的 现代 理论 ， 您 经 
常会 碰见 一 个 函数 ， 它 的 表达 式 是 


ES 
yl rE(r -6x +6) 


10 当 z 各 (aro 一 9) 或 zE(ro+G,D) 


当代 的 数学 家 们 常常 称 此 画 数 为 “软化 子 ” (Mollifer)， 许 多 文 
献 中 经 常 引 用 它 的 有 关 性 质 。 此 函数 恰好 适合 我 们 的 需要 。 现 在 
就 来 证 明 它 满足 上 节 中 提 到 的 三 个 条 件 。 条 件 1 ) 和 条 件 2? 明 显 
成 立 。 只 须 验 证 条 件 3 ) 成 立 ， 即 证 明 此 %(z) 在 fo ,8] 上 有 连续 
导数 。 显 然 只 须 证 明 导 数 $' (7) 在 ro+6 及 r, 一 6 处 连续 就 行 了 。 
现在 只 验证 BW'(r) 在 +,，+6 好 连续 《在 +。 一 6 处 的 验证 方法 相 辣 )。 
于 是 ， 问 题 归 结 为 
94) 求 出 导数 2 (zeo+G) 


5 证明。 ram ,sg (z)=% (ro+6)。 


现在 来 求 (zs+G) 
当 z 之 X06+ 6 时 ，Y%(X) 三 0 


ae 16 。 
3 pT)— pr td) -0 
"Pott rz—(rot6) 
即 ¥(z) 在 点 zs 十 6 处 的 右 导 数 
$+ (ro+6) = 0 
再 求 左 导数 % 《xo6) 


1 
cr~r0 -37 


lm Pz)—BX+6) _ lm 2 
省 人 二 I— (ro+6) -mx0+5-0 zZ 一 Xo 一 G 
i 
0 os 1 (Tz—T0)—6? 
《G 十 G 
lim 2 (在 此 用 了 L “Hospital 法则 ) 


一 ,一 -co 1 > 1 -+ 1 
0 于》 


= lim (一 ‘to:+L 
(— 2)e V6:+ 1 


CO 
,lm Vert+i=6 
,met ,mm 
一 im r=0 


“lim (—2)e'tr Yor+i=0 
“, pilTot+6)= 0 
于 是 p(xrotod)=0 


还 要 证 明 
lim A 0 
ss 
先 求 右 极限 ,因为 当 X 之 zo+6 时 $8(z) 三 0 , 汤 以 当 xX 之 Xx,+8 时 ， 
‘(rz) 三 0。 于 是 
lim B(xz)= 0 
sgt+d+0 


再 求 左 极限 。 


ss 17 » 


当 z 之 z,++6 时 
1 7 1 
BL) = | Cz-z02-82 ] 三 必 ‘x-20 7-62 
一 2{7 工 一 0) 
[(z— rT) 6)? 
1 
lim (x) = lime ‘so -0 
srotrd-0 x0+s-0 
—2(7x— To) 


, ECz— xo) — 6 


=lim— 212B Ye:+ 工 一 0 


可 兄 %'(Cz) 在 点 ze+65 处 的 左右 极限 都 是 零 。 于 是 
(z) 一 0 一 区 (zo+G) 


扬 愉 pr(z) 在 点 zu 6 处 连续 。 则 理 可 知 %'(z) 在 rz 一 6 处 亦 连 
续 。 验 证 已 毕 。 

我 们 已 证 明 %(z) 在 [ea,8J 上 有 一 阶 连 续 导 数 、 用 数学 界 常用 
的 话 来 说 就 是 : “yz) 在 Ca ,5]J 上 具有 一 阶 光滑 程度 ”"。 现 在 我 们 
还 要 捐 出 一 点 , 使 用 刚才 的 论证 方法 ,可 以 证 明 % “(7Y) 在 zo 十 6 
这 两 点 连续 ,而 gc (z) 在 区 间 [a,63 上 其 余 点 显然 连续 ,于 是 函 
数 g《z) 在 Ca ,5] 上 有 任意 阶 连续 导 数 。 我 们 说 这 种 列 数 在 [a ,5] 
上 有 “无 窃 阶 ”的 光滑 程度 ， 这 是 把 #8t7) 称 为 “软化 子 ” 的 原 
因 之 一 (更 重要 的 原因 这 里 不 讲 了 ) 。 

我 们 现在 不 但 填补 了 预备 定理 的 证 明 中 的 漏洞， 而 且 由 于 
2%(z)7 有 任意 阶 连 续 导 数 ， 所 以 我 们 实际 上 得 到 了 比 原来 更 强 的 
结果 。 说 清楚 些 ， 我 们 得 到 了 : 

推广 的 预备 定理 ”如 果 序 数 y 二 f(x) 在 [a,6] 上 连续 ， 又 
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对 任何 具有 如 下 性 质 的 函数 7(z) 成 立 ， 这 些 性 质 是 ， 

1》7(z) 在 Ca,56J 上 有 + 阶 连续 导数 《+ 为 任何 给 定 的 非 负 
整数 ) 。@ 

2》7?(a)=0=TCD) 

3)》 |7(z) <e (为 任意 给 定 的 整数 ) 
那么 ， 函 数 ! 二 f(x) 看 fa,5J 上 恒 为 零 。 

请 注意 ， 在 这 里 若 7 二 1， 就 是 原来 的 预备 定理 。 

3. 最 简单 泛 函 的 尤 拉 方程 

现在 转 入 正题 。 我 们 来 研究 形 如 


QcyCz= 人 PTr ylr),y’(r)dr 


的 泛 函 ， 看 看 它 在 什么 样 的 曲线 y(z) 上 取 极 值 。 
说 得 更 清楚 一 些 就 是 ， 在 具有 下 询 性 质 的 函数 y(r) 的 集合 刀 中 
求 函 数 y,(z)， 人 和 使 泛 函 QEy(x)J 在 yot7) 上 取 极 值 。 这 些 性 质 
是 ， 

1) y(r)ECra, b)®, 

2) vy(a)==yo，y(B5)=y 《yo,y 1 是 任意 给 定 的 数 》 

3》 y(z) 曲 线 位 于 平 画 上 的 有 界 区 域 有 内 。 

读者 还 记得 微分 学 中 关于 函数 极 值 的 费 马 定 理 吧 。 此 定理 指 
出 若 函 数 y= 六 7z) 在 [ca,5] 的 某 内 点 zu 处 取 棚 值 ， 目 关 (zo) 存 
在 ， 则 f(rxo)=0。 这 是 函数 在 内 点 存在 极 值 的 必要 条 件 《 当 时 数 
存在 ) 。 这 种 使 导数 为 零 的 点 叫做 函数 的 驻 点 。 

也 许 您 读 到 这 里 时 会 提出 一 个 疝 题 ，“ 既 然 泛 函 是 函数 概念 
的 推广 ， 泛 了 芳 的 变 分 是 函数 的 微分 的 推广 ， 而 泛 函 极 值 又 是 函数 
极 信 概 念 的 推广 ， 那 么 ， 费 马 定 理 为 什么 不 可 以 推广 呢 ”? 如 果 


全 ”我们 都 定 了 7{z) 的 零 阶 导 男 数 即 了 (z) 本身。 
图 ”我们 治 在 Ce,p3 上 有 1IT 阶 连续 导数 的 丽 数 的 集 含 ， 记 作 CICa,p1。Ceraypb] 
的 含义 类 此 。 当 k= 0 CCayb] 即 Crab]。 
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您 这 样 想 的 话 ， 我 敢 说 您 在 数学 上 有 相当 强 的 想象 力 ! 但 是 我 还 
要 反问 一 句 ，“ 您 打算 怎样 推广 费 马 定理 呢 ? ? 

“既然 函数 六 z) 在 内 点 zo 处 取 极 值 , 则 在 该 点 的 导数 (如 果 存 
在 的 话 ) 为 零 ， 那 么 当 泛 范 QCU(CzDI] 在 曲线 goCz) 上 取 极 值 ， 
QE2(z) 在 go(z) 的 变 分 《如 果 存 在 》 也 应 该 是 堆 了 2” 。 

您 真是 “不 鸣 则 已 ， 一 鸣 尺 人 ”1 由 这 个 想法 就 可 以 导出 那 
著名 的 尤 拉 方 程 来 。 下 面 我 们 来 详细 地 谈 这 个 问题 。 

定理 设 F(x，#，y ') 是 三 个 变量 的 连续 函数 ， 且 当 点 
(x ,y) 在 平面 上 的 某 个 有 界 域 旦 内 ,而 y' 取 任何 值 时 ,了 (x,y ,y ) 
及 其 直到 二 阶 的 偏 导数 〈 指 对 变量 z,Vy 及 2 “的 侦 导 数 ) 均 连 续 。 
若 在 满足 

1) y(x)EC!ra,b] 

2) ya)=yo, yb) = 

3) yl(7) 曲 线 位 于 六 耐 上 的 有 界 区 域 如 内 的 静 数 集 合 中 ， 
泛 函 QC[y(z) 在 某 一 条 确定 的 曲线 y(tx)〉 上 取 极 信 ， 且 此 曲线 
y(Z) 在 [a ,68]J 有 二 阶 连续 导数 ， 那 么 函数 y(7x) 满 足 微分 方程 

FP, ,= 0 
(在 此 ， 毛 谓 QAro(z) 在 y(Cz) 取 极 值 是 指 取 弱 极 值 ， 也 就 是 说 
QCy(z)J 在 y{x) 的 某 个 1 阶 2- 邻 成 中 取 最 太 值 或 展 小 值 。 参 看 
第 二 章 82 中 第 5 节 ? 。 

证 明 是 然 泛 丽 QLY(Cz) 在 曲线 zz) 上 取 极 值 , 所 以 当 yKCz) 
是 y(Cx) 的 某 个 1 阶 e- 邻 域 中 的 曲线 ， 则 QfYi(zZ)]J 空 (或 过 ) 
QEy(z)]。 所 谓 “yi1(z) 在 v(x) 的 1 阶 e- 邻 域 中 ”， 指 y,( 工 ) 满 
足 

(yi (7) ~ yr) Se 
lyi(r)—y’ (rz) Se 《 当 ELa,b)) 
yi(T)E CTa,b] 
现在 ， 任 取 一 个 函数 7(x)， 使 
1) 7(r)EC!'Ca,6) 


se 20 。 


2) 7(a) 一 0 一 7?(P) 
于 是 ， 当 el 充分 小 时 ， 曲 线 yiCz)= 一 2CZ)+aTCZ) 在 yz) 的 1 
阶 s- 邻 域 之 内 。 
事实 上 ， 当 lc| 充分 小 时 ， 
ly CT) 一 世人 | Se 
yi(xz)—y'(2) = lan’ (rx) ee 
现在 ，v(z) 及 ?Kxz) 均 已 给 定 ， 
于 是 泛 请 QC[y(z)+47(x)] 显 然 是 变量 4 的 沙 数 ， 可 记 为 
QLy(r) tan(r) =yp(0) 
当 a 二 0 时 
$60) = QLy(r)] 
取 极 值 。 这 就 是 说 ， 函 数 gta) 在 = 0 处 取 极 值 。 册 费 马 定理 可 
知 ， 若 $8'(0 ) 存 在 ， 则 
#$’(0)=07 
但 是 


pC0) = Ca) -| 地 | Rs Ca 
一 We a Pr HE TI)+tan(z) ,y(t) 
tan’(z)Jdz} 
0 RY 
= (Ff ,入 (ytaD+tF +] 四 
=| CF n+ Fn ede 
= 到 d (Fy'yd 
= 人 (CF, Fr)dz (3(a)= 0 = 7(h)) 
所 以 gq 0) 存在。 于 是 p 上 (0)= 一 0， 即 


中 一般 的 夸 上 都 是 这 样 讲 的 。 坦 是 作者 认为 还 应 当 作 如 下 说 期 ， 斌 以 看 出 2=0 
是 函数 旬 人 2)? 的 定义 城 的 内 点 。 事 实 上 上， 若 & 比 4 略 小 一 些 ， 曲 线 Vz) +ea7Iz) 在 
1 的 一 阶 &- 邻 不 之 内 。 所 以 当 C 睹 小 于 0 时 ， 虽 (2 有 定 尺 。 间 理 ， 当 2 赂 大 于 0 
时 ，9 (ay 有 定义 。 所 以 aa 0 是 定义 城 的 内 点 。 


{crF,~ qd F'n(r)dr= OQ 
a dr 


EEF)= -Fry y= Ft Fy + Poy 


可 见 -FF ,存在 而 且 连 续 。 于 是 下,~- 拷 -中 ,' 是 连续 函数 ， 而 


7(z) 是 满足 条 牟 
1》#(x) 在 Co,b] 上 有 连续 导数 
2) 7(a)= 0 =7(6) 
的 企 和 何 函数 。 由 本 章 $ 3 中 的 预备 定理 可 知 


= 0 F,/=0 (10) 

这 就 是 著名 的 尤 拉 (Euler) 方 程 。 证 毕 。 

这 个 定理 中 有 一 个 条 件 显得 不 自然 ， 就 是 要 求 泛 函 CCYzI 
在 具有 二 阶 连 续 导数 的 曲线 y(z) 上 取 极 值 。 我 们 说 QCCz)] 在 
VCz)》 上 到 极 值 ， 是 指 到 能 极 值 。 也 就 是 和 葛 线 y《(z) 的 某 个 1 阶 
s 邻 域 中 的 C: 类 曲线 相 比 较 ， 泛 函 怠 在 y(2) 上 取 极 值 。 您 想 一 
下 ， 在 yx) 的 邻 域 中 的 那些 曲线 都 只 要 求 是 C! 类 的 ， 独 有 2V(z) 
本 身 却 要 求 是 C: 类 的 ， 这 就 显得 不 协调 ! 重要 的 数学 定理 往往 具 
有 内 在 的 和 谐 性 。 许 多 人 记得 大 诗人 杜甫 的 警句 ，“ 二 名 三 年 
得 ， 一 险 双 泪 流 ”! 但 人 们 往往 不 知道 数学 家 不 把 数学 定理 琢磨 
到 十 分 完美 的 程度 是 不 会 罢休 的 。 后 来 ， 终 于 证 明了 ， 对 于 泛 画 
Qtytz)]= | FCz,y(r),y'(z)Idz, 车 

1) F(z,y,y ) 对 X,Yy 及 y' 三 个 变量 的 二 阶 偏 导数 都 连续 ， 

2) Fr 0, 

3 ) QIry(z)] 在 C' 类 曲线 y(z) 上 取 极 值 ， 
则 此 kz) 有 二 院 连 续 导 数 。 这 个 定理 证 明 从 赂 。 但 是 由 些 可知， 


全 这 里 ?7(z) 记 汤 足 的 条 件 ， 比 3 的 预备 定理 中 的 条 件 少 一 条 。 对 7(z) 时 求 
的 条 件 越 少 ， 这 种 ?4z) 就 越 多 ， 预 备 定理 的 结论 更 能 成 立 。 
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只 要 广 ,'' 妆 0， 可 以 取消 定理 中 y(z) 是 C* 类 曲线 这 一 条 假设 ， 
这 就 自然 得 多 了 1 

我 们 把 满足 尤 拉 方 程 (10) 的 函数 叫做 “ 泛 函 Q 的 极 带 ”。 中 
“ 极 带 ”这 一 概念 与 分 析 中 的 “ 驻 点 ”相当 。 对 于 溺 数 p(x) 而 
言 ， 若 pzo)=0 , 则 称 zo 为 P(x) 的 驻 点 ,在 “ 驻 点 ”函数 p(x) 
不 一 定 取 极 值 ， 同 样 ， 在 “ 极 带 ”y(x) 上 泛 函 QI[y(7)] 也 不 一 定 
取 被 值 。 但 是 ， 在 一 定 条 件 下 《参看 本 节 的 定理 》 若 泛 函 


Qfy(z)]=| Cz,y,y')dz 在 曲线 yCz) 上 取 极 入， 则 y (7) 
是 极 带 。 
现在 来 举 一 些 例子 。 
例 1 求 泛 函 QEy(z)]= 人 (一 27y +y'*)dz 的 极 带 
解 。 “F(z,y yD=~27y ty 
“FF,=0 
PF,/:—=—2r+2y’ 


如 Fy’=—2+2y” 


dr 
于 是 ， 尤 拉 方 程 〈10) 即 
一 2 二 25 一 0 
所 以 极 带 为 


y(z)= 广 ?+Ciz+C， (Ci 及 C, 为 任意 常数 ) 
例 2 求 泛 函 Qty(z)] 一 人 (zy 二 2)dz 的 满足 y( 0 )=0， 


yl 1 )= 可 的 极 带 。 


解 : “P(ry,y = yy 
“= 


名 有 化 世上 把 满足 尤 拉 方 程 的 函数 称 为 泛 函 的 " 航 值 曲线"， 这 样 容易 使 妆 学 阁 
误 认 为 在 这 种 曲线 上 泛 防 一 定 取 极 值 。 所 以 在 本 书 中 采用 “ 极 带 ”一 词 ， 引 自 斯 米尔 
若 夫 所 著 《 丸 高 等 数学 教程 第 4 卷 1 分 册 站 译本 p210、 


d 
dz 


于 是 ， 尤 拉 方 程 为 
zz 一 202 一 0 


Fr = (2) =2y" 


它 的 通 解 是 
yz) = + Cr+tCs, 


-1 -1 
yt Cts 3 


.Cs 二 0，Ci 二 了 。 泛 函 的 极 带 为 


We 
二 -一 “十 一 ; 
YC7) 24 24" 


例 3 求 泛 画 QCy(x)]=| y*dz 的 极 带 
解 ; P(r yr = 


“y=0 
于 是 尤 拉 方程 为 
工 , 一 0 
但 
r=2y 
“Y= 0 


极 带 是 函数 % = 0 (显然 泛 函 QCy(z)=| yrdx 在 曲线 J 人 (x) 
= 0 上 取 最 小 值 )。 

这 个 例子 告诉 我 们 ， 著 矿 (z,y ,y 人) 中 不 含 y'， 则 龙 拉 方 和 
化 为 

F(xz,y)=0 

之 形 ， 这 已 经 不 是 微分 方程 了 (经 党 有 人 说 “ 尤 拉 方 程 是 微分 广 
程 ”， 这 人 句 话 是 有 毛病 的 ， 所 以 它 的 解 不 含 任意 常数 。 在 这 种 
情形 下 ， 尤 拉 方 程 的 解 往往 是 一 条 或 几 条 曲线 ， 因 此 不 能 再 要 求 
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它 满足 任意 给 定 的 边界 条 件 
ya)=C, 及 y(6)=C 

Ts 

现在 ， 我 们 再 来 研究 在 本 章 §2 中 所 提 到 的 泛 丙 Q(y(z)]= 
[wiry3a4r (还 要 求 y( 0)= 0 ,vy( 1)=1) 的 极 值 问 题 吧 。 
当时 借助 于 几何 直观 看 出 QiyCz)] 在 过 400,0) 和 B(1,1) 两 点 
的 直线 上 好 最 小 值 。 现 在 ， 我 们 用 无 拉 的 办 法 ， 看 看 结果 如 何 。 

P(r yy =V 1+y 


尤 拉 方程 化 为 


出 此 式 可 见 C 史 1 


C 
“yy et 


即 y “= 万 ( 妃 亦 任意 常数 。 由 此 式 已 经 看 出 y(X) 是 直线 )。 

yr)= Dr+E 

Xvy0)=E=0 

yt1)=D+ FE=1 
QIy(r)] 的 满足 y(0}=0,y(1)= 1 的 极 带 是 直线 y 一 z。 
现在 用 尤 拉 方程 来 解决 变 分 学 史上 最 著名 的 问题 , 即 约 坦 ' 贝 
努 里 在 1696 年 提出 的 “最 速 降 线 问题 ”。 

最 速 降 线 问 题 的 担 法 是 这 样 ， 设 平面 扩 与 地 平面 垂直 ， .4 各 


。25 。 
召 是 此 平面 上 和 任 取 的 两 点 ， -4 点 的 位 置 高 于 已 点 。 全 站 全 


作用 下 沿 曲 线 有 4 访 由 4 点 降落 到 点 ， 所 Y 
需 时 间 为 四。 疝 4 孔 是 什么 曲线 时 ， 时 间 
了 的 值 最 小 【《 设 质点 在 4 点 处 的 初速 为 
零 ) ? 这 里 需 加 说 明 ， 4 点 当然 不 取 在 喇 
点 的 正 上 方 ， 因 为 那样 问题 就 太 简 单 了 。 
解法 ， 取 坐标 系 如 图 3， 记 质 点 的 质 
量 为 说 ， 速 度 为 w， 又 时 间 为 上 ， 则 


于 2 
—fHiyU = 
了 gy 


>»U = 208y 
曲线 5 =2%(z) 的 弧 长 微分 是 
dS=VITy Ed 


又 一 一 -一 几 


dt Vity dr 
v ~ V2gy 


eh 
-Cv()) 
由 此 可 网， i 使 泛 函 
Grey(zD= 人 Ya 
在 此 曲线 y(z) 上 下 极 小 值 即 可 。 现 在 写 出 尤 拉 方 程 
“Flryy -YH 


Ml 


F,=(1+y:)2 (-3) yi 


VV y+y) 


ea 26 ， 


1 3 7 
yy i ea 
此 方程 很 难 友 接 求 解 ， 北 另 怨 办 法 。 
因为 下 = 了 (y,y’)， 所 以 
Fi/=PF,(y,y') 


EF Py + Fy" 
dr 


于 是 尤 拉 方 程 可 表 为 
F,—Tury’ ~— PF,ry’= 0 
现在 可 以 证 明 
-y=C (C 为 任意 常数 》 


EF-y FF, Py’t+ Fryr—y* FP, 


—y (Prvy + Fr. yy’) 
=y'[F,-F'y’ -Fy y=0 


将 已 = 全 扩 一 代入 此 式 得 


1 
-一 一 一 一 一 人 
Vyt1ty’) 

“ay(1+y’ )=D (万 亦 任意 常数 ) 


令 y’=te9， 则 


D 


dy=— Dsin20d8 


dy _—Dsin206de ~2Dsingcos do 
Ne lgg lga 


— N11 + co 


- - 一 -一 -一 一 -- cs i 


[meade 


。27 。 


人 三 sin29+ 到 ( 巨 亦 任意 常数 》 


A 


be 
| 8 = 1 + Cos20) 


这 就 是 尤 拉 方程 的 解 。 
令 20 二 一 9， 则 以 上 二 式 化 为 


i D Tt 
T =3(P -sny)— FH+ E 


t 


D 
\ y= 1 ~ cosp) 


及 当 Pp= 二 0 了 时， 取 工 一 0 一 2， 于 是 


T=r(p— sng) 


y=r(l1—cosg) (和 =- 号) 


所 以 证 画 Q = 全 *W + 风 ?”dx 的 航 带 是 旋 轮 线 ( 也 称 摆 线 )。 可 


以 证 明证 函 Q 在 过 4、 呈 两 点 的 旋 轮 线 上 确实 取 最 小 值 ， 这 一 点 
我 们 不 讲 了 。 

经 过 这 一 番 探 索 ， 我 们 才 知 道 过 A、B 两 点 的 “最 速 降 线 ” 
竞 是 一 条 旋 轮 线 ! 

在 本 节 的 最 后 ， 我 们 来 介绍 更 便于 运算 的 变 分 记号 。 

我 们 将 泛 芒 QLy(z)J 的 变 分 记 作 56&Q， 通 数 y(z) 的 变 分 记 作 
Sy〔 见 本 章 §2 的 第 (4》 式 和 第 《7 》 式 ) 。 于 是 

Tl(y,6y)=6Q 

又 6Q=9'(0) 《 见 第 (8) 式 ) 
在 此 9p(c)= 避 [二 a6o] 


gp'(0) = 人 [5Fs6y+ Fdy)' dz 《 见 p.20) 
而 (6y)’'=6y’@ 


D y= TT) YT) (BUI 87(T) YT = Ur 


。 0D8 。 
b 
§0= | Plr,y,y’)dr 
=| CF,6vt+ FOy’ Jdr a1) 


请 仔细 观察 (11》 式 。 它 说 明 对 于 形 如 | (zx,y,y )dzr 的 泛 本 


已， 我 们 计算 变 分 6SQ， 就 像 计 算 函 数 的 微分 一 样 中 。 
于 是 ， 让 6 已 = 0， 得 


| Fy+ Fdy’)dz= 0 
丰 二 
又 | Fusy'dr=/ Fud(sy) 
nl 
a 
de 


但 6y 是 可 在 相当 广泛 的 函数 类 中 任 取 的 函数 〈 详 情 见 本 节 上 
文 ) ， 所 以 


这 样 来 推导 尤 拉 方 程 ， 书 写 方式 简洁 多 了 。 今 后 我 们 将 用 这 种 办 
法 导出 各 种 泛 函 的 尤 拉 方程 。 
还 要 说 明 一 点 ， 在 这 种 推导 过 程 中 ， 我 们 用 到 了 


dl(8y)=6(dy) (12) 
即 微分 符号 d 与 变 分 符号 6 可 以 互 换 。 事 实 上 
d(6y)=(6u) dr=(6y’ Jdr=6(y'dr)=6tdy) 
所 以 (12》 式 成 立 。 这 种 互 换 法 则 给 变 分 计算 带 来 很 大 的 便利 。 


中 对 于 函数 下 {fz,y (0),Y' (a)),Q 是 参 变量 。 当 5 国定 而 a 变化 , 则 48 下 = FF ,dy 
*+ Fidy’ 


@ 2。 


S 4 关 半 人 cr， Y， YY ee Y) dx 型 泛 
函 间 尤 拉 方程 


在 上 文中 ,对 于 形 如 QCy(z)]= F(z,y,y')dz 的 泛 两 失 


导 了 尤 拉 方程 。 我 们 详 述 了 有 关 条 件 ， 作 了 相当 严 阁 的 论 证 。 
现在 ， 我 们 将 采用 上 节 最 后 所 说 的 方法 来 导 由 汉 两 Q Cy(z)] 一 
[Fssy,y, ,ye )dz 的 克拉 方程 。 对 于 各 种 条 件 只 作 大 略 
叙述 ， 读 者 可 以 自己 补足 这 些 条 件 。 只 要 懂得 了 前 面 几 节 的 内 
容 ， 要 详细 论证 后 文中 的 结果 是 并 不 困难 的 。 因 为 主要 的 难点 已 
经 突破 ， 下 面 这 几 节 都 是 推广 要 了 。 
首先 还 是 针对 较 近 的 日 标 吧 ! 我 们 来 的 导 形 如 QC y( z )3 
= 人 PCz,o 风 gdz 的 泛 画 的 尤 拉 方 程 。 
设 丽 数 忆 (x,y,y',y”) 充 分 光滑 ， 当 泛 函 在 曲线 Cz) 上 
到 极 值 @ 则 变 分 
65Q=0 
但 是 
8Q={ (FytFdy tdv' dr (与 (11) 式 


对 比 ) 
在 上 和 节 中 已 指出 

bh 也 d 

人 rwayraz= -[ 6y dz 
又 


团 指 所 (zw,y,y' ,4 有 对 zf ,yy? tr 的 充分 高 阶 的 连续 偏 导 数 。 其 充分 光滑 
程度 足以 保证 以 下 各 种 运算 顺利 进行 。 
国 ” 指 QCy(z)) 在 满足 边 信和 条件 yCa) = yyy Cyis yb)= ys y’(5) = 


y ' 且 充分 光滑 的 函数 类 中 取 弱 极 值 。 


。30 。 
{Fay'dzr = 人 Fr (063dz 
b b 

=6y’ Fur | -fay -Fdz 
d 

= 一 sy: dit)dr 

《注意 5y'(0)= 0 =56y’(p)) 

:dd dj» 

=- 上 让 (34) ( 才 Fr)az 


= -6 (Fw)| + 


qa? 
J sy Fra F mds 


= sy -Fri(F ud? 
60= 全 Ce， 到 Ey et 0 
和 $3 中 所 说 的 情形 类 似 ， 因 为 y 可 在 相当 广泛 的 函数 类 中 任意 
选取 ,而 ,一 -外 ,，+ -zr 又 是 + 的 连续 丽 数 《 当 y(z) 
取 定 一 一 我 们 设 QCv(z5)] 在 y(rz) 上 取 极 值 )， 世 以 


2 


d d 
i 十 gmt v= 0 外 (13) 


这 就 是 学 六 Q[y(2X)] | Pn ,y*)dz 的 尤 拉 方程 。 当 


泛 函 Q 在 曲线 y(x) 上 取 极 信 ， 且 f(z,y,y ,y ) 及 曲线 y(x) 充 
分 光滑 ， 则 y(z) 满 足 方程 (13》。 

一 般 说 来 ， 用 这 种 论证 方法 可 以 证 明 ， 对 于 泛 前 Q[y(z)j= 
[CR )dz， 设 P(z yy yy 1 ) 对 
VY ,VY "YY 10 春 言 有 充分 高 阶 的 连续 人 篇 导数 (或 说 斑 充 分 兆 
”全 可 参看 本 章 ?中 的 “推广 了 的 预备 定理 ”。 


s 31 。 


少 )。 若 泛 函 在 曲线 y(z) 上 取 极 人 入, 则 y(z) 满 足 如 下 尤 拉 方 程 


和 
FE F ,+ 人 Fs 


t(D, 一 0 《137 “ 
紫 方 程 又 称 克 拉 一 泊 松 方程 。 由 
8Q=) CP ,dy+ Fdy’ +F ,rdy” 
十 .十 PF, Gy" Jadzr= 0 
再 作 若 干 次 分 部 积分 ， 并 利用 
Gy(a)= 0=6y'(a)= = (a) 
6y(6)=0 =6y’ (8)=.… =6y" (6) 
邑 霜 证 贿 (13》 式 。 请 读者 自己 作出 证 明 。 


§5 ”关于 | FELz,y (x), yt), os 
YE) YI TN Yr) Yr) dz 


型 泛 函 的 尤 拉 方程 


现在 ， 我 们 的 叙述 可 以 更 简练 一 些 了 ， 这 样 读者 反而 容易 理 
解 证 明 的 指导 思想 。 


对 于 形 如 QQ[y1(x)， 2 一 人 F(z, yi Ye Y! se, 


2%5)dzr 的 泛 函 ， 设 PF Ory yy 对 YY 
y2 "94 而 言 充 分 光滑 。 若 泛 函 外 在 函数 组 y,(7),…,y,(x%) 上 了 
极 久 “ 指 Q[y,,… ,yy 在 满足 给 定 边 值 且 充 分 光滑 的 函数 组 的 集 
合 中 到 弱 极 值 ) ， 则 y.(x) 满 足 尤 拉 方程 组 


中 指 Qry(iz)) 在 满足 边 值 条件 


goto=en (k=0,1,mn—l 
=e Cy' "ny) 
有 充分 光 少 的 函数 类 中 到 弦 极 信 。 


ss D2 所 
WE 
vi dartv"'=0 (t=1,2,.",n) 《14) 
证 ，“.“6Q=0 又 
et. a 
osQ=| (Fy +t,, 0 


+ Fdy!+ Fdy tt Posdy’ dr 
站 
=-| sv. (入 Fi)dz 
50=[[(F4 -六 Ft) sy +(F,, 
-ey tet CF, 
-Fi)sy. Jdr= 0 


Sy,{ + 一 1 } 2， ) 可 在 相当 广泛 的 函数 类 中 任意 选取 。 
可 取 6y; 二 0 一 6 一 … 一 64%。 这 时 


d 
f(z- 大 Fs.=0 
由 6y, 的 任意 性 及 玉 ,， 一 -人 -下 ,的 连续 性 ， 得 
F 二 


v1 dzx f= 0 
其 余 各 方程 的 证 明 同 比 。 


$6 多 元 函数 的 泛 函 及 其 克拉 方程 


1. 多 元 函数 的 泛 画 
前 面 所 说 的 Q@[y(x)] 是 依赖 于 一 元 函数 y(z) 的 泛 画 。 现 在 
来 研究 依赖 于 多 元 通 数 的 江 函 。 


| [% 巡 


es 933 = 


我 们 来 看 重 积分 
zcesoazay=eczkro] 
0 


一 了 委 Z 魏 1 
在 此 ， 积 分 区 城中 是 由 不 等 式 | 二， 二 1 所 界定 的 正方 形 ， 
Z(z,y) 是 定义 于 上 的 二 元 连续 画 数 。 把 在 咏 上 连 绪 的 二 元 琴 
数 的 集合 记 作 C(D)。 任 取 z(z,y)ECCD)，Q 的 什 叭 一 确定 ， 
所 以 Q[z(x,y)] 是 依 燥 于 z(x,y) 的 泛 函 。 
再 看 一 例 。 对 于 QLz(z，y)]= Sup Z( z，y)D， 将 在 域 
0 去 z 扫 1 人 
7?(_ 1 二， 二 1 外 的 有 界 函 数 的 集合 记 作 M(G)， 对 于 任 取 的 
ZKz,g)EM(G)，Q5Z(z，2) 值 唯一 确定 ， 所 以 局 是 依赖 于 
2 和 zy) 的 泛 函 。 
以 上 是 依赖 于 多 元 函数 的 省 函 的 例子 ， 至 于 一 般 定义 ， 就 不 
必 说 了 。 
9 9 
2 Weve), 全 ]drdy 型 泛 函 的 尤 拉 方 


程 
现在 我 们 来 研究 泛 通 


Qc)I= /fF (sv, s C58) 13, 5]dzdy 
卫 


在 此 ， 设 刀 是 平面 上 的 有 界 闭 区 域 , 且 刀 的 边界 曲线 9D@ 是 光滑 
的 。 考 虑 在 区 域 疡 上 有 一 阶 连 续 偏 导数 的 函数 。 把 这 种 函数 的 集 
合 记 作 C'《D)。 又 在 区 域 力 的 边界 9D 上 给 定 一 函数 p(z,y)。 
把 C1(D) 中 满足 边 值 条 件 z(x,y) ho=9(z,y) 的 函数 z(x,y) 的 
集合 记 作 R1(D)。 即 

RD)={2(7,y); z(x,y) EC(D), zl(zr,y) lap 


名 ”Sup 表示 数 的 集合 的 上 确 界 ， 也 就 是 数 集 的 最 小 的 上 界 。 
名 在 现代 数学 文献 中 ， 常 把 区 域 D 的 边界 记 作 3D， 这 种 记号 来 源 于 拓扑 学 。 


a 3 生 曙 
=9p(r,y)} 
设 泛 范 QC2 (zx, y)] =)Jjr[s, y, 2{7, ,2 5 dzdy 
1 


在 R1(DD) 中 的 函数 z(x,y) 上 到 极 值 。 任 取 2C7,y)€R!'(D), 


令 


62= 2(7,Y) ~ 2(r,Yy) 
这 是 函数 z(z,y) 的 变 人 分。 显然 
G2 lp=[ 2 CX, y)—2(r, Yy)d ln 
=p{(r,y) — Pr,Y) 


和 
则 
-让 二 让- 估 -人 ( 失 )-a 
(注意 -93 二 -92 是 函数 -92 的 变 分 ) 
同 理 


i 


当 z(z,y) 下 6Gz 给 定时 ， 泛 函 @[z+ cz] 显 然 是 zx 的 函数 ， 
记 作 
pla)=Q[z+a'62] 
因为 泛 函 Q@ 在 z(z,y) 上 取 极 信 ， 所 以 阔 数 9(a) 在 a 二 0 处 
到 极 值 (参看 本 章 §$3) ， 于 是 
p'{0)= 0 


但 是 QEzCz,5)I=)f F(z,y,2,p, gdrdy 
DD 


以 z+a6z 取 代 2 有 时，p 由 p+2*6p 取 代 ，9g 由 9+ar69 取代。 


所 以 
va)= F(z,y,2+0"62,p+tadp,9 
让 
+a:é9q)dzrxdy 
了 w0= 作 FF.6z+ FSp+ F,Sq)drdy (15) 
2 


”和 本 章 旨 2 第 4 节 中 所 说 的 一 样 ， 可 以 证 明 敬 泛 了 前 QLzCz,y)] 
= 上 FG,y,z,p,q)dzdy 的 变 分 6Q 存 在 ， 则 
Pp’(0)=60 


, sQ=)(F .8+ FSp+ F.9drdy (16) 
总 


这 是 一 个 很 方便 的 公式 〈 请 读 老 将 此 式 和 《11) 式 作 一 比较 》。 
由 2 (0) 王 0， 得 


sQ=)|(F.82+ Poap+F69)drag= 0 7) 
但 是 
| FP, 62) = dz-(F)+F, (2) 
二 3 
一 82 3 Fdp 
( 符号 5 与 -人 -可 交换 ) 
Fdn=-37 oo (Fd2 ) -8 
同 理 -F802) = 62 + Fo8g 


| os 
ee a ， 之 Zz)— - 8 


于 是 


| ("Spt+ Fgdrdy 
| (P62)+ oy 


I CF) + (EF) Jezdrdy 


由 (F182) By + FG) Yady 
= | (Fo6ady— Fedzdz) 
= fir dy Pedz)s2- 0 
a (Fdp+ FS9drdy 
=-J[ 友 二 汪汪 -Pu) ]`5zdzay 
Ov 


| F.- 训 F,— 20 F, J ye 0 
此 


DO 人 
因为 -3 了， 一 -6y 7 连续， 且 52 可 在 相当 广泛 的 
函数 类 中 任 取 ， 所 以 


On -9 11 


D 在 这 里 用 了 格林 公式 :人 (3 + v Jdzdy = om- udr) 


* 37 » 


人 
Ox 


> 
-高 


这 就 是 泛 函 Q [F(z,y,2,p,9)drdy 的 尤 拉 方程 一般 说 来 


F,— FPF, F,=00 (C18) 


”这 是 偏 微分 方程 。 此 方程 的 解 ， 也 称 为 “ 极 带 ”。 如 果 泛 函 Q 
=fF Gz,y,z,p,o)dzdy 在 (x,y) 上 取 极 值 ( 指 弱 极 信 ) ， 则 


z(x,y) 应 当 满足 方程 (18》。 


例 考虑 泛 函 
qQcz(z,23=|[( S2 ) + (2 yy ]ezdy 
={|(p:+ gdrady 


2(T,8) ip 一 OCT VD) 
P(x,y) 是 给 定 存 边界 上 的 连续 函数 。 如 果 此 泛 函 QQ 在 次数 
z(x,yV) 上 到 极 值 ， 那 么 z(x,y) 应 当 满足 什么 形式 的 方程 式 呢 ? 
因为 P(x,y,2,7,9)=p*+ a? 
“P=2p, Fo=2g 
代入 尤 拉 方 程 (18》， 得 


$+ 全 了 一 0 (在 万 内 处 处 成 立 ) 
EE . (19) 


方程 -2 + SS -0 是 著名 的 拉 下 拉 斯 方程 


作对 于 瑚 阴 QCutzieyyzo)y] = fe fF (se rst sr 


3 和)aximdz, (DD 是 n 维 空间 中 的 财 区 域 ， 下 充分 光滑 。) 可 证 朋 ， 若 Q 在 满足 一 定 
i “zr) 上 取 级 什 。 则 4 满足 克 拉 方 程 下 ,~ 六 2- 


F 二 Dr， 
人 (187) 


9 38 « 
求 函 数 >(z,g%)， 使 它 在 闭 区 域 万 (万 =DUaD) 上 连续 ,在 
开 区 域 九 内 满 足 拉 下 拉 斯 方程 ， 在 边界 3 万 上 取 给 定 的 函数 值 , 
这 种 问题 称 为 “ 拉 下 拉 斯 方程 的 第 一 边 值 问题 。 此 例 说 明 ， 求 泛 


画 Q= 人 [人 3) + (3 六 7ray 航 值 的 问题 可 转化 为 解 拉 、 


L 


拉 斯 方程 的 第 一 边 值 问题 。 反 过 来 ， 在 一 定 条 件 下 ， 解 拉 卜 拉 斯 
方程 的 第 一 边 值 问题 也 可 转化 为 解 上 上述 泛 函 QQ 的 极 值 问题 。 十 九 
世纪 的 大 数学 家 黎 曼 (Riemann) 曾 利 用 这 种 观点 来 证 明 拉 直 拉 
斯 方程 的 第 一 边 值 问题 的 解 存在 。 他 的 方法 很 有 启发 性 ， 可 是 在 
他 当年 的 论证 中 有 一 个 漏洞 中 。 后 来 ， 另 一 位 数学 家 外 尔 斯 特 拉 
斯 《Weierstras， 此 人 可 说 是 十 九 世 纪 中 最 严格 的 数学 家 ) 指出 
了 这 个 毛病 。 


再 举 一 例 。 对 于 泛 函 Q [zz,y)] = 人:(z,o dzay = 


dz zr,)drdy 又 z(z;, 纪 | =0(7,)=3 问 此 泛 


通 在 什么 函数 2(z,y) 上 取 极 值 〈 在 连续 函数 类 C( 刀 ) 中 考虑 ) ? 
显然 Q[z(z,y)j 之 0， 且 当 z{z) 三 0( 在 但 内 )Q[z(z,y)] 
取 最 小 值 0 。 
但 在 如 =DU3D 上 连续 且 在 品 内 为 零 的 通 数 ， 是 不 会 在 边 
办 6 站 上 取 非 零 值 的 。 所 以 z jso= 3 是 不 能 成 立 的 。 这 就 是 说 ， 
此 泛 函 在 C(DD) 集 合 上 不 能 达到 极 值 。 金 


QD 德 思 数学 家 歼 可 1851 年 在 他 的 博士 论文 中 论证 了 关于 共 形 映射 的 基本 定理 。 
在 这 篇 论文 中 ， 他 不 加 证 机 的 认为 有 函数 wz, 纪 使 天 耿 Qi (w+w 
dr dy 腺 极 值 。 , ? 

多 显然 汉 画 人 2? (zyg) dzdy 可 以 无 限 地 接近 零 ( 当 Cz,9) 在 集合 CCD) 中 


恋 动 】，、 但 不 公 达 轴 堆 。 用 分 析 的 诸 言 来 演 ， 当 z(zyy) ECID),Qtz (Cr,Y)3J 有 下 
确 界 过 ， 阳 在 C 2 上， 巴 不 能 达到 下 确 界 入 。 


em 39 。 


S7 本 章 提 要 


一 、 基 本 定义 

1.、 汪 函 ”把 具备 某 种 性 质 的 函数 的 集合 记 作 娓 ， 若 对 于 任 
何 / 5 DD， 变 量 QQ 都 有 唯一 确定 的 值 各 它 对 应 ， 则 称 变 量 QQ 为 依 
束 于 少数 了 的 泛 画 ， 记 为 Q=Q[ 有 或 9 一 QLf(z)J。 一 元 函数 
/Tz) 也 可 换 成 多 元 函数 ， 如 YQ 一 QCz(r,y)] 

2. 变 分 

《1)》 函数 的 变 分 
对 于 一 元 函数 yf4z)， 规 定 它 的 变 分 为 

Sy= Yr) yr) 


在 计算 时 ， 符 号 4 与 5 可 互 换 ， 即 
d(5y)=6(ay) 
对 于 多 元 娩 数 z(z,y)， 规 定 它 的 变 分 为 
Bz2= ZzZ(r,y)—2(r,Yy) 
人 -分 =p， -和 2 = 9。 在 计算 时 有 
7 oz ” Oy 3 


3 (82) =6()-p 


训 (62)=6( 注 )=64 


(2 ) 泛 函 的 变 分 
u) 证 函 Q5y(z)3 (或 9[z(z,y)] 的 增 量 的 线性 主要 部 分 叫 
变 分 . 记 作 SC。 
5) 若 运 函 QCy(z)]= Cz,y(z),y'(z)Jdzr 的 变 分 存在 
《 没 玉 充分 光滑 》， 则 
$Q=9(0)=| (Fy+ Fy, dy Ydr 
《在 此 ，P(4) 二 QLy(7)t+a6y]，y(z) 及 6y 给 定 》 


- 40 =» 


2 若 汉 也 Q(z(x,5)]=fP (ya(ry), 2 


‘G2 
-3 jazdgy 的 变 分 存在 《 设 严 充 分 净 滑 ， 并 记 -22_ = p， =- 


= 4 ), 则 
sQ=|f FB.t Fdpt+ Fgqdrdy 


3. 万 拉 方 程 
(1 》 对 于 泛 函 Q[y(z)] = [Frz,y(z),y’(z)Jdz， 其 


尤 拉 方程 为 


C2) 对 于 泛 函 QCy(z)]= 有 FE, yn), vy’ (zy 


y (zx)jdx， 其 尤 拉 方程 为 
df "od 
FF Pr tl) drt 


一 性 
(3 》 对 于 泛 函 数 Q[gi(z) ,yz)] = 人 FCz, yy， 
81 人 dz， 其 尤 拉 方 程 为 
Fs -fF 0 (i=1,,n) 


C4) 对 于 泛 画 QCaCz,5)] = 站 FP [x,y (zy) -92 
D 
0 
-5 ]4rdy， 其 尤 拉 方程 为 
9 0 BE GE 
Ds Pr (my der) 
《5 》 对 于 泛 函 OCR(Cziyza yzo) 


9 41 。 


i Qu 
-J 人 Cr， Ox 了 


Ou 
A ]az…-az。 


其 尤 拉 方程 为 


oml 
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第 三 章 “ 变 分 法 在 力学 中 的 一 些 应 用 


在 第 二 章 中 叙述 了 泛 岗 介 取 极 什 的 必要 条 件 ， 慨 括 起 米 就 
是 ， 若 泛 函 日 在 某 函 数 上 取 极 值 《此 函数 满足 事先 给 定 的 边 值 条 
件 ) , 则 该 函数 满足 尤 拉 方 程 。 很 明显 ， 我 们 还 要 解决 另 一 个 问 
题 ， 变 分 理论 才 会 完 束 。 这 问题 就 是 ， 如 果 某 函数 满足 尤 拉 方 程 
和 边 信条 件 ， 它 是 不 是 能 使 泛 函 名 报 极 值 ? 此 外 、 还 应 当 说 明 一 
全 问题， 如 果 Q 在 某 函 数 上 取 极 值 ， 那 么 Q 究 竞 取 极 大 信 还 是 极 
小 倩 ? 从 逻辑 学 的 角度 来 观察 ， 应 当先 搞 清楚 了 这 些 问题 ， 再 来 
谈 变 分 法 的 应 用 。 但 是 ， 对 初学 者 来 说 ， 作 者 认为 在 学 了 尤 拉 方 
程 之 后 ， 还 是 先 接触 应 用 问题 为 好 。 理 由 是 ， 

1， 广大 的 工程 技术 人 员 关 心 的 是 用 变 分 法 解决 实际 问题 。 
对 这 些 同志 来 说 ， 只 要 掌握 了 本 书 第 一 、 二 章 的 内 容 ， 就 可 以 用 
于 实际 了 。 

2. 从 变 分 学 的 发 展 史 来 看 也 是 这 样 。 尤 拉 和 拉 格 兰 日 的 工 
作 只 提出 了 泛 函 取 极 值 的 必要 条 件 《〈 即 本 书 第 二 章 的 办 容 ) 。 尤 
拉 是 在 1736 年 发 表 他 的 成 果 的 ， 其 中 有 后 来 命名 为 尤 拉 方程 的 著 
名 方程 式 。 在 此 后 七 十 余年 间 ， 数 学 家 们 长 期 探索 泛 函 日 到 极 值 
的 充分 条 件 ， 没 有 取得 多 少 成 绩 。1807 年 雅 可 必 〈Jacobi) 在 这 
方面 稍 有 进展 ， 但 他 对 问题 的 叙述 是 不 准确 的 。 直 到 1879 年 ， 外 
尔 斯 特 拉 斯 才 正确 提出 并 严格 论证 了 谤 晒 取 极 值 的 充分 条 件 。 这 
已 经 是 尤 拉 方 程 出 现 后 的 第 一 百 四 十 三 个 年 头 了 ! 尽管 充分 条 件 
的 建立 是 这 样 晚 ， 但 是 并 没有 妨 开 变 分 法 在 这 一 个 半 世 纪 中 广泛 
应 用 十 实际 。 

3. 泛 函 极 值 存在 的 充分 条 件 的 论证 是 相当 困难 的 《本 书 限 
于 篇 幅 ， 不 谈 这 个 充分 条 件 了 ) 。 

本 章 主要 是 介绍 变 分 法 在 力学 中 的 应 用 ， 其 他 方面 略 有 涉 
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及 。 读 者 理解 指导 思想 之 后 ， 自 然 可 应 用 到 宽广 泛 的 领域 中 去 。 


831 哈密 和 尔 顿 原理 和 拉 格 关上 日 方程 


1. 了 哈密 尔 顿 原理 

当 和 牛顿 建立 了 以 三 大 定律 及 万 有 引力 定律 为 基础 的 力学 理论 
以 后 ,无 数 的 自然 现象 都 得 到 了 定量 的 说 明 ,事情 似乎 已 经 很 完 
了 。 但 科学 家 是 永 不 满足 于 已 有 成 就 的 人 ， 他 们 仍然 继续 探索 宇 
窗 的 奥秘 。 岳 来 ， 拉 格 兰 日 在 十 八 世 纪 提 出 一 个 变 分 原理 ， 从 这 
个 原理 出 发 ， 运 用 变 分 法 ， 能 够 十 分 方便 地 解决 力学 问题 ， 并 所 
还 可 以 推 蜡 时 力学 中 的 主要 定律 。 在 他 著 的 《 分析 力学 》 中 叙述 
了 这 些 成 果 。 据 说 他 诬 过 这 么 一 句 话 ，“ 在 我 的 &《 分 析 力 党》 中 
没有 一 张 图 ”。 他 还 创立 了 “ 拉 格 兰 日 运动 方程 ”， 它 比 年 顿 的 
运动 方程 适应 的 范围 更 广 ， 用 起 米 也 往往 更 方 使 。 以 后 ， 喻 密 尔 
顿 (Flamilton》 又 发 展 了 拉 格 兰 日 的 理论 ， 他 在 1834 年 提出 了 
一 个 著名 的 原理 现在 称 为 哈密 尔 顿 原理 》， 大 意 是 ， 

在 质点 《甚至 是 质点 系 或 物体 ) 的 一 切 可 能 的 运动 中 ， 真 实 
的 运动 应 使 积分 


s=/ -Dod 


取 极 信 。 在 此 ， 不 和 UU 分 别 是 动能 和 位 能 ，t, 各 1, 是 两 个 任 取 的 
时 疗 。 

这 个 原理 竞 成 了 力学 中 的 基本 原理 。 以 它 为 基础 ， 可 导出 和牛 
顿 三 大 定律 和 能 量 、 动 量 和 动量 虐 守 恒定 律 。 读 者 车 对 此 有 兴 
趣 ， 请 阅 朗 道 所 著 《 力 学》 一 书 (学 了 变 分 法 以 后 ， 看 此 书 是 不 
会 感到 困难 的 ) 。 现 在 我 们 来 把 此 原理 说 得 更 清楚 一 些 。 

哈密 尔 顿 厌 理 的 精确 表述 是 ， 假 定 在 += 及 t = 时 刻 质点 
的 位 置 已 分 别 确定 在 4 点 及 点， 那么 质点 运动 的 真实 轨道 及 速 
度 ， 使 积分 


ss 44 »* 


Ss=) "(7T-U)dt (5 又 称 “ 作 用 量 ”》 


取 极 值 。 在 此 了 及 过 各 为 质点 的 动能 和 位 能 。 所 谓 真 实 轨 道 请 看 图 
4, 图 中 的 > 二 (zx1,7x,,7X3) 是 三 维 空间 坐标 , + 是 时 间 坐 标 , 带 箭头 
的 曲线 表示 真实 轨道 。 这 是 示意 图 。 

例 ” 试 由 了 哈密 尔 顿 原理 导出 和 牛顿 
第 一 定律 

对 于 自由 质点 ， 动 能 

mr ‘ -。\ 卫 
和 = 了 (2 二 好 二 


位 能 U 恒 为 零 。 于 是 作用 量 


X= {Xi, Hz, A 


Ss=) Bt yt)d 
当 t 二 tf， 质点 位 和 于 Atai,61,01) 点 ， 
t= 质点 位 于 (azy,pa,cz) 点 。 即 
z(t1)=a) r(t,)=as 
y(t1)=b, y(t,)=6, 
2(t1)=e 2(12)= 6, 


显然 ,此 S=| 扣 (za 二 大 二 ze)dt 一 SCz0D ,y(t), 24)] 
是 轨道 


T=x(t) 
: Y=y(t) 
| z =2(1) 


的 泛 项 。 这 是 第 二 章 和 5 中 所 说 的 | 
V1(z), 下，Yi(ZT) Jdx 更 泛 函 。 由 $5 可知， 若 泛 通 S 在 轨 


@@ (Dy(tt) z(t) 表示 出 质点 举 称 (2,Y,2) 与 时 间 ! 的 关系 。Z ，# ，z 分 别 


过 a d. 
是 -9 革 ，-3，-9 荆 ， 这 是 牛顿 引 人 的 导数 《他称 之 为 “ 流 数 ” 》 记 号 。 


s A5 。 


| r=r(t) 
道 ; 4 二 y(t1) 上 有 取 极 值 ， 则 z(t) ,y(t) 和 z(t) 满 足 尤 拉 方 程 
E(t) 
| d 
F Te 一 人 
d 


\ 
但 是 == 守 (2+ V7) 


ey ie dq 
“Fs: 二 0 Fi=mr7 dt F. 


1 


“, 工 二 0 同 理 得 y=0,，2z 一 0 
可 见 自 由 质点 的 加 速度 为 零 ， 所 以 其 速度 为 常量 〈 说 得 更 清楚 一 
些 ， 速 度 是 常 矢 量 ) 。 也 就 是 说 ， 当 质点 不 受 外 力作 用 时 ， 它 或 
静止 ， 或 作 匀 速 直 线 运动 。 我 们 由 哈密 尔 顿 原理 导出 了 牛顿 第 一 
定律 。 

例 试 由 哈密 尔 顿 原理 导出 质点 在 重力 场 @ 中 的 运动 方程。 

在 地 球 表面 附近 ， 可 认为 重力 加 速度 9 的 值 不 变 。 将 坐标 系 
原点 O 放 在 地 平面 上 , 取 z 轴 坚 直 向 上 。 质 点 的 动能 7 了 == 全 (x*+y 
+ z*)， 位 能 U =mgz， 二 是 泛 函 《 即 “ 作 用 量 ”)》 

S =|, (Fa+ y+ 2) —mgz Jat 


Fo (+ Yd —mgz) 


F=0=F, FPF,=—myg 


FF:=mi, Fy=my, FF:=mg 


出 ” 指 地 球 表 商 附 近 的 均 尺 前 力 场 。 


md 用 。 


由 P.-EF:=0 得 2 = 0 


由 FP, fF 得 #= 0 
由 FP.~ pF: 0 得 mg 一 ms 二 0 一 一 9 
于 是 得 质点 在 重力 场 中 的 运动 方程 
| z=0 
#0 
| s—=-g 


例 由 险 密 尔 顿 原理 导出 质点 在 地 球 引力 场 中 的 运动 方程 。 
将 坐标 原点 置 于 地 球 中 心中 ， 质 点 及 地 球 的 质量 分 别 记 为 严 
和 M， 质 点 的 全 标 为 《z(1)，y(t)，z(t)) ， 于 是 动能 及 位 能 
分 别 是 
了 = 本 (入 十 办 十 知 ) ， 


U=-GM"--KT rvVriryta 


作用 量 为 
S -| en 十 东 ? 十 过 3) 十 KY Jat 


= FEE sO YC) ,201) 20) 901), 201)]dt 


nm -2 [4 > ni 
ss 十 光 2 十 记 2) 十 A 


| mK mK 


FP.= OrVriyyite ~ rs 7 es 
F,= -2 Py=my 


合 严格 地 说 ， 此 分 标 系 的 宇 条 纳 强 还 未 能 印 地 球 牢 向 地 连接 在 一 起 ， 关 于 这 一 
参阅 爱 因 括 坦 所 著 改 狭义 与 广义 相对 论 引 说》 中 译本 第 10 抽 “人 其 到 上 略 学 标 系 ” 一 


点 请 
节 。 我 们 不 妨 其 定 爸 术 系 的 z， 纪 ，2 三 提 正 向 均 指向 三 个 充分 遥远 的 天 体 。 . 


»- A7T * 
mK 


3 


五 .= 一 2 FF;=ms$ 


人 


尤 直方 程 〈14) 化 为 质点 在 地 球 引力 场 中 的 运动 方程 。 


wo 

二 一 开 
二 一 

乡 一 一天- 

2 

| #=—K— 


2. 质点 系 的 广义 坐标 和 广义 速度 

确定 质点 在 平面 上 的 位 置 ， 需 要 两 个 坐标 ， 这 是 常识 。 但 是 
请 注意 ， 这 两 个 坐标 不 一 定 是 直角 坐标 。 例 如 ， 在 平面 上 卫 极 堂 
标 系 ， 则 质点 位 置 由 坐标 《向 径 ) 积 P? (角度 ) 所 确定 。 由 


| T=rCO8P 
| Y=rsing 
我 们 把 r 、92 称 为 质点 的 广义 坐标 。 
一 般 说 来 ， 若 有 两 个 独立 参数 q,，9:， 县 
T= TI(d, 9a) 
y=V(qw 9s) 《 即 上 由 (91，9:) 可 唯一 -确定 (XT，、y)) 
网 称 9/，9: 是 平面 上 点 的 广义 坐标 。 
对 于 空间 中 的 质点 ， 它 的 直角 坐标 是 《z ,2y,z) ， 若 有 狸 
立 和 参数 gq,，9;。9s， 且 


| X= XI(gq qr 9s) 
] 2 一 2(C9u 92， 9s) 
| 之 一 全 (qi sz， 93) 


( 即 由 (C9,，9:，94s) 可 唯一 确定 (x ,5# ,2)) 刚 称 91，9;:，9s 是 
空间 中 质点 的 广义 坐标 。 

决定 空间 中 一 个 质点 的 位 置 要 三 个 广义 坐标 。 要 确定 入 个 质 
点 的 位 置 常 需要 3 闪 个 广义 坐标 。 一 般 说 来 ， 能 够 确定 体系 〈 质 


* A 。， 


点 系 或 刚体 ) 位 置 的 KK 个 独立 参数 9,，9:，*…9x 称 为 此 和 体系 的 
“广义 又 标 ”， 而 它们 对 时 间 的 导数 ，9:，…Sx 称 为 广义 速 
度 。 数 到 称 为 体系 的 自由 度 。 

3. 拉 格 兰 日 方程 

现在 ， 我 们 来 导出 描写 质点 运动 的 拉 格 兰 日 方程 。 先 设 质 点 
只 有 一 个 广义 学 标 9 《这 时 质点 的 自由 度 为 1) 。 因 为 质点 位 置 
由 广义 坐标 9( 纪 决定 ， 所 以 动能 全 和 位 能 避 是 9 及 4 的 区 数 。 把 
TT 一 局 叫 微 拉 格 兰 日 函数 ， 记 为 

Tsar 

于 是 ， 质 点 的 作用 量 为 


S= 人 0 gq, gd 
因为 取 极 值 @， 所 以 真实 轨道 9(1) 注 足 
Sd -7=0 《20) 


这 就 是 力学 中 著名 的 拉 格 兰 日 方程 。 辐 样 ， 车 质点 系 〈 在 此 处 ， 
质点 系 包 括 刚体 ) 的 位 置 由 广义 郁 标 91,，9:，…，9x 决定 《这 
时 称 质点 系 有 天 个 自由 兰 ) 且 g (i) 及 q(ts) (i =1,,K) 
均 已 给 定 〈 即 在 + = 所 及 + =1, 两 时 刻 ， 体 系 的 位 置 已 给 定 ), 当 
质点 系 由 机 时 刻 的 位 置 变 到 ts, 时刻 的 位 置 时 ， 作 用 量 


Ss=/ ?Lt,g, gr ,G1 dx)dt 


取 极 值 。S 是 第 二 章 $ 5 中 所 说 的 | FEz，y,(z),，…，y,(7z)， 


jz)，… 思 (rz)]dz 型 泛 函 。 由 此 可 知 真 实 轨 道 9,(1) (=:1， 
…， 上 此》 满足 


0 (i=1,, KK) (21 》 
这 就 是 质点 系 的 拉 格 兰 日 方程 组 。 它 是 在 广义 坐标 系 中 质点 


凶 指 当 ?9(fD) = 9)，9(t4) 一 Qi: 给 定时 ， 质 点 的 嘉实 轨道 9: (位 及 9,(1) (这 是 
洽 足 gtfi) = 44 gf) = 2s 的 一 切 可 能 轨道 中 的 一 个 ) 使 S 取 极 值 。 
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系 的 运动 方程 ， 表 达 了 质点 系 运动 的 一 般 规律 。 
读者 也 许 会 间 ， 拉 格 兰 日 方程 究竟 有 。 

什么 优点 昵 ?我们 还 是 举例 来 说 明 问题 

] 

例 在 长 为 1 的 细 丝 《其 质量 略 去 不 
计 ) 末端 悬 一 质点 ， 其 质量 为 m， 细 丝 另 
一 端 冉 定 ， 置 重力 场 中 。 写 出 质点 的 运动 
方程。 

到 坐 标 系 如 图 5 所 示 。 此 质点 的 位 置 显 然 由 一 个 广义 坐标 ? 
确定 。 动 能 了 = 于 mo 了 miCD)2， 当 p- 0 时 ， 取 质点 的 位 能 
为 零 ， 则 在 2 处 质点 的 位 能 为 

U=mg: AC =imgil( 1 ~ cosp) 
于 是 拉 格 兰 日 函数 为 


L=7-— U 一 -mg 一 IO 1] — cos®) 


些 
吧 


质点 的 拉 格 兰 日 方程 (20)〉 即 


OF df a 
由 
UO 2 A 
0 lnigsing 有 J (ml) = lg 
可 得 
9 + 了 sing = 0 (22) 


这 就 是 单 捍 的 运动 方程 ( 设 摆 杆 的 质量 可 上 略 去 ) ， 这 是 非 线 性 二 
阶 常 微分 方程 。 若 | 2 | 充分 小 ， 用 sinPp 守 89 取 

9 + =0 (23) 
代替 .上 述 方程 来 描述 单 摆 运 动 ， 这 是 线性 二 阶 常 微分 方程 。 它 的 
通 解 为 


。50 。 


P= cicosy T+ + cSin 1 t 


主 


Pp(0)=0 9(0 )=.4 


oAV sn ft 
可 见 当 质点 作 微 小 摆动 时 ， 其 周期 是 2 六 二， 但 当 扫 幅 太 大 ， 


其 周期 就 不 是 此 数 了 。 

网 闭 我 们 对 摆动 作 了 一 番 分 析 。 贝 此 例 已 经 能 看 出 拉 格 兰 日 
方程 的 户 大 优越 性 。 它 适用 于 各 种 形式 的 坐标 系 。 在 解 应 用 题 时 
不 必 画 力 的 分 解 图 ”只 须 写 出 动能 了 及 位 能 尽 关于 广义 坐标 9, 及 
广义 速度 q,， (i 二 1，…，&)》 的 表达 式 ， 然 后 代入 拉 格 兰 日 方 
程 即 可 。 

在 这 里 我 位 还 要 说 咀 ， 因 为 本 书 只 是 简略 地 介绍 一 下 变 分 法 
在 为 学 中 的 应用 ， 所 以 对 于 拉 格 兰 日 方程 成 立 的 条 件 〈 即 对 什么 
样 的 质点 系 ， 它 的 运动 方程 可 用 拉 格 兰 日 方程 描写 ) 没有 作 和 精确 
论述 。 读 考 欲 知 其 详 ， 请 参阅 项 道 所 著 《 力 学 》 或 国内 近年 来 出 
版 的 有 关 专 莹 。 


由 


8$2 变 分 法 在 力学 中 的 应 用 举例 


在 处 理 和 名 种 工程 技术 问题 时 ， 常 常 需要 推导 出 描写 菜 些 物 体 
运动 的 微分 方程 ， 然 后 科 用 微分 方程 的 解 来 找 出 运动 规律 。 这 个 
办 法 是 牛顿 创立 的 。 他 找到 了 描写 天 体 运动 后 来 看 出 还 适用 于 
各 种 物体 ) 的 微分 方程 ， 然 后 断定 天 体 运动 轨道 是 加 锥 曲线 。 另 
一 个 著名 的 例子 是 麦克 十 威 〈Maxwell) 找到 了 描写 电磁 场 变化 
的 微分 方程 ， 然 后 根据 此 方程 断定 电磁 波 存 在 ， 这 是 1868 年 的 
事 。 直 到 1887 年 《也 就 是 在 麦克 士 威 根据 微分 方程 作出 他 的 科学 


» 51 。 


预言 的 十 九 年 之 后 ) 蔡 兹 (Hertz) 才 用 实验 证 实 了 电磁 波 的 存 
在 。 在 科学 史上 还 有 大 最 的 例子 可 以 说 明 微分 方程 的 巨 
大 作用 。 在 我 国 工 程 界 ， 现 在 大 家 习惯 于 用 达 朗 倍 尔 床 
理 列 晶 力 平 衡 方程 ， 出 此 得 到 运动 方程 。 这 种 方法 固然 
有 它 的 优点 ， 但 是 若 用 变 分 原理 来 导出 方程 ， 却 有 更 大 
的 优越 性 。 在 本 节 中 ， 我 们 将 从 各 方面 举例 来 说 明 这 种 
方法 。 

例 推导 细 杆 的 纵 振 动 方程 。 

取 华 标 系 如 图 6 所 示 。 将 杆 上 点 xz 处 的 位 移 记 为 个 
u(xX，t)， 由 虎 克 定律 ， 长 |! 的 杆 在 力 了 的 作用 下 ( 力 下 与 杆 平 


行 ) 所 发 生 的 长 度 改 变 A1= 疡 <-， 在 此 玉 为 杨 氏 模 量 ，S 为 村 
的 横断 面 面 积 。 于 是 相对 伸 攻 


! a RY . Ee ry 
er 2 eae 1 一 LSe 
现在 截取 村 段 如 图 7 所 示 。 


则 
Tt IrT+u(rt+ dr,1)— ru(r, 1) 
六 


u(rT+t IT,t)— Hx, 1) 
rx 


—u,(r+ O47,1) (0<e<1) 
Flr+tAr, 1)=ESu,(rtOsIr, 1) 
在 此 式 中 ， 令 x 一 0， 得 
Fz, 1)= ESu,(z, 1) 
Flr+Ar, 1)= ESu (r+ /Sr, +1) 


现在 求 弹性 位 能 吕 。 由 弹性 力学 可 知 ， 弹 性 栖 ZX+4AX 
在 变形 过 程 中 克服 内 力 所 作 的 功 ， 作 为 能 量 积 鞭 在 
体内 ， 此 能 量 称 为 弹性 位 能 。 杆 段 [x，z 二 x) 
所 受 弹 性 力 为 《〈 见 图 7》 图 了 


[3 


Fixt) 


4 


F(x+4X,t) 


。52 。 
SEu(lrt Ar,t)— Sh (r,t)— SEu,.(r+0sr,t) Ar 
. 《0<6g< 1) 
于 是 ， 在 dt 时 间 内 ， 此 力 所 作 的 功 为 
SEvu,, (r+oAr,t} Ar udt= Su,(r,t}udrdt 
这 是 在 杆 段 Cr，x+ rz] 中 内 力 所 作 的 功 。 至 于 弹性 体 克服 内 
力 所 作 的 功 ， 当 然 与 此 式 反 导 。 所 以 在 [t，++ It1J 时 间 亲 隔 
中 积 葛 在 杆 段 [rz，x + 人 zx] 中 的 位 能 微 元 是 
dU=— SEu,(r, tudrat 
于 是 ， 弹 性 位 能 为 
U= -| | SEu..dz, tudrdt 


SF 全 fo _ 
=— SE ja 全 二 (UH)— Hd,, ]zz 


=— Sk | | nds ja 


全 页 [sda 

a | 

=SE | dz (Hr )at 

-SE | 计 s2d2© (24) 
再 求 动能 。 在 杆 段 [Tz，zx+ xz] 中 ， 动 能 微 元 是 

dT= 广 Plz)ui(z, tJ)Sdr (pl(z) 是 质量 密度 ) 
所 以 动能 是 
T=| 诗 PCr)SuiCz, t)dz C25) 
拉 格 兰 日 函数 为 
L=T- U=| (or)Su 一 于 Sa )dz 


9D 在 推导 过 程 巾 ， 用 到 了 字 # C0, 1)=0, wtl, 机 = 0 及 Noz 0)= 0 


s 虽 人 e 
作用 最 为 
S=| 7-cndi 


一 | 站 1 :1 2 . 
=| | ($e)s SEEa3 ]dzd; 
在 此 ，S 是 zz， 好 的 泛 柄 ， 即 
Suts, D1) | Prenat 


2 {r,t dradt 
这 是 在 第 二 章 8 6 中 讨论 过 的 多 元 函数 的 泛 函 ， 它 的 尤 拉 方 程 ( 参 
看 (18) 式 ) 是 


8 3 
i 


有 Sed :ol 2 一 工 2 _ ? 
, F=3P(r)uiS 3 Et 2 P(r)Sg 3° EP 


f=0 《4 一 六 4 一 9)》 


已 ,一 0 FPF,=—-SEp=—SEru, 
Fo=p(r)Sgqg= pl7r)3u, 
8 0 


Ox 


(SEu.) -OCp(r)Su)= 0 


4 (注意 模 截 面 面 积 S 是 常数 ) 《26) 


当 p(z) 三 P (常数 时 ， 取 汪 二 a?), 得 


Hit = QU, (27) 

这 是 偏 微分 方程 中 熟知 的 弦 振 动 方 程 。 

例 推导 弦 的 自由 横 振动 方 程 | 
( 设 驴 两 端点 固定 》 

取 坐 标 系 如 图 8 。 用 函数 u(z， 
1 表示 在 时 刻 上 点 z 的 横向 位 移 。 

取 喜 发 〔[z，z+ Ix ， 此 段 经 
过 变形 后 的 长 廊 是 图 
dy=Y1+2a3 dz， 所 以 此 耻 的 伸 长 量 是 


和 54 由 
Vitu edzr 一 dr= (VTS ~1)dr 
~( 1 + 去 #2 一 1 ) dr -uidr 


《C 当 zz? 之 1 ， 即 兴 比 1 小 得 多 ) 
在 此 眉 中 ， 弹 和 性 位 能 与 伸 长 量 成 下 出， 即 


dU ds C28) 
于 是 位 能 
} 
二 
v=| 3 uidr (29) 
又 动能 为 


T =| 二) ud (p(x) 是 缠 的 质量 线 密 度 ) (30》 


作用 量 为 
s=[ {doc ur, 0) -Kuilz, 站 ]azdt 
这 是 多 元 函数 u(rz，t) 的 汉 函 ， 它 的 蕊 拉 方 程 《 即 方程 (18)) 是 
CKuilr, 人 一 -SCPCz)a (zt) j=0 
K 


EE 
这 就 是 方程 《26) ， 当 P(z) 择 P( 常 数 )， 取 从 一 ax*， 我 们 又 得 到 


方程 (27》 。 可 见方 程 (26》 和 和 (27) 既 可 以 措 写 细 杆 的 纵 振 动 ， 
又 可 描写 弦 的 横 振 动 。 

例 ”推导 弦 的 受 连 俩 振动 方程 。 
问题 的 条 件 与 前 例 基本 相同 。 叭 一 区 别 是 设 装 上 各 点 均 有 外 力作 
用 ， 此 力 与 弦 的 平衡 位 置 ( 即 z 轴 ) 垂直 ， 外 力 的 线 密度 为 
p(x,t)@, 于 是 作用 在 弦 段 [z,z+ -Az] 上 的 总 外 力 是 p(x, 1)dz。 
在 时 刻 +， 获 上 的 z 点 位 移 到 4u(z，+)， 于 是 在 此 蓄 自 上， 外 力 
作 的 功 可 近似 地 表 为 

一 及 (ZE 要 并 

”“ 面 “投下 力 疝 下 半 划 为 正 ， 反 之 则 区 为 负 ， 


妇 人 一 


所 以 ， 在 t 时 刻 ， 外 力 位 能 为 
-| vr, ft)u(r, tdr 
这 份 能 景 也 积 暮 在 弹性 体 中 。 可 见 弦 的 总 位 能 应 当 是 弹 几 位 能 与 
外 力 位 能 之 和 。 于 是 
i : 
v=| Hur,t) plr, tyulr, 1]zz 
T=| 二 oz)u3(z,0a 
S$ 一 | a » Pu Ku + wu ]drdt 


= 一 中 Kul+yu = 了 PCz)92 一 Kp'+yu 


由 ”所 .-_,--9-F,=0 得 


”oz at 
BETH)— Pru t+ Ku,,= 
。 __£ wT,t) 
oo. Hi 二 PUL) HH i (31) 


| 一 一 sa 所 一 lL KK 2 (r,t) i 
当 P 振 P(z) 三 上 (常数 ) 时 ， 他 Oy plr, 1) 


ss t+) 及 plz) 均 为 已 知 冰 数 ， 故 P(r,t) 为 已 知 孙 数 ) 
Wit=au,, Pr,t) {32) 

方程 (31》 及 (32) 均 为 纺 的 受 迫 振动 方程。 

例 推导 腊 的 自由 振动 方程 

取 举 标 系 如 图 9 。 wu(z,y， 
+) 表示 在 1 时 刻 腊 上 点 (z， 
y) 的 位 移 。 在 膜 上 任 取 一 小 
块 ， 其 面积 改变 量 为 
AS=VITuitru — 1drdy 


~ + us)drdy 


《 设 已 及 ts 比 1 小 得 多 ) 商 9 


。56 。 
于 是 ， 在 此 小 块 上 弹性 位 能 是 
dU (ust ub)dzdy 〈 玫 是 常数 ) 


,-。 vs) a may 

又 T=||o(r, uidzdy (plz,y) 是 膜 的 质量 面 密度 ) 
i 作用 量 为 

S =| ,dt (ae -Cu 让 Jadzay 


be | [zs ts ul Vt) UT Yt), 
和 
n 


ut yt) drdydt 


这 是 多 元 函数 u(xyy,t) 的 泛 函 ，z,y,t 是 自 变 量 。 它 的 尤 拉 
方程 (见方 程 (18’)) 是 


口 9 0 
ee re 
而 者 1 :_K 2 2 
. F=3P(Y, yu 2 (us+ us) 


a 五 .一 四 ， FP. =—Ku,, F, =~ Kut,, F, =pu, 
” EC )+—0 (Ku )- -9_ (p(x y)u)= 0 
Ox 了 Oy v Ot 9 所 


[4 tt 


K 
Pg ee + uy,) (33) 


当 o(z，2 sp( 常 数 )， 令 入 = q， 于 是 


2 一 GIU Fu) (34) u 
这 就 是 膜 振动 方程 
例 导出 梁 的 微 振动 方程 
取 坐 标 系 如 图 10。 仍 以 u(x,) 表 
示 t 时 刻 梁 上 的 点 7 处 的 位 黎 。 由 弹 


® HT « 


性 力学 可 知 其 位 能 与 曲率 的 平方 成 正比 。 于 是 ， 在 梁 段 [x, zx 十 
ZI1x】 上 ， 位 能 为 


q i ,, . ” 六 
如 UL 一 人 (em) drz~KRui,drd 
U = | 大 


又 T= 二 | P(2)42(x,1)dx (p(x) 是 质量 线 密度 ) 
作用 最为 


S = ( 寺 Pw3- Ku ,jszdi 


-|[ | Fru te dt drdt 
tiJo 


FPF=SPp(r) u(r,t) Kus.(z,t) 


对 于 泛 范 3$， 它 的 尤 拉 方程 是 
9 02 


_.0_An 一 
gy Or F., Ot F. po 


在 第 二 章 中 ， 没 有 讨论 过 这 种 泛 函 的 尤 拉 方程 。 现 在 对 方程 
《35》 作 一 简略 的 论证 。 
把 问题 所 得 更 普遍 一 点 。 考 虑 泛 函 
S=||Fcew ul, y)u zuysd ,Uvv ld ray 


(35) 


在 此 ， 区 域 台 的 边界 9 如 及 函数 下 均 充分 光滑 。 由 6S = 0， 可 得 
泛 函 3 的 尤 拉 方 程 中 。 推 导 过 程 大 致 如 下 : 


_ fff oF aF aF aF 
a5=|[( TR TE i 
LD 


OF 
ee 6s Ysdy 


图 “由 弹性 力学 可 知 ， 常数 公 = 二 吾 了 ， 百 是 场 氏 模 量 ， 了 足 人 的 模 断 面 的 转动 
饥 量 。 


as BB « 


在 此 、#, 二 pp， 2 一 8 HT, Uyy=$ 
在 第 二 章 § 6 中 普 指 出 
jjer, Sp+ Fdqg)drdy 


(3 )Sudzay 
又 
or 9 一 Sr +-37 5s ] Jdxray 
+ oC, ER 
I 0+ 处 况 uj 
ee) Ge en 
9 + 总 3 一 ).5u jazay 
用 格林 公 
由 二 Cs )+-35( dou, ))irdy 
-J Cg 
《4 及 Wy 在 边界 上 给 定 ， .64s| sp 一 0 一 Gzyl ap) 
-六 ( 的 om + 高 (8) J 


”给 定 边 伟人 条件 x| an = 史 ，#:13n = 水 ， wv] an = 名。 在 满足 这 种 边 值 条 件 且 
充分 光滑 的 函数 类 中 研究 证 函 S 的 极 值 。 设 S 在 u(x，y} 上 取 极 值 ， 找 出 
uy {Ts g 引 应 当 满 足 的 尤 拉 方 款 。 


[ls 好 


。 859。 


(名 (次 .各 2 
le 人 六 (全 ea 
-这 (aa 

但 
[i (B50) (er) 


-J( -= 


JP Mar yrs 


-人 这 (3 人 ~( 3 )]svdzab 


oF 
| Rs 3 dy dg 


OO oF ao 3F 
FE “了 Os Judrdy= a 


因为 此 式 对 满足 一 定 条 件 的 函数 类 中 的 任何 Sz 都 成 立 ， 所 
以 
oF D 3F 3 0F, 0 oF 0 OF 


Ou or dp dy dg Dr oa dx: Or 和 Oy: Os 


(36) 


这 就 是 泛 函 S 一 || 亚 (z,gin(z,g)，aevavyacsyavsdzdg 的 万 拉 


方程 《此 方程 用 得 不 多 ， 在 第 二 章 中 未 提 到 ) 。 
由 此 可 见 ， 对 于 泛 函 


9 60 » 


t . 
s=|'*| P(r,t ur, t) ,uu uu, drdt 
1J0 


ee a 


品 7 9 、 
求 出 五 ，， i Ss 代入 (35) 式 ， 得 
-$0 (37) 


这 就 是 梁 的 自由 〈 横 ) 振动 方程 。 

例 求 在 重力 场 中 定妆 的 轴线 的 7 
形状 〈 设 染 的 两 端 固定 且 轴 线 的 弯曲 
程度 很 小 ) ,注意 现在 梁 受 重力 作用 ， 
所 以 不 能 用 梁 的 自由 振动 方程 。 又 定 
梁 处 于 平衡 状态 ， 这 时 动能 为 零 。 取 
坐标 系 如 图 11。 以 kz) 表示 梁 轴 线 
上 z 点 处 的 位 移 〈 现 在 * 与 时 间 无 
关 ) 。 在 [z,z+-xz] 谍 的 重力 位 能 = -9p(r)a(z)dz。 此 鞭 


的 弹性 位 能 = 广 已 7z3-dz。 勾 益 位 能 是 二 者 之 和 。 所 以 


图 11 


v=| (3 — Eu — P(r) u(r) Jaz 
作用 量 泛 函 为 
Ss=-|,| ($87u2.— gp(r)ur) Arat 
eh i ey zs— GP(T)u z)] Tz 
F=3EJui.(r) + gpP(x)u(r) 


尤 拉 方 程 龙 


a 人 1 。 


- gpP(z)— FIJIud(r)=0 
#9 (x) = 2 (38) 


这 是 4 阶 常 微分 方程 ， 它 的 积分 曲线 有 4 个 任意 常数 ， 它 们 可 由 
边界 条 件 定 出 。 由 此 可 得 梁 的 轴线 方程 。 

这 个 例子 也 可 用 如 下 方法 求解 。 当 梁 处 在 平衡 态 〈 稳 定 平 
衔 ) 时 ， 动 能 为 等 ， 位 能 最 小 ， 所 以 


U=| (E71, gO(r)u(r) Jaz 


-| PFIix,u(r), uw (rdzr 
0 


取 极 小 值 。 

这 时 UCu(z)] 是 一 元 函数 4u(7) 的 泛 洲 。 若 可 在 u(x) 上 取 极 
信 ， 则 u(x) 满足 尤 拉 方程 (12》 ， 
写 出 的 方程 仍然 和 刚才 的 相同 。 

例 设 梁 的 两 端 因 定 ， 在 粱 上 某 
点 有 一 集中 荷载 。 试用 广义 坐标 求 染 
的 挠 度 。 

取 坐 标 系 如 图 12 所 示 。 设 梁 上 坐 
标 为 C 处 有 一 集中 力 已 作用 。 这 时 梁 


图 12 
轴线 可 表 为 富里 哀 (Fourier) 级 数 
y= gine + ausin E+ sin -2 


也 就 是 说 y= 了 (al，as， ……，4,.，…)。 
现在 把 〈cli，a:，…，as …) 视 为 此 质点 系 〈 即 梁 ) 的 广义 坐 
标 。 于 大 此 系 有 无 穷 个 广义 坐标 。 过 去 我 们 只 讨论 过 具有 有 限 个 
广义 举 标 的 质点 系 ， 但 是 现在 我 们 不 妨 继续 讨论 下 去 。 
C 点 处 的 位 移 为 
.~ 27C 


y(c)=asin e+ oa,sin 


I! I 


十 Se 十 a.,Sin 2 二 人 全 四 


# 62 » 
外 力 五 产生 的 位 能 为 个 


nixre 


U,=-Py(c)=-Poa. A 


tw] 


弹性 位 能 是 


Us=3E1 | 


= 方 EJ | D> 一 一 CSin 3 qz 


4 
二 Wan 


nm 


是 总 位 能 为 


S a 
U=U,+U,=— PHasin +E Sos 


ml i 


一 Ca CE 7) 
因为 粱 处 在 稳定 平衡 态 ， 所 以 忆 取 极 小 值 ， 于 是 


A 0 (n=1,2,..) 
. a MTC 
oe D173 nia.~— Psin 7 二 省 
A nae 
人 “1 f 
由 此 得 轴线 方程 


中 可 套 铝 徐 芝 给 所 编 《 弹 址 力学 简明 教 窟 中 的 人 $5 一 ? 

Gf mm 夺 RR 

名 ”在 这 自用 到 fF Sin Tsin TAIT dr = ti 
9 i ! py 


当下 于 好 


ss 63 。 


“I 
nnec 好 区位 
y= SET 2 这 Sin ~ 一 一 Sin 7 


由 此 方程 可 见 ， 当 c = 0，! 时 ， 上 4 三 0。 这 就 是 说 ， 当 集中 力 
作用 在 端点 ， 梁 不 发 生 弯 曲 。 
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第 四 章 ” 变 分 法 的 新 发 展 一 一 有 限 元 法 


在 第 三 章 中 我 们 举 了 许多 例子 ， 对 于 它们 我 们 都 是 写 出 有 关 
泛 重 的 尤 拉 方 程 就 止步 。 但 是 在 解决 实际 问题 时 ， 还 要 求 把 这 些 
方程 解 出 来 。 令 人 遗憾 的 是 这 些微 分 方程 多 数 是 不 能 积分 成 初等 
函数 的 。 于 是 人 们 只 好 另 想 办 法 。 所 谓 “ 解 变 分 问题 的 直接 法 ” 
就 是 适应 这 种 需要 而 产生 的 。 尤 拉 开 始 研究 变 分 问题 时 ， 就 用 过 
“有 限 差 分 法 ” ， 这 是 “直接 法 ”中 的 一 种 重要 方法 。 所 以 尤 拉 
也 是 变 分 问题 中 站 按 法 的 创始 人 。 后 来 又 出 现 了 里 兹 (Ritz) 法 ， 
个 过 金 (I axrepxxH) 法 等 重要 方法 。 在 这 些 方法 的 基础 上 ， 产 生 
了 “有限 元 法 ”。 这 种 方法 最 早 见 于 柯 朗 (Courant)1943 年 的 论 
文 。 他 在 这 征文 章 中 ， 首 次 用 了 一 组 “三 角形 元 素 ” 来 解决 力学 
中 的 一 个 握 转 问题 。 但 这 种 思想 当时 并 未 引起 重视 。 在 柯 朗 这 一 
论文 发 表 后 大 约 十 年 ， 才 有 人 做 类 伺 的 工作 。 许 多 人 (包括 工程 
师 和 物理 学 家 ) 在 不 同 的 领域 中 独立 地 发 展 了 有 限 元 法 中 。 到 了 
60 年 代 ， 有 限 元 法 正式 登台 @。 时 至 今日 ， 它 的 应 用 范围 越 来 越 
广泛 ， 理 论 上 也 还 在 发 展 。 对 今天 的 工程 技术 人 员 来 说 。“ 有 限 
元 法 ”已 成 为 必 备 的 工具 了 。 在 本 章 中 我 们 先 扼要 地 谈 一 下 里 兹 
法 〈 从 数学 关 度 看 里 效法 与 有 限 元 法 的 指导 思想 极为 相 侯 》， 然 
后 介绍 有 限 元 法 。 


@ 在 有 限 元 革 创 始 人 的 名 单 中 ,还 应 当 提 到 工程 师 Argyris, Turner, Martin 
及 Topp 等 人 人， 还 可 提 到 物理 学 家 Prager，Synge 及 MeHahon。 这 就 像 一 种 
加 实在 许多 果园 中 同时 成 熬 一 样 ， 罩 要 的 数学 思想 也 有 它 的 成 熟 生 节 ， 

多 ”中 国 数学 家 冯 康 在 1965 年 首次 证 明了 “ 线 姓 元 2? 的 收 艇 伺 。 这 是 当时 的 世界 
记录 | 


9 65 。 


$1 从 时 次 法 谈 起 


1. 从 具体 问题 说 起 


对 于 下 述 的 常 微 分 方程 边 值 问题 
了 工 1 一 -2 (Piz) )+ rr)ru(r)~ f(r) 
(or<b) (39) 


u(a)= 0 =u(b) 
在 此 ， 设 ，P{x)eECira,by, r(xz)eCra,b] 
Prfz)>>0，riz)>0 凡 rz) 不 恒 为 沉 。 
我 们 把 “ 求 函 数 w(z) E cz[a,53 昌 满足 (39)* 称 为 解 问题 (39)。 
取 泛 函 
Qru(z)]=| CP Cx)) uCr) + rz) u(r) 
a —2f(r)u(r) dr (C40) 
当 u(x)E cfa, 65] 时， 显然 Qfu(z)J 有 意义 。 
现在 我 们 来 证 明 ， 若 4,(z) 是 问题 (39) 的 解 ， 则 wkz) 使 活 
计 外 [Cu(z)] 取 最 小 信 Q@[uotz)]。 皮 之 ， 车 W(xz)Eci[a, 站 使 
泛 函 QCu(z)] 取 最 小 值 QUuo(z)]， 则 xfz) 是 问题 (39) 的 解 。@ 
证 ，1)》 设 xoCz) 是 (39) 的 解 。 任 取 环 (z)ECira,6]， 显 
然 六 二 ut+WECS(a,b), 上 且 
QV Ca)= {CCPV Yr 4 ra Id 
=—PVV| + Pr rr 2 dz 
但 Fr(o)=0=Frb ”一 
所 以 
QV Cz)I={ CPV + rt 一 27J]dz C41) 
(显然 当 Y (XT)EC8Ca,b]， 《41》 式 就 成 立 ) 又 


@ “将 所 有 在 fa，b】 上 有 二 阶 连 续 导 数 且 在 冰 端 点 为 零 的 函数 的 集合 记 作 
Cota, b), 
见 Céras bY= {fw f(r EEC:a, by foam 0=f(b)) 。 


。66 。 
QrFr(z)]=Orus(ry+JP(z)] 
={ CP t WY tr ut WY) 
au + 了 az 
={ CPusttrus—2fudz 


十 2 [ce ui +trul dx 
-2 [fwaz 
+ | CPW’ + rv de 


但 是 
{cpusm + ru dr 


= PuW [+ {f= Put)’W+truMW dr 
={ Cc-(Pus) +ruoliy dz 
=| waz 
所 以 
QLVCX)) 


=f CPuts+ 9 {cPw +r]de 


= QCuoz)I+ CPW’ +rmWijdr 
即 


Oro Fw is Or + {CPW :+t rm de (42) 


显然 ， (42) 式 对 任何 夯 (z)EC 红 ,53 成 立 。 和 但 是 著 玩 (z) 不 
恒 为 零 ， 则 


站 67 时 
{EP tz) rr) dr> 0 


实际 上 若 ) LP(z)i cz)+ rz)]az 一 0， 则 由 
POC)W' (zr)+r(z)W2(z) 是 [a,6J 上 的 非 负 连 续 函 数 ， 司 知 
P(r (rT)+rr) (rz) 0 
于 是 PLX ’':(7)= 0 
但 PLrz)>0 在 [a,6J] 成 立 所 以 
W(x)=0 
于 是 WV'z)=0 
所 以 
WW(z)=C (常数 ) 
又 Vo)=0=W(6) 可 见 人 W(xz)=0 
这 就 是 说 ， 只 要 搬 (z) 不 但 等 于 useoGz)， 由 


QV (r= QLuo(r) +) CPW + ri dx >Q[u (rz)) 


可 见 泛 函 QLx(Cz7] 在 (39) 的 解 zo(z) 上 达到 最 小 值 Qiue(Cz)。 
2》 设 泛 隙 QLulz)J 在 u(xz)ECiCe,b] 上 达到 最 小 值 
Qtz)j。 现 在 ， 任 取 凡 tz)ECi[c,p]， 对 任何 实数 4 来 说 是， 
显然 有 
CW[ao(z) 十 4 了 《zz) 


=WLuo(r))+ | PW 六 十 rd 


十 2 | CPousCaW’) + ru a) 


— f(AW )Jdr QL Wor) 
当 zo(z) 及 下 (zz) 加 定时 ， 


ftP Ow’ +r(Am ydz +2) CPousCAW’) 
truoCAW )— fi ) dr 


中 指 对 任何 a Cz) & Cicayb] 来 说 ， 有 Qtuwkz)] 艺 Quocz)] 


. 68 [3 
可 视 为 4 的 国 数 ， 将 它 沁 为 PC4)， 于 是 
=) PW + ri dz 


i 24| CPousW’ ri 了 一 Jazx 


>0 
对 于 二 次 三 项 式 A4:+ BATC 来 说 ， 淳 
Ai:+BATC20 
对 任何 实数 成 立 ， 则 28: 一 44C 志 0。 这 时 ， 若 C= 0 则 8B:<0， 
所 以 8 二 0 。 册 此 可 部 
J Pout + rudy ~ fy dx=0 
所 以 
ECPout)’ + ru I dz =00 
对 任何 大 (z)EC 红 Le,b] 成 立 。 由 此 可 以 证 明 儿 
一 oa) +ru -f=0 
于 是 xue(z) 是 问题 (39) 的 解 。 
事实 上 ， 阁 在 某 一 点 tS (a,6)， 有 
— Prudlzo)) rz) ro) -fro)>0 
我 们 就 取 一 个 函数 7(z)， 它 在 避 点 附近 大 十 零 ， 症 其 余 位 置 上 
为 零 ， 且 7CX)E CiL4,6J， 于 是 
.ECPerYutCe) rr) uz) /2)I ndz>0 
这 就 引起 矛 征 。 辐 样 ， 若 
— PTOu To)) —rtroOu (tro) — zr) <L0 
也 导致 证 盾 。 可 多 


@ 在 这 里 出 到 人 poufW' dz= PuutW -fePouy’ Was 
= -全 《人 
儿 参 着 本 书 第 二 章 $ 2 之 2 中 的 “推广 的 预备 定理 。 


s 89° 


—(P(r)us(r)) —r (ru r) f(r) 0 

这 就 是 说 ，4o(f5) 是 {39) 的 解 。 

我 们 证 谓 了 ，wo《z) 是 (39) 的 解 的 充分 且 必 要 的 条 件 是 : 泛 
溪 Q[Cwutz)] 在 u(x) 上 取 最 小 值 。 

这 就 是 说 ， 求 uu(z)EC3Cta,bJ]， 使 汉 函 (40) 取 最 小 值 与 解 
问题 (39) 是 等 价 的 。 

由 于 这 个 发 现 ,事情 发 生 了 逆转 ,刚才 提出 的 任务 是 解 (39)， 
现在 我 们 掉 过 头 来 求 u,(X)， 使 泛 函 QCuolz)] 取 最 小 值 。 

里 次 法 是 解决 这 个 问题 的 一 种 近似 方法 ， 大 致 如 下 ， 

为 求 u(xz)EC3Ca,b]， 使 泛 函 
QCu(z)I=f CP Ce) ur) tr Cr) uz) 
—2f(r)u(zx) lar 

取 最 小 值 ， 我 们 选取 基 范 数 
PLZ), PT PAT)Y piCZJEC8Co，Do] 
(i=1,2,.",7) 


设 它们 是 线性 无 关 的 。 令 
Sap(r)=u,(7)D 


f=l 


Ql To ,9 为 任意 实数 .于 是 ,省 画 C[us(Cz)] 成 了 alias a, 
的 函数 。 
为 简化 书写 方式 ， 记 积分 


人 fa)etzyez=(F9， 
把 它 称 为 f 和 9 这 两 个 函数 的 内 积 。 于 是 
QCuz y= C-(P HW )' utriu—2fuldr 


名 这 就 是 说 ,和 若 a,9 1 (CT) + 898(T) +. 玫 ,(Z) 宇 0 当 z€Co, 了 DJ 时 成 立 ， 
则 常数 a 二 [中 (f=1,2,3R) 
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一 (了 1 四 一 2( 太 0) 
由 上 4 的 定义 可 知 
LCbidt ot) =m b+ 6.Lu 对 任何 数 b1.8, 成 立 。 因 
此 ， 我 们 说 算 子 工 是 线性 的 。 又 
(f ,bg + b.9)=b.(f,9)+ b,(f, 9) 
(bfit+bf,,9)=b,(f1,g9)+ 6,.(fs,9) 
所 雇 
Qu (r= Lu, uu,)—2(f.u,) 


=(L( 2a.9.), Dap)) 


~ 2(f, a,91) 


= > a,a(L9,,P;) —2 Ya,(f,,9,) 


dd 
tf» f=l 了 上 了 


一 天 (aiq ga) 
在 此 ， 刻 (a,,…,4,) 表 示 QR[Ew,(7)] 是 01,…a, 诸 变 景 的 函数 。 
现在 选取 21,42,… ,a 使 QCw.]== 耻 (a1,…,4.) 取 最 小 值 。 
用 数学 分 析 中 的 老 办 法 求 极 值 ， 即 由 


9 =0 (F=1,2,.,n) 


得 之 ,201( 工 po 站 一 201 ,91)D 


了 1 


名 上 uw 的 讲义 见 (39) 式 。 


@ 注意 > qrgqtLP P= oo LP PI aa LP, Pp) + 
ds- 
+ agtLP PP Yt a PP +arar (LP ,Ps) 
中 we 十 dsao, (LP PI tmta di (LP ,Py 
+ gr0 LP oT) +t arda (LP,pP,) 


重 ?1 四 
所 以 DaLo PH)= ,0) (f=1,.,n) 
(43) 
这 是 关于 a, ,4s,… ,a, 的 线性 方程 组 。 解 出 (43)， 得 到 ab 
42,* ,0, 的 值 ， 也 就 得 到 as.(z)= > api(z)。 我 们 利用 


ed 


(Los sp = CP) p+ rp, pd 
[3 由 
=-—Poip,| + CPY!P! + rpp, dz 


b 
-| CPPIPDI+roP,ldr 


(注意 ¥ (0)= 0 =29,(45)) 


再 记 
{tpPoiws +r9,pdr=B(,,9)) = hs 
又 令 
(f ;Pi) 一 b, 
于 是 方程 组 (43) 可 表 为 


Ei 


| Rl1G1+ Rady+ + hi 9 一 五 | 
| RoiG1 ft kids + 十 Ri.4, = bs 

(44) 
| RL a, 盾 及 .Gas 十 ov kia0,=b, 


或 用 和 矩阵 记号 ， 写 成 
£,A=B 
企 此 
K =[k,/) 
4 一 [aiqa,…y， as] 中 7 一 [D ,bb,.ID 


他 ”47 表示 所 阵 过 的 转 时 矩阵 。 为 便于 排版 ， 才 这 样 书写 。 


。72 全 
现在 证 明 甜 阵 久 是 对 称 正定 的 ， 于 是 方程 (44) 存 在 唯一 解 初学 
者 可 以 暂时 不 看 以 下 几 行 文字 。 娄 搞 清 楚 的 话 ， 请 参 出 线性 代数 
课本 。 
因为 
Re 一 人 5Peipf+ rp de = 
所 以 抵 阵 玉 是 对 称 的 。 青 证 明 抑 阵 玉 是 正定 的 ， 即 当 
CC7 一 (ciicaca) 办 0 
了 时， 
CTKC> 0 (注意 CTKC 是 一 个 实数 ) 
事实 上 
(es 本 里 
CIKC=Fcic[RD | : l= > oe, > hiio 


= > ocikis= > ciciB(P, Pp) 


Fn 


了 S eic] CPoip! +t ryder 
1s» 了 一 上 
-fp( co) (Dew) 
下 j=l 


十 了 (5 op,)( Dep Idr 


tml 1+:1=i 


-a(Be0, Sov) 
fi 


1 


=BP(u (7Z) ur)) 
在 此 处 


加 在 此 处 0 表示 零 和 打量 (0 0,",0》。 
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HC 雹 ) 二 pO A 


1~=| 
当 5 严 0 〔 即 cc 不 全 为 0) 时 ， 这 时 
#.《 工 ) 不 代为 零 D， 所 以 只 要 能 证 明 有 5 之 0 存在 ， 使 
crKc={ CP + ru ) dz>6) wid 


由 由 
| Be 
订 知 
CrKCS>0 
故 轨 阵 氏 正 定 。 


现在 来 证 这 种 正 数 5 确实 存在 。 因 为 

人 (Dae= er) (注意 xE Ci[e,6] wu(a)=0) 
又 

pcs as|<¥Y foc ade yf pyar G5) 
所 以 


fo ux) dz 一 fu ddr < 


I 
{Fh EAE 
(fra) 


= dt 


po—a): 


3 2rd 


DD 因为 加 1 (2) ,m+, 恩 ,。(X} 线 性 无 关 。 
网 通 尚 招 册 不 每 式 称 为 许 吕 人 尔 兹 不 等 式 。 请 参 疙 斯 米尔 诺 夫 所 著 《 高 等 数学 孝 
程 》 第 五 卷 第 一 分 册 第 二 章 第 55 节 定理 3 。 


于 是 
| za sy), rd 
ee 
这 T min Plrx)u’ rz 
a KELar bb) 
5 
>| min Pz) 人 Et ru2]adx 
>| [minPkz) jidz 
=6| ez)dz 
在 此 


6= mim Re 


zeEra — 9 


由 年 阵 天 正定 ， 可 方太 4) 有 唯一 姑 (09,。 a9)。 
再 将 函数 (ai,… ,4,)= S a. aj 


19 J=l t=l 
求 导 两 次 ， 得 
oF 
Qaida 
可 见 完 阵 
[ OF 
Qa.00) 
亦 正 定 ， 所 以 函数 玉 (a1,… ,a,) 确 实在 (a?,…,a9) 处 取 极 小 值 。 
但 如 (al …，aw) 是 二 次 函数 ， 所 以 此 极 小 值 也 就 是 严 的 最 小 值 。 
当然 ， 当 x 取 定 后 ， 泛 函 QEu,(x)} 的 最 小 值 不 一 定 就 是 泛 
闭 Q[u(z)]j 当 w(x)ECi[ae,6] 时 的 最 小 值 .但 是 ， 在 数学 上 证 明 
了 一 个 重要 结论 ， 就 是 ， 对 于 每 个 n， 我们 用 上 述 办 法 (了 即 里 兹 
法 ) 求 出 (a?,…,03)， 得 到 


=2k,} 


和 


把 它们 称 为 “ 变 分 问题 的 近似 解 2> 中 ， 把 这 种 

HT) HL) ,aT) ,ee 
序列 称 为 这 一 变 分 问题 的 “ 极 人 小 化 序列 ”， 记 此 变 分 问题 的 精确 
解 为 so(z)。 有 


fs) ur) dz 0 (Mn) 


也 就 是 说 ， 由 里 兹 法 得 到 的 极 相 化 序 济 {x,(x)} 平 方 平均 收 襄 于 
变 分 问题 的 精 痛 解 zxo(z)。 
这 样 ， 我 们 用 里 效法 解决 问题 才 有 了 稳固 的 理论 基础 。 
我 们 在 理论 上 稍微 走 和 下 了 一 点 。 虽 然 我 们 已 经 声明 初学 者 在 
读 第 一 遍 时 可 以 不 理会 这 些 ， 但 总 是 一 种 心理 上 的 负担 。 为 了 组 
和 紧张 气氛 ， 还 是 来 看 一 个 例子 吧 。 
例 求 函 数 x(z)EC50,1]， 使 它 满足 
人 rE(0,1) 
u(t0)= 0 =2(1) 
这 方程 可 在 (39) 中 令 记 (7) 圭 1 二 r(x2)，f(z) 寺 一 % 而 得。 此 方 
程 的 满足 边 值 条 件 at0)=0==u(1) 的 精确 解 是 
ux) (ee ")—z® 
现 柚 用 里 兹 法 求 近似 和解 ， 以 资 比较 。 
取 P (x)=2tl 一 7x), PTI)=r(1 ~ 2), PI)=I(1 一) 
显然 ,0)= 0 =%(1) 
且 Y,4z) 在 [gj 有 二 阶 连续 导数 ， 也 就 是 说 
PT)ECILa 6b] Ci i,2,3) 


@ ”我们 把 “ 求 z(z) E C 拭 ac， 使 泛 图 Qutz) 1 双 晤 小 值 ”了 称 为 《 变 分 间 
题 > 。 


名 此 淫 可 用 Lagrange 变 易 常 是 法 求 和 但。 请 参 纲 营 油 分 方程 教 本 。 
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还 可 以 验证 91,8,,P; 是 线性 无 关 的 DP。 所 以 它们 可 以 选 为 基 范 
数 。 
用 


Ai 一 人 cei(z)p1(z) 二 mi(z)Pi(z)]az 


(7 一 1,2,3) 
可 得 
Ri 一 0.3667 kl, =0.1833 kis=—0.1763 
ksi = Bi, Rs =0.1338 Rss= —0.8940 
hsl= Rls Rs, = R,s Rss=0.0897 〈 取 近似 值 》 
又 


人 1 1 
2 Sodz={ -= 《1 一 zj)dz 一 一 了 


1 1 
ba={ Spdzr=) sd = -站 


1 l 
| /vdr=| ER 
0 


30 
由 此 得 到 线性 方程 组 
0.3667a,+0,1833a: 一 0.1763as 一 -十 


f 
| 0.1833a, 二 0.1333a, 一 0.8940a, 一 -而 


( 
Ee -区 
EE, 
qi 一 一 0.4425，ay 一 0.4635，as 一 0.0343 
由 此 得 到 近似 解 
Us(T)—=—0.44257(1— 7X)+0.46357x°(1 ~ x) 


+0.03432°(1— xr) 


全 着 CP,+CP:+ Cpe=0, 则 -Csxrtt (Cs -Cr+ (Cr — Cr 
+ CT 二 4 于 是 ~Ca=6O=C-C=C1~ C=C, MiC= 0=C,=C, 


7 。 
一 一 0.034374 一 0.429273 十 0.9060zz 一 0.44257 
试 选 点 == 亏 ， 算 出 精确 解 u(z) 在 z 一 二 处 的 值 4( 】)， 又 算出 


ui( 到)， 作 一 比较 ， 以 观察 近似 解 误差 的 大 小 。 


一 一 0.0566 


os( 三 ) = .0343(2) —0,， 4292( 二 ) 


一 一 0.0505 
1 2 3 4 5 6 
现在 对 于 点 ZT 一 8， 8， 祷 ， Bt 8， 8 8， 分 别 求 出 近 
似 值 wa Zz ) 及 准确 什 u(z) 呈 并 列表 如 下 ， 
表 


准确 信 Wiz) 近似 值 ws(z) 


~0.0191 E ~™0.0420 


机 


一 0.0350 -0.0609 


3 
中 


~0.0484 ~0.0618 


nH 
hh 


~0.0566 | -0.0505 


3 


~0,0579 -0.0327 


我 
tt 
= = |> “| “| “| 


一 0.0503 一 人 .0142 


时 
4 


H 
1 
| 


my 


-0.0317 ~0,0017 


看 来 误差 还 不 算 小 ， 但 是 要 注意 ， 我 们 只 取 了 三 个 基 画 数 (=3) 
阿 ! 
加 ”这 里 的 准确 值 当然 有 会 人 误差。 


« T& » 


2， 里 兹 法 的 计算 步 了 又 

上 面 已 经 详细 地 讲 了 一 个 特例 ， 现 在 把 里 效法 的 计算 步 又 总 
线 如 下 。 

设 M 是 一 个 函数 类 。 函 数 类 MM 中 的 任何 活 数 u(7) 均 满足 给 
定 的 边 值 条 件 和 一 定 的 光滑 性 《在 上 文中 ， 上 由 二 C3[a,4b]。 在 各 
种 问题 中 ， 术 可 取 成 各 种 类 型 》 。 要 求 出 uolz)EM， 使 泛 汉 
Qtia(z)3 取 最 小 值 。 这 是 一 个 变 分 问题 。 

我 们 用 里 效法 求 此 变 分 问题 的 近似 解 。 步 又 是 ， 

1)》 在 函数 类 好 中 ， 选 取 线 性 无 关 的 几 个 基因 数 p1 ,Pp;,…， 
wp.。 作 


ty Sarpy) 


记 M,=! Tu.p,(z) | 了 一 切实 数 ! 
这 冉 泛 函 
Qtu(z=Q[ Dar.lr) ] 


= (A,0s,", G0,) 
也 就 是 说 ， 它 是 变 元 alias,… ,au 的 2 元 本 数 。 
2) 求 出 ay ,a at, 使 洱 数 C4,,… ,a,) 在 《2 ,a9,. 
4s》 处 取 最 小 值 ， 即 
严 (ag al，Gy = min QLu(tzr)ld 
u {x) EAM 
oF 
办 法 是 取 oa = 0 
得 到 线性 方程 组 (44)， 即 
天 ,过 一 如 
解 此 方程 组 ， 求 出 《GD a2, 0 9 在 一 定 条 件 下 ， HZ ) 一 


也 显然 mia QuacTl] 之 min Qtutz)D] 
区 us) EN 


» TI rr 


之 a49,(z) 使 玫 函 QTu(z)7] 在 M。 中 取 最 小 值 


CGfa(z) 一 min. QIru(r)J, 且 


{ 2) — wr) dr 0 (Bn 00) 


CuolI) 是 变 分 问题 的 精确 解 》 
这 就 是 里 兹 法 的 基本 轮廓 。 


$2 有 限 元 法 


在 用 里 赣 法 求 泛 函 QEu(z)] 的 极 入 前 数 Duo(z) 的 近似 解 时 ， 
我 们 先 选取 试探 孙 数 P,(x),… ,P(x)《 即 上 文 所 说 的 基 菠 数 )， 


求 x.0z)= > 20.9.(7)， 使 泛 孙 QCu,(z)] 在 试探 函数 类 M, 上 


取景 小 值 。 把 这 个 wu,(z) 作 为 x,(z) 的 近似 解 。 这 个 方法 的 指导 
思想 是 容易 理解 的 。 一 般 说 来 ， 也 收敛 得 比较 快 〈 指 ?不 太 大 时 ， 
误差 比较 小 》 。 但 是 ， 里 效法 对 试探 函数 P,(Cz) 的 光滑 性 要 求 较 
高 。 对 于 边界 曲线 比较 复杂 的 二 维 区 域 (高 于 二 维 的 区 域 更 不 竺 
言 》， 要 构造 出 这 种 试探 性 函数 往往 是 困难 的 ， 有 财 还 是 不 可 能 
的 。 此 外 当 半 增 大 时 ， 精 度 提高 很 慢 〈 注 意 这 上 与 上 文中 所 说 的 
“ 当 # 不 太 大 时 误差 比较 小 ”并 不 耶 盾 》 。 面 对 这 种 情况 ， 有 些 
人 对 P(x) 的 光滑 性 降低 要 求 ， 干 脆 选 肥 分 片 多 项 式 《〈 见 下 文 ) 
作为 2,(r)。 这 种 基 国 数 仅 仅 连 续 ， 在 有 些 点 连 导数 也 没有 ， 这 
就 扩大 了 试探 函数 类 的 范围 。 当 时 ， 人 们 难以 估计 到 的 是 ， 这 种 
做 法 居然 引起 了 计算 方法 的 但 火 变革 ! 产生 了 现代 的 适用 于 广 阅 
领域 的 有 限 元 法 。 本 节 打 算 从 变 分 法 的 角度 介绍 一 下 有 限 元 法 的 
指导 思想 。 

再 来 看 一 下 问题 〈39) 。 我 们 已 经 知道 ， 只 要 求 出 u(x)， 


虽 ” 指 泛 诡 Qtw tz) 在 we ix) 上 取 梓 值 。 


使 〈40)》 式 所 表示 的 泛 画 QCu(z73 取 最 小 值 ， 则 wz) 是 (39) 
的 解 。 但 泛 函 QCv(z)] 可 按 下 式 进行 变换 : 
Qru(nD3=| 5-CP(z)a Cn)) az) 
十 8 一 2 (rT)u(rT) dr 
=— P(x)u’(r)u( rz) 


+| CP Cs) usr) +raa(z) —2f(r) u(r) dr 


=| CPu'ttru-2fulde 


(注意 4#(a)=0=4u(6)) 
现在 ， 我 们 用 有 限 元 法 来 求 u,(x)， 使 泛 通 


QCuCz)=H| Pu str 2fudr® (46) 


在 xz) 上 到 最 小 值 。 
1 》 选取 分 片 多 项 式 作 为 基 函 数 〈 这 是 关键 问题 ) 。 我 们 要 
作出 一 个 由 线段 连 按 而 成 的 函数 (7), 使 它 在 Ce,8] 上 满足 
V(r,)=u(7T,)=u, (=0,1,2,.",7) 
《 取 a 一 ro 一 ZW) 
《 见 图 13》 
{zi—r)fh MrECro, T) 
取 Pol zx)= 4 和 入 
(0 当 z ECLzxo rz) 
(rT—7T0)fh 当 TE[TXo, Ti] 
Pz)= Ge/ 有 当 z EL[zI, Tz: 
0 当 工 后 [zs I 


”为 便于 推导 有 关公 式 ,我 们 在 原来 的 Qiu(z) ] 表 达 式 前 乘 以 过。 若 wukz) 使 
人 Pa+ ru- 373dz 取 最 小 什 , 则 utz) 也 使 二 人 5Puc 十 后生 一 
2fujdzx 取 最 小 重 。 反 过 来 ， 有 类 似 结 论 。 


s 员 1 ， 


(zz 一 DA TE[TIL1 TN) 
PF) mT hr ECT Tn) 


0 TELT TU LT ty) 
及 


we 当 xzE[z,-i 7) 
人 (Cr 一 lh, MTELT 1 Cn) 
Vo(z) ,P(X) 及 V(X) 分 别 如 图 135、c、d 所 未。 


Us VX) 


< 

| 一 一 

人 

| 3 
| Ea 

| a 


Xs Yi 
| 
1 | 
上 
9 (A) 
A 
SR | 
Xo Ns Xz 1 
1 l 
: | 
: 1 
4 
| 1 
t J 
1 gz) 1 
[7 
Trt Ny 


sa RI 。 
将 图 13a、b、c 、d 联系 起 来 观察 ， 可 以 看 出 


V(x)= Sup,(7) 
j= 


刚才 选取 的 Po(T)，P1(F)，…，P,(7z) 都 是 分 片 多 项 式 的 最 简单 
情形 《中 分 片 线 性 多 项 式 ) ， 在 有 限 元 法 中 ， 这 是 很 重要 的 基 函 
数 ! 现在 将 有 限 元 法 文献 中 关于 %,(z) 等 的 适用 术语 介绍 一 下 : 
区 间 (zx;.1,7X,) 称 “单元 ”, 记 作 8;。Zo，X1，"…，X, 都 
叫 “节点 2 。 又 定义 如 下 函数 
NIN(r)=(T1— x) /hi 


(47) 
六 DZ 一 (人 一 2 ) fh 
把 它们 都 叫做 “ 形 画 数 ”。 于 是 ， 上 述 共 函数 可 表 为 ， 
(六 ‘(rz) rEe, 
P(r)= < (48) 
‘0 当 X Ee 
又 当 0< 所 1#， 有 
NIN(r) 当 r Ee, 
PAT)=E NIN(r) WrE en (49) 
“0 当 z EejyU en 
又 
NI (z) 当 z Ee, 
P(r) = (50) 


0 当 z €e, 


这 些 式 子 虽然 难民 些 ， 但 只 要 用 得 多 了 ， 也 就 部 了 。 
家 定 z 一 疡 0， 其 余 #，2，…，2. 任 意 变 从 ， 风 项 数 斑 (z) 


= > uv2i(Z) 的 全 体 组 成 试探 函数 类 1,。 用 集合 记号 表示 ， 就 
是 
MM 各 wn.(7) | = Pa 取 任意 实数 


se 总 3S 。 


任 取 (I)E 寻 .,， 显 然 (q)= 二 全,。 又 因为 
0 当 i1 守 72 


所 以 


V(x,) = > up (r= uP (rT) = 
这 样 ， 我 们 已 经 用 节点 
ee 
将 区 间 《aebD) 划分 为 2 个 单元 2 一 《ZJ 7) 《7 一 1 1》 ， 
选 定 了 基 函 数 {po(z)，9Pz)，…，9waCz)}， 并 构造 了 试探 丫 数 


(Zz)= CRASS 


2 ) 按 单元 写 出 单元 护 阵 
和 里 效法 一 样 ， 下 面 的 工作 是 要 在 试探 函数 类 好 .中 ， 求 出 


一 个 函数 xs(z)= > wi9,(z)， 使 泛 了 QLV(z)I 在 这 个 wu(Cz) 


上 取 最 小 值 ， 即 
min QV (7))= QLu.(z) 


Vx} 下 
和 


当 信 (7)E 以 ,， 这 时 


Qu (zy)=0 unr) ]= Fmt,) 


i=0 


注意 到 uo 二 了 ,已 经 事先 固定 ， 所 以 当 (z) 在 于 ,类 中 变化 时 , 泛 


辐 由 此 可 以 证 明 P9 ,tz)， PCZ) ,Pp,(z) 线性 无 尖 。 事 实 上 若 》 cP,(7) 


i= 人 0 


0 ， 取 zz=zi， 则 》 cpi(zi)=ci=0(fe0il ,1) 
了 和 


ee 共和 


函 QCV CX) 是 wi，w1，*…，U, 的 mn 元 函数 。 凡 要 再 用 (这 也 和 里 
兹 法 一 样 ) 
oF 
On, 
写 出 关于 ,ws，… ,4w, 的 线性 方程 组 ， 再 解 出 u?, x2,，… ,wu?，, 就 得 
到 近似 解 


一 0 (1 =1,2,°,7n) (51) 


au(z)= > up,(z) (ud = = PP,) 
但 在 单元 。; 较 多 时 ， 为 使 计算 工作 有 条 不 率 ， 人 们 总 是 先 按 单元 
写 出 单元 矩阵 ， 再 组 合成 总 矩阵 ， 才 得 到 方程 组 〈51) 的 。 上 基体 
说 来 ， 因 为 


QC- CPV str afr ds 
- 1 全 7 本 
= 之 工 | CPV ?+rV 2 2fV dr 
ji i-1 


-SH Pr tar ds = Do 


j= i=l 


在 引 ， 
Q,= 工 | CP pp a 
2 * 
注意 yls)= Suw, (2% A 
J-0 /1 1 ~ 
在 ej 单元 ，zj-<IZ<Y， 这 时 只 gw” 1 1 < 
有 Zi;(x》 和 8;-,(T)》 不 为 零 ( 见 Nj- FE Ti 
图 14) 。 靖 说 详细 点 ， 在 ej 单 
元 上 ， 


BizZ) 一 AZ) 一 (之 一 Z1-1)787 
PI)=NEN Cz)=(7— 7)h, 
所 以 在 单元 ej 上 


* 85 * 


V(r)=u p(T) uP (rz) 
=uN DT) + u Nn (x) 


Qs,=21/. CPV +ry: 2fr] 


-HP td) 


triusN Dn tu ND 2f (uN ED +t uNe dz 


0 NEY tronsny Ye 
十 


wa (PE i + rNinNn Jaz 


+1u 上 [P( 安 过) 二 CN#n): az 


i 


-| VD (rx)dr— 人 | NID (x) dz 
1 4 


令 
dN dN 
Ri. PE 一 二 rN EDNSD Jaz z 
《52) 
:=|. fN nCzr Ydr 
i 
Ca,P=1,2 f=1, ,7) 
则 
Ose t hf + hi b{us— bEu, 
我 们 要 写 出 导数 
2- 0Q, 
sD 
=ml 
来 。 因 为 


1 
1 一 六 1 二 人 


es 860 。 


+ hd din blu 
所 以 
Qin FAL 
+ bf bftiton 
可 以 看 出 其 余 Q, 中 不 含 4 ， 所 以 


00. so ao 


du, 一 oz Au. Ou, 


=hisniit Rin — bi + hitiu, + histu — bi 
现在 ， 只 要 令 呈 一 0， 就 得 到 关于 4 的 线性 方程 组 了 。 但 是 ， 


我 们 打算 在 每 个 单元 e, 上 求 出 单元 矩阵 ， 再 组 合 成 总 矩阵 。 记 
以 ， 对 于 Q, 而 言 ， 我 们 党 求 -3 一 及 3 如 下 


站 =hkiit_rt hi — bt 
(53) 
3 =hids Ri — bi 


未 为 夭 阵 形式 。 令 


和 人 | 
有 | [| 
于 是 
_6Q， 
Du ,-1 EE kt | He b! | 
0Q ， CR bi | 


ss SB7 » 


太 : 又 称 为 单元 e ,的 “单元 刚度 宗 阵 ”。 开 :和 B81 痢 携 带 着 单 
元 e, 贡 献 的 信息 。 我 们 按 单 元 逐个 算出 矩阵 乒 ' 和 妨 : 来， 就 完成 
了 第 二 个 步骤 。 

3 ) 将 刚才 得 到 的 单元 矩阵 玉 :和 号 :， 组 合成 总 矩阵 ， 写 出 
线性 方程 组 。 

这 一 步 又 在 许多 书 上 都 是 通过 例题 来 说 明 组 合 规则 。 但 只 死 
” 记 组 合 规则 ， 是 学 不 好 有 限 元 法 的 。 在 此 拟 说 明 有 关 组 合 规 则 的 
原理 。 

由 上 文 可 知 ， 我 们 要 写 的 线性 方程 组 就 是 


OQ _ 
Bo 09 (: 二 1 ,2， 得 
也 就 是 

GO Ou | | 
.OE 
du du | 
Qi ao | 
Ous Qus 
0Q， 
Ep Wi 


* 人 和。 


Ss ht + hau 6 
Se hu + hu 6 
Shyus th usb 


Ss =h97ty,., 十 R22H 2 b2? 


有 一 卫 


Ge- rit st hi 1 
二 人 


GO。_， =R21' ,a+ Rilu, 1— bi! 


n~1 
~- =hR?t,s + hisu,—b? 
hi 十 有 at 一 


将 它们 代入 前 面 的 方程 组 ， 得 


(Ris + ht du thisus+ hlsu — bi}—bi=0 


pa 


31 Us (R22 + RITI)U, 2 十 天 ?1 


— br?2— bls 0 (54) 
R2iluns t Chizt + hil)u, + Ri2u, 
— bi-t— b=0 


友 iy2 十 及 22 一 0 一 0 


由 表 为 矩阵 形式 


。 89 。 


?oy tzy 


Tl fT 有 
" ui TY 


TT- 


1! iy toty lty 
《99) 


0 ly ty 
3 Hy yt 


?3 ry t+?y 


二 90 
天 即 总 刚度 和 矩阵， 又 取 
b+ bf—hktsuo 
8 十 号 3 
bi+bt 
B= : (56) 
bss+ br-! 
bi I+? 
b3 
UT=Cu, Ws, He) (57} 
虽 方 程 组 表 为 
KU=B 
我 们 设 2&(zo) 一 如 为 已 知 ， 得 到 关于 未 知 量 4 ,9 的 线性 方程 
组 (57) 。 
如 果 u(26) 二 4 及 u(x,) 二 44 均 为 已 知 ， 那 么 未 知 时 训 只 有 


suas ea 了。 而 且 在 列 方 得 组 时， 当 一 0 这 一 方程 也 失去 
意义 了 。 于 是 ， 在 这 种 条 件 下 我 们 得 到 关于 未 知 量 4w:，u，… 
2 的 方程 组 

(Rss t+ RI uthIus = b+ bt— kisu | 

Riouit (hd, +t Ri)u, + Rl,ts =b?+ bob? 


多 


及 2128 8 十 《 玉 232 十 大 1T1)2s 2 十 下 ?21 人 -1 《58》 
一 消 2 2 十 器 1 

Riilw st (R22: + Ri), 
=6"1+ 6b — hr,p, | 


将 (58) 式 表 为 矩阵 形式 
K*U=B (59) 
在 此 


[| 


i pp i 


= 91 
b+ 63— kisuo | 
bi+b? 
Caro ‘8 .C60) 
bs 2+ bi-! | 
6 二 站 一 | 
又 总 刚度 矩阵 为 
| R32 + RE kt 0 0 
nl 0 
| hi: hs at 
\ hi kit 
(61) 


观察 《55) 及 (61) 式 ， 可 得 总 刚度 矩阵 及 和 天 * 的 “拼凑 法 ”， 
这 是 进行 有 限 元 法 实际 计算 时 必需 知道 的 方 法。 现在 对 和 矩阵 
太 * 的 “ 拼 资 法” 说明 如 下 〈 移 阵 攻 的 拼凑 法 与 天 * 的 摘 凑 很 相 
似 ， 仅 仅 在 矩 笑 的 布下 角 有 些 不 同 。 请 读者 自行 分 析 ， 基 好 结合 
等 例 ， 不 难 作出 有 关 结 论 ) : 

对 边 异 单元 Ee1， 算 出 82。。 对 单元 2;, 算 出 和 矩阵 


PD Pe 
-| 旬 名 | 咯 A 为 间 元 的 机 了 
21 2 


一 般 说 来 ， 对 单元 e ;来 说 ， 算 出 贡献 矩阵 六 二 
| 
| 二 8 | (j=2,3,"* ,7 ~1) 
《注意 :一 Ai) 对 边界 单元 e,， 算 出 R10。 

再 将 这 些 矩 阵 帮 ' 按 对 角 元 素 人 位置“ 焊接” 起 来 ， 就 得 到 和 矩 
隆信 *。 焊 接 方 法 如 下 表 所 示 ， 


ki | 


五 


中 若 拼 淡 矩 阵 兵 ， 则 需 算 出 


M1 R43s 


*。 92 


s(T9) 


一 一 一 一 一 一 一 -一 -一 一 一 一 一- 


-一 一 
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ya 


本 me 


rr rr 


‘oodnsor nsometor cag 


; H 
: ; 
H [3 总 
: 
; : 
了 ?总 了 
! 3 
: } 


+ 
?+ 


~ 
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nn 


tet 


DTT 


i n 


CTT 


nn 


ET TT 


9 3 。 
至 于 矩阵 8 的 组 合 ， 可 先 按 单元 ey 求 出 矩 陈 


pol 


2 


再 控 《〈60》 式 算出 忆 。 

4) 解 线性 方程 组 《57) 或 (59) 。 由 于 电子 计算 机 正在 普 
友 ， 如 何 解 线性 方程 组 在 这 里 就 不 必 多 说 了 。 只 说 明 总 刚度 矩阵 
除了 具有 对 称 性 、 正 定性 之 和 外、 还 具有 稀 羽 性 。 人 矩阵 的 每 一 行 至 
多 有 三 个 非 零 元 素 帆 ， 这 就 大 大 减少 了 计算 量 。 这 一 点 也 显示 了 
有 限 元 算法 的 优越 性 。 

用 有 限 元 法 解 问题 《39) 的 计算 步骤 已 经 韵 代 完毕 ， 现在 举 
一 例 ， 详 细 说 明 计 算 过 程 。 只 要 把 这 个 例题 的 计算 过 程 全 部 差 履 
了 ， 也 就 进入 有 限 元 法 的 大 门 了 。 

例 求解 

一 2 人 T) 二 MT) 一 一 了 XE 0, 1) 
US 


《注意 这 就 是 在 上 节 中 用 里 区 法 解 过 的 例题 7 

我 们 在 讲 有 限 元 法 时 ， 从 始 至 终 都 以 变 分 法 为 指导 思想 〈 实 
际 上 有 限 元 法 正 是 在 变 分 法 这 棵 古 树 上 开放 的 新 花 ) 。 为 了 使 初 
学 者 对 有 限 元 法 的 指导 思想 更 加 了 解 ， 我 们 在 解 这 个 例题 时 ， 先 
不 按 上 文中 总 结 出 的 四 个 步 又 来 计算 ， 而 是 更 直接 地 从 泛 函 写 出 
线性 方程 组 ， 再 求 其 解 。 

我 们 已 经 知道 ， 只 要 求 出 zeCz)E ci[C0,1], 使 泛 函 QQ[u(z)] 


| .ce (rz)+a(z)+2zu(z)]dz 取 最 小 信 ， 则 ws(zy 是 本 例题 
的 解 用 节点 0 4， 计 ， 守 。 工 将 区 间 [ 0 ,1 ] 划 分 为 4 个 单元 ， 
即 ， 


外 总 刚度 矩阵 的 每 一 行 的 零 元 分 布 ， 与 单元 的 篇 号 ， 排 列 有 关 。 这 一 点 ， 在 多 
维 问 题 中 独得 更 消 楚 。 人 得 总 阶 度 矩阵 总 具有 血防 此。 


® 日 4 。 


ur)= up.(z) 


:= 


但 zx= 0 = 已 经 给 定 ， 所 以 


WT)= > up.(7) 


| 


在 单元 e, 一 ( 0 ， 卫 ) 及 其 两 端点 ， 显 然 试探 函数 vi(z) 是 满足 
u(0)=0 及 4,( 二 )= 忆 的 线性 函数 。 由 图 15 可 知 ， 当 z EL0， 
1 

二 ] 时 ， 


x 
WAT) = = AUT 


4 


Hi4(T)= Au 


us 
: pa 
ws 全 | 
1 or 
~ 
| ” 
cg ~ 
0 于 一 一 一 
过 
图 15 图 15 


当 z E€ “= [村 ， 诗 ] 时 ， 由 图 16 可 见 


co 人 -re 


一 8iK2 一 477 十 22(47 一 1 中 
ur)= 4 — #1) 


—_{fl 3 | 
当 z Ee,= {=z 字 ] 时 ， 由 多 17 可 见 


uelz) = (3 — = ) 二 + si( -3)/ ] 
一 zz(3 一 47) 十 Ha(47 一 2) 


Ui(X)=d(us — 4,) 


当 工 és,=[(3, 1 ] 时 ， 由 图 18 可 见 
wz) 一 Ma(1 一 >) fo =u(1—7) 
u(r)= — 4us 


泛 函 Q[u(z)]= 人 Q， 


Q,=| ee(z)+ u(r)+2rulz ) dr 
4 


-|. (ut uit2ru)dr 
现在 按 单元 算出 Q,， 


! 
Q,= tieut TAZ) uit2T(4T) Ud 


QD 其 看 (45) 式 


sa 96 * 
1 
~(4 + 29 


1 
Q:= [1C((2— 4 + (47 —1)u,) +t 16{ wu, —u) 
+27C(2— A477) H+ (dr 1)u ddr 


1 
wm 十 HH) :+ 452 一 #1)* 十 (ut us) 


3 
Qs= jiK(ws(s-4z)+ 13( 4 一 人 7 十 16(8s 一 72 
+27((3—4T)u, + (47 —2)u) dx 


~ st Wa) + 4d{ 28 一 ts 2 + Cut us) 


ri 
OQ.= |3C16(1— 7)us+ 164u3 +2z.4(1 —z)us dr 
4 


由 此 看 出 ， 现在 泛 函 QCwakz)] 是 关于 变量 4 #2 ts 的 二 
次 函数 。 再 到 


-0 (+ 一 1,2,3) 


就 得 到 关于 4,，#,，，4 的 线性 方程 组 。 从 刚才 的 式 子 中 再 次 署 
出 ，Q， 是 2 的 疯 数 ，Q， 是 21， #1 的 遂 数 ， Qs 是 4， Ww 的 国 数 ， 
QQ, 是 4 的 函数 (在 上 文中 曾 指出 Q， 是 2 -ly ii :的 函数 ) ， 于 是 

oOQ O00 + 0Q: 
du 二 a =(s + 二 于 2 


+ 言 (由 +u) 8(u,— 4) > 


Sa 8 )us+ 4 


外 为 统一 起 见 ,在 计算 定 积分 时 ,一 律 用 梯形 公式 或 矩形 公式 求 近似 值 ， 不 管 能 
不 能 求 出 原画 数 。 现 在 用 的 是” /zy def (942) 6 -四 


|1o 构 


e 97 。， 


22 - Se+ Set) + (us 1) 
+ 汪 十 泪 -(ze 十 2s) 一 3(ws 一 uo) + 六 
=( 云 - 8 ) (去 - 8 jxs + 六 
1 


_ 1, .12 
=(3 + 16- 2 十 35 


再 人 29 = 0 (i==1,2,3) 得 线性 方程 组 如 下 ， 


16,25u1—7,873u, = 二 一 0,.125 | 


—7.875u,+16,25#,— 7.875us 二 一 0.25 
—7.875u, + 16.25us= 一 0.375 
再 表 为 罕 阵 形式 ， 即 


16.25 -7.875 0 | wu) [—o0.125 
I 16.25 -7.875 | za | 一 0.25 


0 一 7.875 16.25 | 2 一 0.375 

也 就 是 

K*U=B 
KK* 开 总 刚度 矩阵 ”， 正 如 上 文 记 措 出 的 ， 它 是 一 个 对 称 、 正 
定 且 稀 琵 的 矩阵 ! 

解 此 线性 方程 组 ， 得 
r ud=—0.0354 
ud 一 0.0571 
2 一 一 0.0439 


于 是 ， 我 们 用 有 限 元 法 得 到 本 例 中 微分 方程 边 值 问题 的 近似 解 
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外 
W(x)= uip (x) 


nw} 


详细 写 出 u(x) 的 表达 式 ， 就 是 

f —0.1416Z7 

| 当 z e[0.,] 

| 一 0.0354(2 一 4z) 一 0.0571(4z 一 1) 
| se 


一 0.0571(3 一 4r7 一 0.0439(4 工 一 2 


sc( 3 


一 0.17586(1 一 了 ) 
se 


人 
re 


1)= —0.0354 


A 
Cy 
加 


和 
F 
OO 

中 
1 

心 

OO 

Cn 

| 

jp. 


i ofa Ni oj 中 


~ 
An 


人 

[3 

Am 

wl | 
> 
| 
© 
心 
心 
ow 
二 


)= 一 0.0220 


再 将 此 有 了 痕 元 解 与 里 兹 法 解 及 精确 解 列 表 比 较 〈 表 2 ) ， 以 便 观 
察 误差 分 布 状况 (参看 本 章 § 1 中 的 天 1) 


a 099 。 


: 血 兹 续 解 近似 值 | 有限 元 解 近似 什 


如 果 节 点数 昌 增多 ， 像 刚才 这 样 先 写 出 泛 函 人， 再 对 4, 求 导 
就 太 麻 类 了 中 。 所 以 我 们 还 要 学 会 先 被 单元 求 出 单元 矩阵 ， 再 拼 
成 总 刚度 扎 阵 人 《或 天 *) 及 B 的 方法 。 为 使 读者 能 掌握 这 种 方 
法 《就 是 现在 通用 的 有 限 元 算法 ) ， 我 们 对 于 刚才 这 个 例题 ， 用 
上 文 指出 的 四 个 步骤 再 计算 一 次 。 


第 一 步 ， 取 节点 0， 坟 地，-，1， 将 区 间 [0，1] 划 


oa (Rn 
第 二 步 ， 按 这 4 个 单元 逐个 计算 单元 和 矩阵。 为 阅读 方便 ， 把 
按 单元 计算 时 要 用 的 公式 汇集 如 下 : 


dN dN£n 
po=| [as 4 


Te + rN Ng laz 


5 一 | FNecDdz 
上 


中 在 实际 问题 中 ， 节 点 往往 有 成 千 上 万 的 情况 。 


。 100。 
NED (Cz)=(7;~ 7)/h; 


(hj=r;— 72;-1) 


NID(r)= (rr h, 
在 单元 e; 上 ， 要 算出 


Ric| ki ki | sz-| | 
\ ki) R2, | 


因为 (0)==0=4(1)， 所 以 我 们 要 用 (60) 和 〈61)》 拼 出 到 * 和 
B。 


现在 z = 了 /4， hi= 于 《了 二 1,2,3) 

所 以 

No 一 (二 -z)/ 全 -4 CN Cz) = 4 

te (i 

CNED (rT))’=4 
x 

P(r)= 1 rr)y, f(r)=—z 

所 以 


4 
= CNED HCN ED )9dz 
=| ce+(-4z)59dz 
二 
~ 二 [16 + (7 —4r)’) 过 
ed -di 16 
(3 而 ) 


了 
= 人 CNN CN ED YH NED NAD Ydz 
wr" 


"101。 


了 
一 | 5-16+(7 一 4z)(4z 一 t1D]dz 


了 
2 ,C16+(47—j+ 1)*Jdz 必 4 地 


又 
4 
bi=|" (一 z)07-47)dz 
1 1—27 
2 Cr 
上 
7 i 1 一 27 
2 rz)(d4rT—j+1l)dr 二 人 
4 
所 以 
4 -名 | | a-2n /6 
Ki= : Bi=: 
63 1 
一 语 个 | (0-21) /04 
由 《60) 式 可 知 
1 1 _63 
46 + 416 16 0 
站 一 _63 1 有 _63 
人 16 416™ 16 16 
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63 i 

We | 
由 《60) 式 得 

" 

bi+b? | | 

,64 

| 8 

=| 5b2+4+53 1=| ~ 

B 2 十 上 i | BA 

12 

# de 

看 了 十 所 | 6 


} 


1 _63 4 1 
Sa ,ee S 
63 1 _ 6， 2 
16 “8 一 16 | 
63 1 3 
之 二 = | 一 -一 
16 88 ll™] | 


化 方程 组 显然 与 《61) 相同 。 

第 四 步 ， 解 这 个 方程 组 。 

这 些 就 是 用 有 限 元 法 解 1 维 问题 的 实际 步 又。 初学 者 如 把 这 
一 节 的 内 容 掌 握 了 ， 就 不 但 学 会 了 关于 1 维 问题 有 限 元 法 的 实际 
计算 (特别 是 关于 总 刚度 矩阵 的 拼凑 法 》 ， 而 且 懂 得 了 为 什么 要 
这 桩 算 。 否 则 ， 如 果 只 知道 机 械 地 按 规则 计算 ， 遇 到 复杂 的 实际 
问题 就 难以 处 理 了 。 

本 书 是 在 变 分 法 的 基础 上 来 介绍 有 限 元 法 的 。 试 图 使 读者 在 
从 得 变 分 法 以 后 ， 进 一 步 领会 有 限 元 法 的 数学 思想 。 在 这 种 数学 . 
思想 指导 下 ， 透 彻 地 掌握 1 维 问 题 的 计算 步骤 。 关 于 多 维 问 古 的 
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具体 计算 及 实际 应 用 ， 在 这 里 就 不 谈 了 。 读 者 在 看 完了 本 章 之 
后 ， 应 当 参 阅 其 它 关 于 有 限 元 法 的 书 ， 使 白 已 能 够 在 理解 的 基础 
上 掌握 多 维 (特别 是 2 维 ) 问题 的 计算 方法 。 

最 后 ， 还 须 指出 ， 现 在 已 经 有 关于 有 限 元 计算 的 多 种 通用 程 
序 。 读 者 在 熟悉 有 限 元 计算 以 后 ， 再 学 会 使 用 这 些 程序 ， 就 能 在 
电子 计算 机 上 迅速 而 准确 地 解决 各 种 实际 问题 了。 


Cl 


~ 


[2) 
C3] 
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第 二 部 分 外 推 法 
第 五 章 外 推 法 


外 扒 法 是 一 种 简便 易 行 而 精度 又 很 高 的 数值 解法 ， 在 国内 党 
握 了 这 种 方法 的 人 还 不 多 。 但 这 种 方法 实用 价值 很 高 ， 应 当 介绍 
给 我 国 工程 技术 界 。 在 本 章 中 对 这 种 算法 及 其 应 用 作 一 较 系 统 的 
介绍 。 


8 1 外 推 法 的 实例 


我 们 打算 先 谈 两 个 最 简单 的 问题 ， 从 这 两 个 例子 中 引出 外 推 
法 的 指导 思想 。 

例 1 求 + 的 近似 值 
我 们 知道 直径 为 1 的 贺 的 内 接 正 # 边 形 的 周 长 为 


d,= nsin 
我 们 知道 
d, =xsin® fF—>7 (n>00) 
现在 求 出 qd.， 作 为 x 的 近似 值 。 
先 取 1 二 12， 得 du 一 12sin15 一 3.1058286 再 取 + 一 24， 
得 d2,==24sin27 一 3.1326285 这 两 个 值 误差 都 比较 大 。 可 是 ， 如 


果 我 们 用 下 面 这 个 算式 
d* = 4d da) 


» 105« 
对 上 述 d1z 及 qd,, 进 行 “ 外 推 >》， 就 得 到 
de 一 (4x 3.1326285 —3.1058286) 


一 3.1415618=3.1416 
是 不 是 碰巧 这 样 呢 ? 再 试 一 次 。 这 回 ， 我 们 先 求 出 
。 。 轩 
de 一 6SmT = 3 a di 一 12sinjz 一 3.1058286 


再 用 刚才 的 办 法 ， 得 
de 一 于 (4d 一 ae) 一 3.1411048 


也 得 出 比 d。 和 qd. 精确 得 多 的 x 的 近似 信 。 
一 般 说 米 ， 我 们 先 求 出 x 的 近似 值 4. 和 d,,。， 再 用 如 下 外 推 公 
de = (4ds, —d,) (62) 


就 得 到 比 d, 和 d;, 精 确 得 多 的 近似 值 d&*。 这 种 由 较 煌 焰 的 近似 全 
“外 推 ”出 较 精 确 近 似 值 的 方法 ， 洋 在 太 简 便 了 。 只 要 对 原来 的 
两 个 近似 值 d, 和 dd 进行 三 次 算术 运算 就 行 了 ! 这 样 好 的 算 法 ， 
我 们 在 实际 工作 中 为 什么 不 用 呢 ? 

现在 来 论证 外 推 公式 〈62) 。 由 熟知 的 Taytlor 公 式 得 


COS OX 
a (0< 9<1) 
31 51 
所 以 
Tt 
x x 1 a CoS 1 三 mT 
Snirn mt 0 "nm 
eR 
Tn 6 nt vn 
在 此 
TT 
cos 87 
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名 
C EL 
于 是 
df 一 一 一 一 (E> 1 
nS KH nn 6 于 ngs 
国 nm 1 | 
EE 十 Can 
Tt? 
令 百 =C， 则 
全: 
可 《63) 


在 《〈63) 式 中 ， 常 数 C 与 4 无关 。 常 数 C, 昌 然 与 4 有 关 ， 但 


并 5 


C. 有 界 ， 即 存在 与 ? 无关 二 
， 即 存在 与 无 关 的 正 数 M= 茸 ， 合 
IC,l <M 
所 以 
Ci 1 人 
aa。 和 《27)3 + Coane = 一 An: 二 了 (64》 
于 是 
Cc 
dds,—d,=37+ pp a 
1 1f€ 
了 d _ 各 ~ “2n 1 
3 (4 2 d,) 郊 十 3( a a 2 
令 
1 Cs 
3 )=e. 
出 


1f Cd _17M 
[as IC ) (+W) -5McU 


ll 1 
ad*= 3 (4 一 gd ) 一 区 十 Gy (65) 


。107。 
由 《〈63》 及 (64) 可 知 
df 一 了 十 o( 六 


aa ,=x+O( 立 ) 
这 就 是 说 ， 误 差 Id, 一 zi 和 jd 一 | 均 不 低 于 坪 阶 《在 许多 文献 


中 说 成 “误差 为 汪 3 阶 ”》 
《65)》 式 可 知 


这 就 是 说 ， d* 一 计 (4d。. 一 qd,) 具 有 不 低 于 去 阶 的 误差 ， 可见 d* 
比 d .和 d,, 的 误差 小 得 多 ! 
例 2 求 定 积分 的 近似 值 
对 于 定 积分 
7 =| rzyaz 
| 13 分 为 N 个 


等 长 的 小 区 间 ， 分 点 是 zo 一 0 ，xz1 一 六， Ts , 


Ns od TAN 一 六 
=1, 记 


fro)= Fo, fOr)=f, 2, i Ey 
取 | 


Tvw=- 苑 A ftet n+ -这 |) 


~{ f(r)dr=1 
这 种 公式 是 大 家 郭 知 的 。 境 在 ， 我 们 要 把 用 梯形 公式 得 钊 的 近似 


! B., | 1 
人 D 着, 注 尼 工 -Hf (1 与 "无 关 》， 则 记 8,=0(ar) 
对 
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值 1n 及 1,ws 用 外 扒 法 处 理 ， 得 到 党 好 的 近似 值 。 
令 方 = ， 又 记 1 wy 为 1;,。 用 著名 的 尤 拉 - 表 克 游 林 《Euler- 


Maclaurin)D 公 式 可 证 明 如 下 结果 ， 即 ， 闭 被 积 函 数 y 二 f(z) 在 
[0，13 上 的 2 +2 阶 导数 部 *” 连续 ， 则 
IT,=C HC htC Rt CR Cnh N+ Cyrhin’s 
66) 
1 
oN Se WALLE: 
o 


CoyCioCos ;Cn 都 是 与 天 无 关 的 常数 ， 系数 Cir 有 界 《 即 存 
在 与 六 无 关 的 正 数 邮 ， 使 ICwsil 委 好 对 任何 和 成立) 吧 。 所 
以 
= 一 了 CC 天 十 加 (六 4 (67) 

《66) 及 (67》 式 指出 1 与 定 积分 精确 值 的 关系 。 它 们 都 称 
为 六 的 渐 近 展开 式 。 对 外 推 法 来 说 ， 关 键 问题 就 是 求 出 这 种 渐 近 
展开 式 。 有 了 这 种 展开 式 ， 就 可 以 用 外 推 法 求 出 更 精确 的 结果 。 
做 法 如 下 ， 

由 《67》 式 可 知 


l= +o(2) ro 


二 了 十 hi + OCRID 


所 以 
4 有 一 4 二 Cd OCR 
和 


一 47 十 CR 二 OCR] 
B41s -1,)=1+00h) (将 〈67) 式 代入 此 式 堪 端 ) 
这 就 是 说 ， 如 果 用 7*= 寺 Us 一 1) 作为 1 一 | f(x)az 


加 关于 Euler-Maclaturin 公 式 及 其 证 明 , 可 参看 球员 生 及 黄 友 详 所 落 丸 数值 
下 近 》p.157~160 
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的 近似 值 ， 那 么 ， 它 的 误差 是 不 低 于 如 阶 的 。 原 来 用 梯形 公式 求 
出 的 近似 值 1 和 7 只 有 不 低 于 外 阶 的 误差 《 见 (67) 式 》。 

我 们 己 经 找 出 了 比 梯形 公式 更 好 的 求 定 积分 近似 值 的 公式 。 
现在 把 求 1* 的 算法 直接 号 出 来 。 因 为 


nor 人 -+ 区 
+ AO A den) 
-60 +27 ( 忆 ) + 21( 各 ) + 2/( 襄 ) 
i gt pe ] 


= (0+ f(r) + 7 
7 


4 三 3 270)+ 4f (a ) 二 4/( 方 ) + af ( 二) 
ep eg) 
于 是 


1*= 诗 1414-1,) 


-an + (zr) + aa 


@ 吧 然 0(- 生 )= O 〇 (kb 


*。110。 


+ 局 本 2/ (3) 小 十 4 2 让 1)+ Fa] 


《68》 

这 论 是 熟知 的 卒 下 生 (Simpson》 公 式 ! 现在 我 们 知道 了 用 辛 下 

后 公 式 求 定 积分 的 近似 值 时 ， 误 差 不 低 于 如 阶 ， 即 - 太 + 阶 。 可 见 
辛 卜 生 公式 比 梯形 公式 精确 得 多 ! 

我 们 从 梯形 公式 出 发 ， 外 推 一 次 就 得 到 精确 度 高 得 多 的 广 下 

生 公式 。 不 过 外 扒 方 法 也 是 多 种 “ 样 的 。 例 如 ， 由 《66》 式 可 知 


1,=I+Cht+ Ch+ OCh') 


A 1 2 .了 一 二 a 
全 =7 + Ch + 16Cah4+ OCh?) 


a 1 i 4 本 
全 一 了 十 18 CC 下 十 a56 Ce + O(hs) 
于 是 
afl,+ bil, + CTs 
2 4 
六 人 
={ai+b oT+(a 一 二 全) CR 
要 加 b C a 
+(a [Cn 二 (六 
令 
[a | 之 + 之 = 作 
| 4 ”16 
1 
得 
| Q en 
64 


人 
2 7 64 
所 以 
1 4 6864 ， a 4 
py ol 7 人 由 == 。 6 ” 
圳 7.- 二 1 时 1=7+Oi/ 多 (到 C 2) 


可 网 ， 若 用 
二 50 一 2073 + 647 41 
《69》 
来 求 定 积 分 了 的 近似 值 ， 则 
fi** OCh') 
误差 将 不 低 于 到 阶 。 如 果 把 


六 一 = (+ ft fr 十 fw-i 十 总/ ) 


代入 (69) 式 ， 可 得 出 比 辛 卡 生 公式 更 精确 的 波 尔 〈Boole) 公 
式 。 但 是 这 公式 太 繁 琐 、 用 起 来 不 方便 ， 在 这 里 就 不 写 册 来 了 。 

从 以 上 叙述 中 可 以 看 出 ， 将 近似 算法 进行 “外 挫 处 理 ”， 御 
往 得 出 更 精确 的 公式 ， 这 也 算是 一 种 “ 精 加 工 ” 吧 ! 在 数值 计算 
的 许多 方面 ， 用 外 推 法 常常 使 我 们 找到 更 好 的 算法 。 外 推 法 是 产 
生 新 公式 的 一 个 涯 泉 。 


§2 外 推 法 在 常 微分 方程 数值 解 中 的 应 用 


在 本 节 中 ， 我 们 将 说 明 怎 样 把 用 差分 法 得 到 的 常 微分 方程 数 
值 解 进行 “外 推 处 理 ”， 从 而 大 大 提高 精确 度 。 

1， 初 值 问 题 的 外 推 

对 于 如 下 常 微分 方程 初 值 问题 


i 


y(0)= yo ‘9 
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有 如 焉 结论， 如果 方程 右 端 函数 六 >,4)》 在 区 域 
DD, a 人 lb, 一 co<2 所 十 cp 


中 连续 ， 且 储 导 数 -9 有 异 ， 妈 有 名 之 0， 使 


2 x, y) |<MY 


在 区 域 D 上 成 立 ， ee a 在 x ECa,0] 上 有 唯一 的 连续 可 
徽 解 由 。 
现在 ， 设 f(z,y》 在 区 域 
门 ，0 二 z 近 11， 一 co 所 y 扩 十 co 


由 连续 且 信 导 数 -9 三 有 界 。 求 《4) 在 =E [0,13 的 数值 解 。 


(1 ) 首先 列 出 与 《4) 相应 的 差分 方程 。 将 [D，1] 分 成 
个 相等 长 的 小 区 间 ， 它 们 的 长 h 二 站。 分 点 为 


To 二 Uy Tz 二 有, Ts=2h, a Tz, = 二 = 1 
因为 


所 以 方程 
y(t)=f(r,y(r)) (Cy(0)= #0) 
， 权 用 方程 
E+) yr) 
CA )! 
yt0)= vo 
近似 好 代替 。 
为 了 方便 ， 将 点 z 一 0,8,26…, (2 一 1)7 的 集合 记 为 
Qi;,={r|r=0,h,2h," ,1— —1)h} 
但 《A4’》 的 解 应 当 与 原 方程 的 解 有 区 别 ， 将 它 记 为 vz) 比较 合 


=f/(r,y(r)) 
《当空 一 0 大 ;2 《一 1) 丰 时 》 


DW 参 略 ].Steer 及 R .Bulirsch 所 著 latroduction to Numerical 
各 aaiysis》 第 7 总 第 1 节 
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和 送 。 于 是 得 到 差分 方程 
SD rEn, 
yi(0)= yo 
《2 ) 解 差分 方程 (41) ， 求 出 y,(z5) 《rzER,) ， 把 它 作 
为 问题 《4) 的 近似 解 。 
方程 《4;) 很 容易 解 出 。 记 
yA KhA)=YE (K=0,1,2.,n) 


CA.) 


即 

Y,(0)= Yo Yh)=Y, ,0 Ynh)=¥ 
由 《4)》 可 知 
Y=yo 
Y=Y,+hf(0,Y,) 
Y=Y,+hf(h,yY,) 
了。 —Y,+ hf(2h,Y,) 


Very + hi(Kh,Y xr) 


ee ee :+ hinom 1)h,Y,.) 

这 是 一 种 递 推算 法 ， 一 下 了 于 就 把 站 1, 了 ,,…,Y 了 ,全 解 出 来 了 。 出 
此 看 出 方程 《A1) 的 解 存在 县 唯一 。 

在 理论 上 ， 可 以 证 明 ， 
若 帮 rz，29) 在 区 域 局 (ac 委 z 委 8， 一 co<g 所 +co) 中 ， 有 六 士 2 阶 
连续 有 界 偏 导 数 ， 则 方程 

yr+h)= yr hf(r, y(t)) 

| 


ZE 人 一 人 Iz=atih,i=0,1," nn)} 
网 Vi(a)=yo | (5 =) 
的 解 yVAz)? 有 如 下 的 渐 近 展开 式 


yr)=y r+ A (T+ hi (Cr) + 
+ A (Tt RMBEys Crh) (70) 


在 此 式 中 ，2%(z) 是 方程 


"114* 


| rE(a,b) 
y(2)=ye 
的 精确 解 。VCr), 玉 (z) ,JVw{z) 是 与 4 无 关 的 汞 数 。 它 们 
均 在 Ca，613 上 有 定义 且 连 续 。 又 
Wi(a)=0 K=1,2,.°Ny 
又 余 项 ,u(r ;8) 满 足 
To 


Easti(zii1 < 生 M(z) 此 式 对 任何 产 = 扩 一 二 二 
n =1,2，…: 成 立 ) 
(70) 式 给 外 推 法 提供 了 理论 基础 。 

(3) 现在 可 以 进行 外 淮 了 。 由 (70》 式 可 知 ， 若 /(xx,y) 
在 区 域 万 

D: 0<rE1, COLL+ 
中 有 4 院 连 续 有 界 偏 导数 ( 取 广 = 二 2》， 则 其 分 方程 (41) 的 
解 y,(x} 有 渐 近 展开 式 

yA YT thW (r+ (7x) + hs Btr:h) 


=y(r) + AW (r+ O(h)® (71) 
于 是 
yi(T)= yr) 十 Ey, (7)+ Och?:) 
所 以 
2ys(T) YAT) ET OCF’) 
肥 


yrT)=2Y (TI) YT) {72) 


全 人 参半]. Sitiosr 及 RR。Rulirsch 所 营 《lntroduction to Numerical 
Analvsis 为 第 ?7 音 第 2 节 
加 开 :(zr) 在 [ap] 尝 续 有 界 。 又 瑟 。 (zs) 在 LasJ 有 办 于 是 
Rl | | 了 + | 忆 村 (常数 ) 


ve hi rT) + hE Th) = OU!) 
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则 
yr{r)—=y(r) + O(h’) 
出 《71》 式 看 出 
yA(r)= yr + O(h) 
所 以 ， 和 由 外 推 得 到 的 近似 解 z*z 比 (4,) 的 解 要 精确 多 了 1 
举 一 个 例子 。 对 于 初 值 问题 


1 (Cz) 一 以 并 ) XE(O0,1) 


(A) We 〈 邑 六 zy) 三 2 


求 它 的 数值 解 。 
第 一 步 ， 别 出 差分 方程 
Yar+h) yr Hhy (rr) (r=0,h,2h,.,(n=1)h) 
yi(0) =1 (4 = 元) 
再 解 此 方程 (4 )。 令 
Yr=y, KA),R=0,1,,f 


CA) { 


将 (有 41) 表 为 
人 Yr (K=0,1,2,%,(n—1)) 
Yo,=1 
于 是 
Yo=1 
Yi=(1+h)Y,=(1+#) 
Yo=(1tH)Y = (+h)’ 
设 i 
了 一 (1 十 上 )” 
则 
Yeprn=t+th Y=(+h)°*! 
所 以 
Yxr=(1+A)* (K =0,1,°" ,7) 《73) 
这 就 解 出 了 方程 (4;)。 


现在 进行 数值 计算 。 取 
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8 
由 公式 (73)， 得 差分 方程 (41) 的 解 也 ec 一 (1+ 二 


了 ,=1.125 

Y,=1.265625 

Y ,=1.4238281 

Y ,=1.6018066 了 
| Y ,=1.8020325 

Y=2.0272865 

了 ;一 2.2806973 


Y,=2,5657845 


我 们 已 知 (.4) 的 精确 解 是 
y(r)=e” 


分 出 职 二 0 证 和 
y(z) 二。 在 这 些 点 处 的 值 如 下 
-i= 1 
| y1=1.1391485 
ys=1.2840254 
ys=1.4549914 
y4=1.6487213 
yo=1,868246 
ys=2,1169 
7 一 2.3988753 
ys—2.7182818 
对 照 这 些 数值 ， 可 见 差 分 解 yr 与 精确 解 yr 之 差 四 。 
为 了 外 推 ， 再 求 出 ys(z)。 由 《73) 式 可 知 


全 ”在 计算 yYk 时 ， 还 有 由 了 四 含 五 人 阁 成 的 误 蒂 。 在 本 书 中 ,就 不 讨论 这 种 “ 舍 人 
误 盖 ? 5 
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ed) 0 


于 是 


由 此 得 


Y 

Yt=1 .1289062 
Yt=1 .2744293 
Vs= 1.4387111 
=1.624170 
Yo = 1.8335356 
Vos 一 2.0698898 
一 2.3367114 
一 2.637928 


第 二 步 ， 将 上 面 求 得 的 两 组 差分 解 进行 外 推 处 理 。 我 们 用 公 
式 〈72) 进行 外 推 。 把 〈72)》 式 写 成 
Y=2Y gg Yr 
在 这 里 
Yr,s=y,( KA), Y «r=Ys(Kh) 


由 此 得 
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‘Yi=2x1—-1=1 

| Vt*=2x1.1289062—1.125=1.1328124 
,Vi=2x1.2744293—1.265625=—1.2832336 
| y $=2x1.4387111— 1.4238281=].4535941 
/Ji=2x1.624170 一 1.6018066 一 1.6465334 
Ys=2x1.8335356 — 1.8020325=1.8650387 
y=2x2.0698898 2.0272865 一 2.1124931 
/# 一 2x2.3367114 一 2.2806973 一 2.3927255 


Y=2x2,.637928 —2.5657845==2,.7100715 
现在 将 上 面 求 出 的 差分 解 Yx, 及 Yx,s， 外 推 解 Y* 和 精确 解 
#z) 一 er 的 值 列 成 对 照 表 ， 以 便 读者 观察 〈 取 小 数 4 位 ) 
表 8 


L0000 
ls.1250 
1,28656 


1.4238 
1.6018 6465 
1.8020 : 1.8650 
2.0273 | 2.1125 
2.2807 | 2.3927 
2.5658 | 2.6379 2.7101 2.7183 


在 此 衣 . 上 上， 我 们 又 一 次 看 到 外 推 法 的 入 大 优越 性 。 我 们 把 误差 较 
大 的 差分 解 了 x.。 和 了 x,s 用 外 推 公式 《74) 处 理 一 下 ， 只 纶 费 了 
很 少 的 劳动 ， 却 得 到 了 精确 得 多 的 结果 ! 

2. 边 值 问题 的 外 推 

对 于 如 下 的 常 微分 方程 边 值 问题 


—y (rT ay(r)=/(r) rE(0,1) 
yO) 二 yo，y(1)= 二 y， 4 为 常数 ( 设 4 守 0 


全 当 4 半 0，f (x) 连 总 ， 可 以 证 朋 (B8) 的 解 友 在 ， 唯 一 。 


cB) { 
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《这 是 线性 常 微分 方程 的 两 点 边 什 问题 ) 了 D 
我 们 打算 用 外 推 法 求 它 的 数值 解 。 仍 取 二 个 步 又 
《1) 列 出 与 《已 ) 相应 的 差分 方程 。 
将 区 间 〔0，1 3 ， 等 分 ， 令 有 = 十， 分 点 为 
To 一 0z 一 Ra 一 284oo37T ==1 
用 二 路 差 商 
[yr ) 2 (FE) Fyre? 
来 近似 地 代替 二 阶 导 数 y”(x;)， 得 到 差分 方程 
[ — [yriri) — 28. CT) TT YR) RS 
(BY) ( =—ay(r)+f(z,;) 
(yO yo YA)=Y, CF =1,2,° ,1) 
《2 》 解 差 分 方程 (8,) ， 求 出 它 的 解 y,(z;)， 把 它 作为 
问题 (8》 的 近似 解 。 
方程 组 (8,〉 是 线性 的 ， 它 的 解 存 在 且 唯 一 ， 这 一 点 只 要 
观察 它 的 系数 矩阵 就 可 以 肯定 ， 不 详细 说 了 。 
在 理论 上 可 以 证 明基 (8B) 方程 右 端 函数 /(x) 在 [0，1] 
上 有 4 阶 连 续 导 数 / 3 ， 则 《8,〉 的 解 y,(x) 有 如 下 的 浙 近 展开 
式 
yr)=yr) + (72) + hE(r,h) (75) 
在 此 式 中 ，WWV(z} 是 与 无 关 的 前 数 ， 它 有 4 阶 连 续 导数 。 函 数 
过 (人 z) 有 界 。 我 们 在 本 节 的 最 后 来 证 明 〈75)》 式 。 初 学 者 和 只 关 
心 外 推 法 的 应 用 的 人 ， 可 以 不 看 这 个 证 明 。 
《3 ) 外 推 。 有 了 《75 这 个 渐 近 展开 式 ， 就 可 以 外 推 了 。 
和 过 去 的 做 法 一 样 ， 因 为 


J RV AY A 
y= ys) + (3) Wr) + (2) F(z 2) 
QD 方程 中 y*(7) 的 系数 取 为 1， 是 现代 数学 工作 者 为 讨论 一 些 理论 问题 {如 


4 广义 解 *” 问题 的 方便 市 设 的 。 如 果 您 不 习 刁 这 一 点 ， 可 将 此 方程 的 请 一 
项 均 药 以 ~- 1。 
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yr) + Ez) + OC) 
又 《75》 式 可 写成 
yr yr + hr) + OC 
若 取 
yr) = CaysC) yr)) (76) 


则 
yrT)=y(r)+ OF’) 


还 是 举 一 个 例子 。 对 于 边 值 问题 
人 rE(0,1) 
yt0)=0, y(1)=1 


用 外 推 法 求 数值 解 。 
第 一 步 ， 写 出 与 (8) 相应 的 差分 方程 。 


将 C0，17 区 间 “ 等 分 ， 令 有 二 二 ， 分 点 是 
To 一 0 T= hth T=nh=1 


差分 方程 取 为 
CO pe 


yi(0)=0, yt1)=1 《fj =1,2,.."…,n—1) 
也 就 是 
UACIRDDACETAD LACIDER /ACIND 一 站 
(B,) { 
yi(0)=0, y,(1)=1 ,名 | =1,2,'" 7 一 1) 
第 二 步 ， 解 方程 组 (8B,) 


取 n 二 2， 则 户 二 地。 这 时 只 有 三 个 分 点 ， 即 


(CB) 


Z0= 一 0， zi 一 上 一 未 ， rs=2h=1 
已 知 
yslT0) = 0, YilT2)=Y(1)=1 
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只 有 yz1)=y,( 吝 ) 的 值 未 知 。 
在 (B,) 中 ,到 j 一 1， 有 hh 一， 得 
yzo) — (2+)y (rt yr)= 0 


所 以 


9 
-Vr) +t1=0 


VC) = 地 
于 是 


y(t1)=1 
为 了 淮 备 外 推 ， 再 求 y§(Ti)e 这 次 ， 取 


n 二 4, h' = 地 


将 区 间 [0，17 4 等 分 ， 分 点 有 5 个 ， 即 


记 
1 一 全 = 寺 ， yi (Ti) = Ys 
则 方程 组 《8,》 化 为 


yi-i— (2+h y+ yr =0 
Vo= 0, Ws= 1 


详细 地 写 出 ， 就 是 
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| 1621 一 2 = 0 


.33 
-yt16 Ys=0 
| i 
i Ya 18z= 一 
于 是 
1 
| y=0.21511777=y; (J) 
ys =0.443680402 = ys(2) 
\ ys=0.6 eR 
ys=0.6999662 ys( 
三 步 1 
第 三 步 ， 针 推 。 用 (76) 式 ， 在 7? 二 处 外 推 ， 得 


(2)= 3 4 二 -ww 


4 
3 G x 0.443680402 -二 一 0.4434257 


与 精确 解 对 比 一 下 。 问 题 《<B》 的 精确 解 是 
y(x)= 2 2 一 6 一 ) 


在 z = 过 处 ，%(z) 的 值 是 


v3)=a( ex = e-7) =0.443409442 
外 推 的 效果 是 明显 的 。 


由 《76》 式 可 知 ， 用 这 种 方法 只 能 在 点 + 一 寺 处 作 外 推 (在 
其 余 点 处 因数 据 不 足 ， 不 好 外 推 》。 这 就 使 数值 解 在 点 > 一 士 及 
z 一 -4 处 误差 较 大 ， 造 成 数值 解 在 各 点 处 误差 不 同 阶 , 最 近 ,我 国 
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学 者 陈 传 斑 和 林 群 建立 了 一 个 很 有 用 的 新 公式 @， 用 这 个 公式 可 
以 在 其 余 各 点 处 同时 进行 外 推 。 现 在 就 来 介绍 这 个 公式 并 且 用 它 
在 分 点 z = 访 ， 闻 处 进行 外 推 . 
用 分 点 


wo=0,r1=h ,=n 一 1 


把 区 间 50，173 1 等 分 。 再 出 分 点 


i hk 月 
0 二 0 ， Ti 二 2: 2 Ls=3 


z= 2 人 1 
把 C0，12 2n 等 分 。 把 "等 分 的 分 点 的 集合 记 为 
9. = Ax; jz ,= = ,f=0,1, ,ny 
把 22 等 分 的 分 点 的 集合 沁 为 


如 果 
TIE Ld TirniE ND, 

但 是 z; 和 zj 的 中 点 
zj=(x)+ En) /2€ AN 


这 时 ， 我 们 从 差分 方程 (8,) 中 可 解 出 
YTi)s Yrsr1) 
和 
ys(x), Yari), yal Tir) 
这 几 个 值 。 于 是 ， 用 《7?6) 式 外 纵 ， 得 
1 
y*(z) = PACERETACS)Y 
@ ”用 此 公式 还 可 以 对 有 耻 元 解 作 外 淮 。 此 公式 已 又 引起 计算 数学 界 的 重视 。 
加 ”Qi 入 自家 示 集 合身 ,74 与 集合 但 ;多 其 集 ， 带 点 二 局 于 sj 但 不 妃 于 


Mis 
忆 表 示 对 两 集 台 施行 的 减法 运算 。 
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vy*(z1) = 4yi(z1+1) — Yi(Tit1) ] 


但 是 y,( 瑟 ?的 值 不 知道 ， 所 以 不 能 用 (76) 式 在 (z; 处 外 推 。 陈 
传 狂 一 林 群 公式 是 
yr CZ = YI + (Ya Yt) + yslCrin) 


—ylTiri )| (77) 


现在 用 《77) 式 在 分 点 + 一 十 及 zx 一 写 处 进行 外 推 。 得 


v(t)= 由 (全 )+ (460) -wo0) 

+ v3) -3) 
=0.21511777 + 去 [0. 443680402 专 ] 
=0.2149903 

(3)= A + 并 vs(3 )=y ys( 雯 3)+ ys(1) 


-v01) |= 二 0.6999662 十 [°°: 443680402 一 了 ] 


一 0.6998388 


(e* 一 e-*) 的 值 是 


在 此 二 点 处 ， 精 确 解 4(z) 一 二 
‘a 
"(4)- 


将 以 上 数值 列表 以 便 对 比 


下 
( ex 一 e 1)=0,2149523 


立 
e “4) 一 0.6997241 


Ts 


| 4 2 
本 
Wa tT ER 
| EE EE 9 
Us (ray 0.21511777 0 .4436804{2 0.6999662 
Yr,) D21499038 0 .4434257 0.6998388 
YT 0.2149523 0.443409442 0.6997241 


从 这 个 例子 可 以 看 清楚 如 何在 常 微分 方程 中 使 用 外 推 法 。 
现在 对 (75〉》 式 给 出 证 明 。 
对 于 问题 


—y (zr)+ay(r)= f(r) rEC0,1) 

y0)=y0, YL)=Y 4 为 常数 。 现 在 设 4 之 0 
这 是 线性 常 微 分 方程 边 值 问题 。 在 微分 方程 理论 中 已 经 证 明 、 落 
f(z) 在 [C0 ,1] 连续 , 则 (2) 在 fo,1] 中 有 了 唯一 和 解 y(z)。 这 时 因 
为 


(B) 


y (I)=ay(r)— f(r) 
而 (2z) 及 jz) 均 连续 ， 所 以 yz 连续 。 击 此 可 见 ， 车 f(x) 有 
4 阶 连续 导数 ， 则 
yr)=ay (IT) f(r) 
所 以 yr) 连续 。 又 
yx)=ay tr) (rz) 
可 见 y (7Y) 世 连续 。 这 就 是 说 ， 攻 (8) 中 方程 看 端的 帅 数 {7) 
有 4 阶 连 续 导 数 ( 当 ZzE[o, 四 )， 则 问题 (如) 有 浴 一 解 y(x)， 它 
在 Ca ,b J 有 6 阶 连 续 导 数 。 
现在 ， 设 jz) 确 实 有 具有 4 阶 连 续 导 数 ( 当 TEfa,b5]) ， 我 
们 来 证 明 (75) 式 ， 即 
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yATI=y rt (rt hE(r;h) 
=y(r)+t + F(z) 
《为 便于 书写 ， 将 万 (zf; 有 有) 下 记 为 BE,(7X)) 
在 此 式 中 ，y,(z) 是 (上 8,) 的 解 《 即 凑 分 解 )》 ， yl(x) 是 (28) 的 解 
《 凤 宇 确 解 )， 上 &,(z) 有 界 。 
用 “待定 函数 法 ”， 先 设 (75) 式 成 立 ， 县 设 Y(z) 在 [a, 8] 
上 有 4 阶 连续 胖 数 。 和 将 (75) 式 代入 (B83,), . 即 代 入 方程 
— Lys rn) — 2y .rT + Yr) 天 
Tuylz) = f(r,) 
记 
Pylri)=[y (xT) — 2 CT TF YT RR 
于 是 
— Py — RDN) PE (rT) + ay(t;) = f(r;) 
E (78) 
但 yz) 行 6 阶 连续 导数 ， 所 以 


晤 E31 
YTy#1)=YL;)+ hy’ (Cz)+iy? (zx;) 十 Sy" (Xi) 


十 4 (XT ) + < Vr) + oh) 
4 51 
h? hs 
Yr)_1)=y(r;) hy (tT {zr;) = arY (3) (zx;) 


hs 二 hs 5’ 6B 

ty (TY (T+ OC ) 
4! 5! 

Py(lr)=[Ly(r ,~ 2 ) 4 yr hh 


一 -二 [ hy*(x;) 十 Tah YO)+ 0Chs) ] 


hn 


I 


y* (7x) + hy 2s) + Oh) 
又 f(z ) 有 4 阶 连续 导数 ， 所 以 
Wr Wr +h (rN) + ) 
1 


+ 了 ha Cx) + oCht) 


WW (rj {rr,) hI (Cx;)+ sh (x)) 
I 


一 二 hai ea (Cz)) + On) 
31! 
于 是 
PWCOrD)=CW (rs ) -2 (zr) + Wr /8 
= 站 [AC5) + 0(h4) 
= (rx,)+ 0(h’) 
所 以 (78》 化 为 
CD) YN EtOH?) 4 0(h:)) 
— PEL)+ ay(r;)+ a (r+ aah 
“Elr,)=/(r)) 
时 
zf toy ed +t i -Sy (zh 
Wr(ri)+ aW (xr,) ]+ COC1) — PE,(z),) 
+ak(r) = 00 (C78)’ 
由 此 可 见 ， 若 矿 (x) 满 足 
Vy 1 
CCY [ -Ww (x)+ aW(r) = Ty (x) 


W(0)=0 
羡 太 (zz) 有 4 阶 连续 导数 ， 又 
DY ~PE(z)+aB(r)=0(1) (j=1,2 ,1) 
E(0)= 0 =E,(1)® 


志 “4(1 表示 有 界 县 

加 80z = YT) + iH ry +htE,trl, By Nn) = yy (00), y= yr(1) 
(0) + 用 忆 ,(0) = 0 对 任何 有 = 二 成 立 {n=2,3,%) =0 
“上 (0) = 0 向 理 ，WWV(1) = 0， E,()=0 
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则 (78) “ 式 成 立 。 在 这 里 ， 注 登 

—y" (zx)+oy(r)= f(r) 
所 以 

一 5 (T+ ay(ri)= fx) 
自然 玻 立 。 


观察 方程 《C)》 ， 这 是 常 系数 级 性 常 微分 方程 。 方 程 右 端的 
已 知 丽 数 j5y (9 (z) 有 2 阶 连续 导数 。 所 以 它 的 解 洲 Cz) 存在, 且 
由 方程 (C) 可 知 刻 (z) 有 4 阶 连续 导 续 。 

再 观察 方程 《D》， 和 看 它 是 否 有 有 界 解 玉 ,(x;)。 把 CD) 详 细 
写 出 来 就 是 


—[CE(r;r) —2E,Cx;)+ Ex /a+ aF tri)=0(1) 


即 
Er) t+ 2B rT) — Er tanbE, (ri)=001)h? 
=0(h?) 
令 
E(x;)=E, 
则 
-Ft+2E,—- E+anrk,=0(h) (D’) 
现在 考察 如 下 方程 
人 《Ey 
£,= 0=E, Cf 1,2 一 
因为 
(Biri—2B,+ EN)E,t+tohE}=9E, 
又 
a Bl 
SDE-2E+ENE, = OE -BE-EtE,.)E, 
i=-l f=l 


m=l mn-1 


= XE-E)E,- EIB 
j=1 


jel 
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全 一旦 一 


= (EE, NE,,, 3 DE-E)E, 


f=n | 
和 一 人 
DDE,~ENNCE EN)+(E, -EVE, 
j=i 
(Ei—E.)E,.. 
由 一 多 
=— SE-En) +t(E,- E)(E,~E,) 


1-! 


(Ei~E)(E,_~E,) 


由 


=— SE,-En) 《注意 Eo= 0 = EE,) 

ji 一 

所 以 
rm 一 人 章 一 下 nal 
PE-En) 和 DRE}= 9,8, 
j=40 ;一 I 1 一 ! 

于 是 
n=-l 了 一 了 对 一 1 
SE En): 十 DnE}= SpE 
了 一 自 fm0 fen 

又 


J 写 已 )| 


1 一 人 
R11 ns—l 
< 2 VR TE Ey 
fe 从 De 


全 这 里 用 到 著名 的 Catcefhy 不 等 丈 卫 吕方， <Y a? Ys 


i==1 1 


所 以 
sl 3 一 1 | 
I sn (有 人 "( 用 仙 内 Sw,5,) 
bn '=0 i=0 
于 是 
"1 人 
~ -一 1 2 rn 
| 二 | i 0 (> og 人 >) 


nV nmax lp nmax Ej 
车 i 


(注意 nh 二 1) 


可 见 
max lIE ,| <nmax [9,! :max [E | 
i 1 i 
所 以 
(max | 五 站)2<<nzmax |P; .max |E sl 
了 i f 
由 此 得 
m 五 < 扫 1 
RE Ej max | 人 中 (79) 
0% jn-i 


《注意 nh=1》 

《79) 式 是 关于 方程 组 〈 五 》 的 解 的 一 个 先 验 估 计 式 。 至 于 
方程 组 《五 ) 的 解 存 在 且 唯 一 ， 只 要 观察 〈 巨 ) 的 系数 矩阵 就 可 
以 肯定 。 

在 方程 (E》 中 ， 取 %, 二 0h 有)， 就 是 方程 (站 )。 由 (789) 
式 ， 可 知 (已 站 的 唯一 解 达 并 了 = 1 ，2 ,"…,n 一 1) 满 足 


hk 
max IE < 
0 号 :号 "一 1 


+ 13Y* 
之 0) 时 ， 满 足 
max i&; <M 


US -1l 


好 卢 ;= 二 世 ;(z;) 有 界 。 于 是 (75》 起 得 证 。 


8 3 外 推 法 在 偏 微分 方程 数值 解 中 的 应 用 


在 上 一 节 中 介绍 了 常 微分 方程 数值 解 的 外 推算 法 ， 这 方面 的 
工作 所 经 相当 成 熟 ， 有 了 许多 好 结果 ， 在 国际 上 也 早 就 为 计算 数 
学 家 熟知 了 。 人 但是， 如何 骨 外 推 法 求 求 偏 微 分 方程 的 数值 解 呢 ? 
在 国内 的 教科 书 中 还 没有 谈 到 过 这 个 问题 。 国 际 上 已 有 一 批 成 
困 ， 目 前 中 、 外 数学 家 还 在 紧张 地 研究 车 售 微分 方程 中 的 外 推 
法 。 这 个 领域 中 的 外 推算 法 ， 在 理论 上 难度 较 大 ， 在 实际 计算 方 
面孔 比较 复杂 。 在 本 节 中 ， 只 对 一 个 具体 的 偏 微分 方程 进行 外 扒 
计算 ， 并且 结 合 这 个 具体 问题 介绍 一 点 现代 结果 ， 使 读者 初步 了 
解 偏 微分 方程 外 推算 法 的 一 些 特点 。 

对 于 

一 Au(r,y)= 二 了 f(z,y) (zy)}E (和 0 为 平面 上 的 开 姑 
(4) 1 形 城 ) 

.#| so 一 0 
在 此 ，<1x6 一 5 十 zx 是 二 维 Laplace 算 符 )， 符 苇 9 只 表示 区 
域 包 的 边界 。f(z,y) 是 在 区 域名 中 有 定义 的 已 知 疯 数 ,问题 (4) 
就 是 ， 求 函数 &(z,y)， 要 它 具 有 二 阶 连 续 仿 导数 ， 且 在 矩 形 医 
域 2 内 满足 方程 一 人 4 二 了 (x,y)， 在 矩形 的 边界 3 人 2 上 满足 边 值 
条 件 x(z,y]jso= 0 《当然 要 求 x(z;g) 在 闭 插 形 只 = U98 上 
连续 ) 。 方 程 

—Au(r,y)=f(r,y) 

是 大 家 熟知 的 Poisson 方 程 . 间 题 (4) 称 为 Poisson 方 程 的 第 一 边 
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利 问 题 。 我 们 用 涛 分 法 求 央 〈44) 的 数值 解 ， 并 进行 外 推 处理 。 
最 后 说 明 外 推算 法 的 理论 人 根据。 
现在 ， 将 问题 (A4) 具体 化 为 
(CA) A 当 (T,y)E LN 
u| ao 一 0 妇 : 0 所 2 所 Z (现在 42 是 正方 
到 也 
OA 形 ) 。 


第 一 步 ” 写 出 与 (4) 相 应 的 差分 方程 ， 并 求 出 此 差分 方程 的 
解 。 


将 [0 ,z] 区 间 分 成 ”个 相等 的 小 区 间 ， 分 点 是 
To= 0, X= h, ry =2h,: ;Ts 二 
又 将 [一 所， 区间 也 +n 等 分 ， 分 点 是 


x 
Wi 小 = 一 本 二 hh Y= -2h 


# ,二 一 于 十 Pe A ( 2 


现在 正方 形 吕 被 分 成 姑 个 小 正方 形 格 子 ， 它 们 的 顶点 《又 称 “ 格 
点 2 或 “节点 59 分 唱 是 

(ria) jy 上 二 0,1 ,2 ,#8 《共有 (na 十 1 个 节点 ) 
记 


Uri YI) = ty 
因为 篇 导数 有 近似 式 
HT Va He Un) fh uri_t, ys) 
: Ui) hh 
又 二 阶 和 何人 导数 有 近似 式 
WT VLU TY — BT Yh 
[uin A | 


= Ljrs — 2875+ Wj/h’ 


局 理 ， 有 
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Uy Tp YH 
所 以 
A TY) Hy Tj Ys) 
Ent Het Hj He 1 A A? 
(80) 
在 (80》 式 右 端 出 现 的 是 在 节点 (zj,y) 及 其 上 上 、 下 、 左 、 右 由 
车 点 处 的 函数 值 。 这 个 式 子 ， 常 称 “五 点 差分 格式 ”。 记 
[jroat Wing dparit in-1— A445 Ou 
把 今 称 为 一 差分 算 子 。 与 《A4) 问题 相应 的 差分 方程 可 以 写成 
OR (j,R=1,2,,7n—1) 
Vi 一 0 《 当 7 ,二 数 中 有 一 数 是 0 或 4#) 


现在 米 解 方程 (41)。 先 取 1 = 2， 则 4 地。 这 时 


CA,) 


一 一 四 一 
立 0 一 0， TT? Ts 二 丰 


二 一 0 二 六 
yo 二 2 Vi 二 40， ys 一 也 
又 
OU t Ht Hot His ~ Au 
方程 (A1》 和 化 为 
TN ， 
Gut ttt tot His—4un =(— 2)(3) sin cos0 
fo 一 上 0 一 He 一 8 一 Ht 
实际 上 只 有 2 是 未 知 量 。 由 此 方程 解 出 xls 得 


x? 
二 可 二 1.2337 


方程 《4) 的 精确 解 是 


zt) 一 Sinzcoos2y (这 可 以 直 控 验证 》 
上 由 此 得 


充 .NX 
mt. = Sin 一 一 
(5, 0) = sm oo- 
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可 见 用 #0 )= wa 一 1.2337 作 为 x( 二 ，0 )= 1 的 近似 信 ， 
误差 太 大 了 。 人 怎么 办 呢 ? 
我 们 很 容易 想到 一 个 办 法 ， 就 是 在 求 出 xi 人 亚 ， 0 ) 之 后 ， 再 


求 出 (了 ，0 )， 然 后 试用 公式 
总 4 -= 人， 
进行 外 推 。 也 许 这 样 可 以 得 到 更 好 的 近似 值 。 
虽然 我 们 现在 还 没有 把 握 ， 但 不 妨 一 试 。 反 正 问题 ( -4 的 精 


确 解 已 经 知道 ， 可 以 直接 看 出 这 样 算出 的 xx( 革 ，0 ) 的 误差 ， 比 
ui 人 ( 瑟 ，0 ) 及 ws( 于，0 ) 的 误差 究竟 是 大 些 还 是 小 些 ? 这 是 一 
个 有 趣 的 试验 。 数 学 定理 往往 是 通过 大 量 观 察 才 发 现 的 。 现 在 ， 


我 们 来 束 xi( 于 ， 9 )。 为 此 ， 写 出 步 长 为 6/ = 志 的 差分 方 各 
( 4) ,也 就 是 


《-41 Pe HR THU gt Hr duir=t~—2). 
l i (j,k=1,2,3) 
Wx 二 0( 当 7 ,不 二 数 之 一 为 0 或 4 时 ) {j,K =1,2,3) 
这 相当 于 将 正方 形 人 2 等 分 为 入 二 16 个 小 正方 形 ( 取 %n = 二 4) 。 可 
以 按 下 式 和 迭代 求解 。 

Cm) | tm—1) tm-1 tm—1) tm i 

| Ur tt Hit ug VR 

nz? /fT ， ff x 立 
+ Sin (7)-cor( -二 + 天 ] 


ufrs=0 Ee 0 i (fj,K=1, 
2 ,3 


初 值 /0% 可 任 取 (j, 太 =1,2,3) 
当 | 一 uf?x ?| 所 e(j ,= 1 ,2 ,3)， 即 停止 送 代 运算 ， 以 
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xz 作为 nx 的 近似 值 。 钙 
现在 ， 取 ee=0.0001， 得 到 (41) 的 解 如 下 ， 
ui =0.5265 ui,s=0.7446 Ws=0,5265 
W211=0.7446 Hz,2=1.0530 ts=0.7446 
Ws —0.5265 te. —=0.7446 Was =0.5265 


不 
“(F 5 
现在 采取 第 二 个 步骤 ， 试 用 公 
3( "(T°)"( 


在 点 ( 态 ，0 ) 处 进行 外 推 ， 得 


于 是 
,~1.0530 


5 -全 


*{ 并 一 工 一 大 
"*(5, 0) SC4x1.0580—1.2337]=0.9928 


与 精确 值 4 (5, 0 ) =sin 于 cos0 二 1 比较 ， 发 现 x*( 荆 ， 0 ) 比 


UH 的 5 0 ) 及 机 的 有 0 ) 的 误差 都 小 多 了 


这 种 算法 有 什么 理论 根据 呢 ? 许多 数学 家 探讨 了 这 个 问题 。 
苏联 学 者 证 明了 如 下 重要 结论 : 
如 果 问 题 CA》 中 方程 右 端 的 函数 f(z,y) 满 足 
f(r,y) Ec (QR) (0<a<1) 
昌 在 的 四 个 角 点 处， 有 
jz 一 0 和 六 .一 fo 一 0 
则 问题 (4) 有 唯一 解 u(x, y)Ecs**(D), 且 方 程 组 (A1) 的 解 
uz,y) 有 渐 近 展开 式 
urT = V(r y+ Er,y) 
在 此 式 中 
@@” 当 炒 知 量 个 数 增多 时 ， 这 种 选 代 过 程 收 全 得 越 来 越 慢 。 记 以 ， 当 未 知 量 个 数 
太 多 时 ， 不 要 用 迁 代 法 。 可 在 电 子 计算 机 上 用 高 斯 (Ga uss) 消 元 法 解 方程 组 ( 4， 
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IPPCz,g) 委 cpate， 了 Jr(zy)Ecsto(O) 
这 是 用 简练 的 现代 数学 语言 表述 的 结果 。 在 这 里 作 一 解释 。 
只 是 矩形 开 区 域 如 加 上 它 的 边界 3 人 2。 所 谓 
A 
是 指 f(z,y) 及 它 的 所 有 一 阶 ,二 阶 和 三 阶 偏 导 数 均 在 上述 续 D， 
且 fiz,5) 的 任何 三 阶 偏 导数 PCz yy)CPCr 是 ose frr fev， 
fr 中 的 任何 一 个 ) 满足 
(PCT yi Pr yo) ECMO (ri— ro) + (yi — ya))’ 
《zy 及 (zya) 是 台中 的 任意 两 点 ，M 是 某 个 正 数 。 
u(x,y)E cr (0) 


的 含义 可 以 类 推 。 
这 些 条 件 看 起 来 比较 难 掌握 。 但 是 ， 我 们 指出 ， 只 要 
f(z,y) Ec DB) 
《 指 f(z,y) 在 上 有 直到 4 阶 的 连续 偏 导数 ) 风 
f(r) Ec (NR) 
所 以 ， 我 们 可 以 写 出 一 个 比较 实用 的 结论 ， 即 
若 函 数 (+,y) 满 足 
“01) flry) Ect(g) 
(2 ) 在 矩形 马 的 四 个 和 角 点 处 ，f(z,y)=0 
(3 ) 在 矩形 的 四 个 第 点 处 ，f -一 上 ,= 0 
由 
(1) 问题 (4) 有 唯一 解 u(x,y)€ c'**(@) 
《2) 间 题 (4,) 的 解 x,(z,y) 有 渐 近 展开 式 
UT) =u rT yt hr y) +OChst®) 《81) 
在 此 式 中 的 这 (zx,y) 满 足 
; W(xr,y) Ec*°(0) 
这 个 推论 就 是 外 推 的 理论 根据 。 事 实 上， 由 (81) 式 ， 得 
全 < 亿 导 数 在 闭 域 可 上 连 线 ” 这 一 概念 的 严 赂 定义 在 这 里 不 讲 了 。 读 者 如 对 此 


有 兴 夏 ， 请 讽 Gilbaeg 及 Tradinger 所 著 《 二 阶 杰 科 型 俩 微 分 方程 》 第 一 
章 的 最 后 部 分 。 或 啊 吴 新 谋 所 蓝 《 数学 物 埋 方 程 讲义 》 第 -- 意 3。 
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hy 四 
uaz u(r, + 【全 ) Wr,y) + ohsr’) 


所 以 
u(r, y= hs zs) uz )]=uCr,y) + oh") 
3 3 


总 结 一 下 ， 就 是 ， 车 问题 (4) 中 的 已 知 函 数 f(z,y) 潢 足 上 
述 条 件 (1),(2) 和 (3), 则 对 于 ( 41) 的 差分 解 4,(z,y) 和 waC， y)” 


可 用 公式 se(z， 乡 ) 一 于 Rd4HA(Z DB) 一 2 人 3) (82) 


进行 外 推 。 由 下 所 得 u*(z,y) 满 足 
lav(r,y)— u(r,y)) =0Ch+") 
在 刚才 的 例子 中 ，f(z,y)= 二 2sinxcosy， 它 有 任意 阶 连 续 偏 


导数 。 又 ff: 及 f ,在 四 角 点 ( 0， -5). (0, 二 (= - 司 


和 ( =， 艺 ) 处 均 为 零 。 可 见 /z,2 ?满足 条 件 (1)、( 2 ) 和 (3 )， 
所 以 可 用 《82) 式 作 外 推 。 
刚才 只 在 点 (也 ，10 ) 处 进行 外 推 。 在 其 余 点 处 因 缺 少 xx(z， 


y) 的 值 ， 所 以 无 法 用 (82》〉 式 作 外 推 。 但 是 ， 在 这 些 点 上 ,也 可 
用 陈 一 林 公式 进行 外 推 。(77》 式 是 关于 一 元 函数 的 陈 - 林 公式 ， 
对 于 二 元 函数 有 同样 公式 ， 在 此 就 不 详细 说 了 。 

我 们 还 要 指出 一 点 ， 车 (4) 中 函数 1(z, y) 不 满足 上 述 条 件 
(1), (2 ) 和 (3 )， 还 要 用 (82) 式 进 行 外 推 ， 可 能 越 推 误差 越 
大 ! 所 以 我 们 在 解决 实际 问题 时 ， 决 不 能 不 管 条 件 性 意 外 推 。 

我 们 叙述 了 吞 形 域 上 的 Poisson 方 程 的 外 推算 法 ， 实际 上 讲 
的 还 是 这 类 问题 中 的 一 个 转 例 ， 即 正方 形 域 上 的 零 边 值 问题 。 对 
于 更 复杂 的 区 域 ， 它 的 外 推 法 也 更 复杂 ， 而 且 还 在 不 断 改 进 。 关 
于 “ 偏 微 外 推 ” 这 个 园地 ， 现 在 好 像 还 挂 着 “施工 重地 、 闲 人 免 
进 的 牌子 。 苏 联 数学 家 马尔 丘 克 和 谢 多 洛 夫 所 著 的 《 差分 及 其 外 
推 》 一 书 的 英文 本 83 年 在 西方 出 版 ， 在 此 书 中 对 “ 偏 微 外 推 ” 给 
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出 了 一 些 很 好 的 结果 ， 影 响 颇 大 。 但 书 中 对 有 的 关键 问题 的 证 明 
未 作 说 明 。 总 之 ，“ 偏 微 外 推 ” 工 作 还 远 远 不 如 “ 常 微 外 推 ” 工 
作 做 得 完善 。 希 望 有 兴趣 的 读者 能 在 这 个 领域 中 进行 工作 。 


$4 万 贝 格 (Romberg) 算 法 


讲 外 推 法 而 不 提 龙 贝 格 算法 是 不 合适 的 。 

Romberg 在 1955 年 提出 一 个 很 好 的 外 推算 法 。 用 这 个 算法 可 
以 逐次 消去 从 低 阶 到 高 阶 的 各 误差 项 ， 写 出 一 个 由 近似 值 组 成 的 
“外 推 表 ”。 在 这 个 算法 中 , 先 算出 一 些 初步 的 近似 值 , 然 后 求 出 
第 一 次 外 推 值 。 将 这 些 值 再 外 推 一 次 ， 得 到 第 二 次 外 推 值 。 再 将 
这 些 值 外 推 一 次 ， 得 到 第 三 次 外 推 值 。 如 此 继续 下 去 ……。 一 般 
说 来 ， 如 果 先 求 出 # 个 初 近似 值 ， 可 以 得 到 n 一 1 个 第 一 次 外 推 
值 ， 再 得 到 nw 一 2 个 第 二 次 外 推 值 ，…， 最 后 得 到 一 个 x 一 1 次 外 推 
值 。 最 后 得 到 的 这 个 近似 值 ， 它 的 误差 会 相当 小 了 。 

为 了 说 明 这 种 算法 ， 还 是 回 到 求 z 的 近似 值 这 个 问题 上 来 。 

我 们 打算 作 4 次 外 推 ,为 此 ， 先 要 求 出 5 个 初 近似 值 ,还 是 用 


i 
d, =nSin— -7 
1 


分 别 取 # 二 4，6，8，12，24， 得 


, i 
det =4SIn =2.8284271 


dr =6sin 志 = 3 


d £2 =8sinT-—3.0614674 


7X 
df = 12sin 15 =3.105826 


df* = 24sin 21 =3.13262856 
第 一 次 外 推 值 是 


(3Y * 139* 
1 /IViX3 -7I /LV 
(3 (3 (3) 人 
x 2.828427] =3.1372512 


号 ww 全 - 


x 3.0614674 (01 , ( 


GG) 0) 9) -ls 


| (s) (a) 
df = 8/ .x3,1058286 — ‘1% 

(0 Js 
8 12 8 2 

x3.064674=3.1413171 


加 3) 


t J x 3.1326285 一 Gy- 0 


xx 和 .1058286 一 3.1415618 


第 二 次 外 推 信 串 
ly pe 
| dt =. (z) _ x3,1404968 - 3 


人 -可 


x3.1404968 二 3.1415905 


1 1 
| de i eae ( 动 ee 


| “0 G7 


x3,1413171 一 3.1415923 
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第 三 次 外 推介 是 
人 
| go) S05 i 
1 T 2 1\: /1 
| 0 CG. (2 (3) 
x 3.1415786=3.1415919 
1 2 1 2 
(6) (x 
el RE 
© -® TF -® 
x 3.1415905=3.1415923 
第 四 次 外 推 值 是 


1 > 
动 


Are 二 = 


(2) 


x3.1415923 一 pe 
1 2 1} Nz2 1 2 1 汪 

(4) -a (3) -a 
x3.1415919=3.1415922 

最 后 这 个 近似 值 4 1”， 已 经 相当 精确 了 。 由 此 得 到 关于 x 的 近 


似 值 的 外 推 表 如 下 
者 5 


>3H415919 


D3 e15905< | D71415922 
, 失忆 1 15 923 
PRALILEER: | | 


如 果 您 认为 4r” 还 不 够 精确 ， 那 就 多 算 一 些 初 近似 值 ， 青 


(DD 如 时 不 用 外 推 法 ， 只 计算 站 径 为 1 的 贺 内 接 正 多 边 形 局 -长 。 紧 算 到 区 内 接 正 
2800 边 形 ， 它 的 轴 长 址 是 3.1415918。 由 此 可 见 外 推 拔 的 作用 。 


让 
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多 外 推 几 次 。 用 电子 计算 机 进行 外 推 ， 是 很 方便 的 ， 这 种 程序 很 
简单 。 不 过 ， 外 推 次 数 过 多 也 不 好 ， 因 为 运算 过 程 中 的 舍 入 误差 
将 成 为 一 大 障碍 ， 使 最 后 的 外 推 值 的 误差 超过 理论 上 的 估计 ( 当 
然 这 种 情形 是 可 以 控制 的 ， 我 们 不 去 讨论 它 了 ) 。 

现在 来 说 明 上 述 外 推算 法 的 理论 基础 。 


令 二 一 ni, lid, 一 nsin -=d(h), 于 是 


d(hy=nsin®_ =.1.8in hn 
好 到 


Ens 了 过 员 六 3 十 vy “一 … ] 
| Th 了 { 区 产 3) 可 (xh) 


3 5 
| 和 | A 
31 51 


即 
d(h)=r+ eht+ cht+ chet ec hs ... (83) 
这 是 df(8) 的 渐 近 展开 式 ， 它 有 无 穷 多 个 含 天 的 偶数 次 者 的 项 。 
《83》 式 是 我 们 进行 多 次 外 推 的 理论 根据 。 
车? 分 别 耶 正 整 数 piyzrsy ny) 
则 相应 地 有 A ,= 元- 的 序列 
hi Rs, 下 下 
于 是 ， 有 近似 值 序 列 
C( 太 qd( 太 GCC 下 
现在 设法 作 一 次 外 推 。 因 为 
dh )=rt+ chit+ cht+ ch it chit 
d(h ri)=n+ ch tt coht st cash t+ chs, te 
所 以 
LAE 1 十 Ca 正太 3 十 Cs 声 * 正 3。 
十 chap? it 
hidCh rd)—=ahit ehh t+ crh?ihtrt cah2ha, 
+ ch? ht 


于 是 
hid(Ch rn)— hi dh )=n( hh?,,)+ cs hht,, 
hI col heh — heh + ch hs hh?,.) 
Fe 

由 此 得 到 
hd(h,)— hiadch) 

六 —h?,, 


th — cih? Chit hi)—. (84) 


Te hin, — cshih?, (hh? 


现在 令 

d(h,)=d{" 

Chridlh ra) hs dh N(R he) = 
其 

df ot hid mh? de (hi— hh?,) 《85) 
内 (83) 式 可 见 


de =rt echt ehit cshst+ chs tue 


所 以 
ds?” =Ax+0(h?) (86) 
由 (84) 及 (85) 二 式 可 兄 
di? =xw+0(hih?,,) (87) 


这 就 是 说 ，df” 比 dq” 的 误差 小 得 多 !1d:n” 就 是 第 一 次 外 推 
什 。 

我 们 再 设法 求 出 第 二 次 外 推 值 45””( i = 0 ,1 ,2 ), 使 它们 
比 &1? 的 近似 程度 更 高 。 


把 (84) 式 写成 
dV nachth?s ~ csh?h? (hi+ hry) 
— Chih (ht hi) Oo (88) 
于 是 
df =n ch hs — Cahir hi hi +t h2,,) 
— Cah ho hi he) — (C89) 


为 消去 此 二 式 右 端的 低 次 项 ， 以 有 3 匀 (89) 的 各 项 ， 以 h2,, 科 
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《88) 的 各 项 ， 再 相 减 ， 得 
hd hi dr ={hi— hi) 
— cshih?s ,hi hi, mh?:) 
Ohhinh (hh 
— Chi hah 
所 以 
[ChidE Dh df I hi hs) r+ ch hh?,, 
eg 
令 
(一 
则 人 和 一 下 十 chehts hiyat 《90) 
于 是 
人 一) (91) 
qf 即 第 二 次 外 推 人 入。 比较 (91) 和 “(87》 两 式 ， 了 可知 d:? 的 
近似 程度 更 高 。 - 
再 求 第 三 次 外 推 值 。 由 〈90》 式 ， 可 知 
人 一 二 Cs 
让 


必 


A C92) 
dds” 是 第 三 次 外 推 值 。 
一 般 说 来 ， 令 
本] 一 hrk) (93} 
可 以 证 明 
dA£? =xw+ Oh? hi ,hi ,yr) (94) 
《93)》 式 告诉 我 们 如 何 从 第 一 1 次 外 推 值 4 2, 得 到 第 KK 次 
外 推 值 4 £'?。 
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把 各 次 外 推 值 排 成 碟 。 就 得 到 外 推 表 。 

关于 求 z 的 多 次 外 捧 法 攻关 如 沦 基站 。 己 经 讲 完 了 。 这 种 多 
次 外 推 筑 法 ， 就 是 龙 贝 糙 算法 。 这 姑 算 法 已 经 用 刘 许 允 方 面 。 

党 来 也 是 滩 巧 ， 我 们 讲 外 推 法 ， 既 是 从 求证 井 始 ， 又 是 以 
求 T+ 值 告 终 。 中 间 涉 及 背 干 与 微分 方程 有 闫 的 问题 。 伺 是 外 推 法 
的 活动 范围 决 不 止 紫 。 它 存 积 分 方程 、 代 数 方 程 ……， 各 方 曾 ， 
有 广泛 的 运用 。 我 们 在 第 四 章 中 讲 到 的 有 限 元 法 ， 也 是 一 种 近似 
解法 。 所 以 用 有 限 元 法 得 到 的 解 ， 岂 可 以 外 推 。 次 80 年 以 来 ， 中 
国学 者 在 有 限 元 解 的 外 推 方面 ， 让 得 了 巨大 进展 必 。 本 书 因 篇 申 
所 限 ， 不 能 淡 这 些 工作 了 。 

大 们 看 到 ， 电 子 计 算 栅 的 迅速 发 展 ， 并 没有 降低 对 计算 方法 
的 要 求 ， 反 而 刺激 了 计算 方法 的 发 展 。 广 泛 运 用 和 发 展 外 推 法 ， 
是 当前 的 一 项 迫切 的 任务 。 
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